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ABSTRACT 

Experimental and theoretical ion and electron temperatures in the top- 
side ionosphere are investigated. Experimental results come from an analysis 
i t ,ti coherent -scatter data taken at Arecibo, Puerto Rico. Consideration of 
the energy balance equations gives the theoretical ion and electron tempera- 
tures. 

Two ions, 0* and H + , are assumed to be the only ions present in the 
topside ionosphere. Experimental evidence shows that a third ion, He* is 
relatively unimportant compared to 0* and H*. Photoionization of neutral 
particles, charge transfer, and diffusion of these charged particles are 
the important physical processes which are used to calculate their number 
density profiles. 0* and H + number densities are shown to affect the 0 + , 

H + , and electron temperatures. 

Analysis of theoretical autocorrelation functions of the incoherent- 
scatter power spectra shows that number densities and temperatures of 0 + , 

H*, and electrons can be deduced from experimental autocorrelation functions. 
Experimental temperature results indicate that the maximum difference between 
H* and 0* temperatures occurs at 550 km. This agrees with temperatures 
calculated from the energy balance equations. The experimental and theoretical 
H* and 0* temperatures agree fairly closely from 200 to 1200 km while the 
electron temperatures agree from 400 to 1200 km. 
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1 . INTRODUCTION 

1.1 The neutral Atmosphere 

The goal of aeronomy is the knowledge of all physical parameters of the 
upper atmosphere. This implies the formation of models to describe the neu- 
tral atmosphere since most physical phenomena deal with the upper atmosphere's 
neutral constituents. Two important physical parameters which are necessary 
in any reliable atmospheric model are the number density and temperature of 
the neutral species. A number of different neutral atmospheric models exist 
depending on different conditions in the upper atmosphere . The most commonly 
used models are the ones which specify a thermopause temperature (temperature 
at the upper boundary of the earth's atmosphere). These models are based on 
observations from a number of different experiments. CIRA [1972] gives a re- 
view of the atmospheric models in the region from 110 to 2000 km which is the 
region being studied in this work. CIRA [1972] atmospheric models are used in 
this work and are the models that are commonly used today. CIRA [1972] models 
are reproduced in Appendix A. 

In the remaining sections in this chapter the ionospheric regions of the 
atmosphere of the earth will be discussed. Also, collision processes, trans- 
port of charged particles, and thermal processes of the upper atmosphere will 
be discussed. After incoherent-scatter radar is introduced, the statement of 
the problem will be given in section 1.9. 

1.2 Ionospheric Regions 

The atmosphere of the earth can be divided into several different regions. 
When studying charged particles in the earth's atmosphere, the D, E, and F 
regions are discussed. These regions are identified by layers or ledges in 
the electron concentration profiles. 

The D region contains charged particles formed by radiation whose 


2 


-19 2 

absorption cross section Is smaller than 10 cm . The D region is between 
40 and 90 km. Radiation which produces ionization in the D region includes 
solar X-rays having wavelengths less than 10 A and ultraviolet radiation of 
wavelengths greater than 1750 A. Also important is the absorption of Lyman -a 
radiation. The gases 0 2 and N 2 are ionized by solar X-rays while NO is ion- 
ized by Lyman-a. 

The E region contains charged particles formed by radiation whose ab- 

-18 J 

sorption cross section is less than 5 X 10 cm“. The E region is between 

90 and 160 km. Solar radiation in the range 100-31 A and that greater than 

800 A are important while solar Lyman 6 line at 1025 A and the C(III) line at 

977 A are important since they ionize 0^, Ionization of O is produced by 

Lyman continuum at wavelengths less than 910 A. 

The F region has charged particles produced by radiation whose absorption 

-17 2 

cross section is greater than 10 cm . The F region is above 160 km. Ra- 
diation whose wavelength is in the range 100 to 800 A ionizes N 2 , 0 2 and 0. 

The F region has two layers, the Fj and F 2 layers. The Fj layer is a layer 
controlled by ion interchange processes. The F 2 layer contains atomic mole- 
cules and charged particles; diffusion is important in determining the ion 
profiles. The Fj layer is below the F 2 layer. 

1.3 Ionisation in the Region 

In the Fj region the major ions are NO*, O*, O,* and N 2 *. NO* is the 
dominant ion at altitudes between 140 and 170 km. Above 170 km 0* becomes the 
dominant ion due to the much lower loss rate for 0*. To calculate equilibrium 
number densities for NO* and 0* one must consider all the important reactions 
that bring about formation or loss of NO* and 0*. Therefore, ion chemistry 
rate coefficients are important in determining which reactions must be used 


3 


and which reactions can be neglected. For N0 + the important sources of ioni- 
zation include charge transfer of NO with 0*( 4 S), N 2 * and 0 2 *. The most im- 
portant sink of NO* is dissociative recombination. For 0* the important 
source of ionization is photoionization while charge transfer with N 2 and 0 2 
act as the important sinks of ionization. In Figures 1.1 and 1.2 are shown 
some theoretical and experimental profiles of NO* and 0* number densities. 

Both Figures 1.1 and 1.2 show that NO* and 0* densities depend on the solar 
zenith angle. Therefore, a model of F^ ionization must take into account 
variations that depend on the time of day. 

1.4 Ionisation in the F 2 Region 

The discussion of the last section concerned the solution to the ion 
continuity equation where chemistry is the important process. The present 
section gives an introduction to the current theory of the F 2 region. 

The ionospheric Fj region is formed through the competition between ion 
transport and chemical process up past the electron density peak which is the 
beginning of the topside ionosphere. The F 2 region begins at about 225 km 
and the principal ion is 0*( 4 S). At high altitudes, around 1000 km, H* be- 
comes the principal ion. The predominant neutral gas is 0, although N 2 and 
0 2 are important in 0* loss processes. The oxygen ion density in the F 2 region 
is determined both by plasma transport and chemical reactions. In the steady 
state and in the absence of diffusion or convection the ion density would be 
expected to increase with increasing altitude since the effective ion recom- 
bination coefficient decreases rapidly with increasing altitude. But at suf- 
ficiently high altitudes ion recombination is balanced by ion transport and 
a peak in the ion density is reached. So, ion transport is responsible for 
for a decreasing density distribution with increasing altitude in the topside 
ionosphere. 
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Figure 1.2 Experimental profiles of NO , 0- , and 0 for 
daytime, sunset, and nighttime conditions 
[after P'nks and Kockavts, 1975]. 
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Atonic oxygon ions produced in the upper Fj region diffuse downward 
with typical velocities in the range 10-20 ■ sec -1 [BouMll, 1962]. A series 
of electron density profiles at Arecibo, Puerto Rico are shown in Figure 1.3. 

Hie production of hydrogen ions in the Fj region is brought about both 
through photoionization and through accidentally resonant charge exchange 
[Hanson and Ortenburger , t 1961] 

0 + ( 4 S) ♦ H( 2 S) 5 0( S P) ♦ H + ♦ tE (1.1) 

The energy defects of (1.3) are snail relative to themospheric energies 

(kT 'u .1 eV) so the reaction proceeds rapidly in both directions. Using 

rates of photoionization and charge exchange given by Banks and Kookarts 

[1973], one finds that photoionization is important oily when the oxygen ion 

2 3-3 

density is less than 10 - 10 cm which is far less than -region oxygen 
ion densities. Calculations of hydrogen ion densities for chemical equili- 
brium (rate of increase in ions per unit volume per unit tine ■ zero) show 
that the hydrogen-to-oxygen ion density ratio increases rapidly with altitude 
because atomic oxygen decreases nore rapidly with altitude than atomic hydro- 
gen does [Banks and Kookarts t 1973]. Thus, hydrogen ions eventually become 
the dominant ion species in the topside ionosphere. 
l.S Collision Processes in the Topside Ionosphere 

The study of collisions is important in the study of electron and ion ther- 
mal balance. Ions, neutral particles, and electrons exchange energy through 
elastic, inelastic, and Coulomb collisions. In this section, collision pro- 
cesses of binary mixtures in thermal disequilibrium will be considered. 

A gas of particles of type i is assumed to have a Maxwellian velocity 
distribution given by 

f(v.) • n.(m./2*kJ\) 3/2 exp(-m.v . 2 /2kT.) 

t* % % > % V % 


( 1 . 2 ) 



t profiles of n(0 + ) and n(l 
over Arecibo, Puerto Rico 
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where is the velocity of the i-th species, 
is the number density, 
m. is the aass, 

is the temperature, and 
k is Boltzmann's constant. 
f(v.) is normalized so that 

X* 

► 

m n . (1.3) 

One of the most important concepts of kinetic theory used in aeronomy 
is the concept of the average collision frequency between particle type 1 
and particle type 2 given by 




ff o( J g)fiv 1 )f(v 2 )dv 1 dv 2 

* 


where g is the relative speed between particles of type 1 and 2 


(1.4) 


g ■ Vj - t> 2 0*5) 

0(9) is the velocity dependent collision cross section which depends on the 
nature of particle interaction, i.e. Coulomb force, induced dipole attraction, 
rotational and vibrational excitation, etc. 

For example, the collision n equency of a particle with atomic oxygen 
with the hard sphere approximation of ^ ■ 5 x 10"^ cra^, m • 2.8 x 10"^ gm, 
and T ■ 300K is: 

v iQ ■ 1.1 x 10” 10 n(0) sec“* (1.6) 

Other types of collision processes will be discussed in more detail later. 
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1.6 Transport of Charged Particles 

The ion composition is an important parameter in the study of ion con- 
duction in the topside ionosphere. The ion composition is controlled by 
diffusion in the topside ionosphere. This diffusion is not ambipolar diffu- 
sion because in the topside ionosphere several ion species are involved and 
the scattering cross section for ambipolar diffusion is veiy small compared 
to the Coulomb scattering cross section of thermal protons with oxygen ions; 
[Hanson and Ortenburger , 1961]. 

The ion composition and ion number density profiles vary from night to 
day. Typical profiles of ion composition are shown in Figures 1.3 and 1.4 for 
nighttime and daytime data. The daytime data shows the case of outward plas- 
ma flow while the nighttime data shows the case of inward plasma flow. The 
maintenance of the large nighttime number density peak is due to thir inward 
plasma flow. The transport of charged particles will be discussed in more 
detail later. 

1.7 Thermal Processes 

The transfer of heat between charged particles and neutral particles and 
other charged particles is important in the study of the topside ionosphere. 
While all particles have the same temperature in the D and E regions, in the 
F region ions, electrons, and neutrals each have separate temperatures. The 
neutral temperature increases rapidly with altitude above 120 due to the 
absorption of solar ultraviolet radiation by the neutral atmosphere. The 
increase of temperature gradually becomes small at altitudes above 500 km due 
to the decrease in absorption of solar ultraviolet radiation with altitude 
and a large neutral thermal conductivity. The ions and electrons have almost 
the same temperatures as the neutrals at altitudes below 200 km due to large 
thermal -electron heating rates by photoelectrons and large heat transfer 



ALTITUDE (km) 
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rates between ions and neutrals, and electrons and neutrals, through elastic 
and nonelastic collisions. The electron temperature increases rapidly above 
the neutral temperatures with altitude above 200 km since heat transfer to 
neutrals and ions decreases wi > altitude. The ion temperature increases 
above the neutral temperature above 200 km but it is lower than the electron 
temperature due to larger ion-neutral heat transfer rates than electron-neutral 
heat transfer rates. At high altitudes the electron and ion temperatures are 
controlled by thermal conduction. Since the thermal conduction is different 
for electrons, positive oxygen ions, and positive hydrogen ions, these parti- 
cles can all have separate temperatures. This will be discussed in more de- 
tail later. 

1.8 Incoherent-scatter Radar 

Incoherent scatter (Thomson scatter) is a process where incident electro- 
magnetic radiation is reradiated by charged particles with random phase. When 
the exploring wavelength is large compared to the Debye length given in equa- 
tion (1.7), the scattering of radiation can be thought as arising from density 
fluctuations caused by plasma waves. The principal component of the returned 
signal are scattered waves introduced by the presence of ion-acoustic waves 

and electron- induced waves at the gyrofrequency (u^ ■ e.B/cm ) and the plasma 

2 1/2 

'requency (p = (n g e q) ) (<? is the speed of light, is the permittivity 

of free space, B is the earth's magnetic field, e, n , n are the electronic 

6 0 

charge, electronic mass, and electron number density, respectively). 

Although it was first thought [Gordon, 1958) that the spectrum of re- 
flected incoherent-scatter signals would he Gaussian in shape with a center 
to half power width of .71 times the Doppler shift of an electron approaching 
the radar at the mean thermal speed, later observations [Bowie e , 1961] proved 
that the bandwidth was considerably less than the Doppler shift of the elec- 
tron. Subsequent theoretical work [Fe,jer, 1960; Dougherty and Farley, 1960; 
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Salpeter, 1960a; and Uagfove , 1961] has shown that the Influence of ions can 
effectively narrow the spectrum width when the exploring wavelength is much 
larger than the Debye length. The Debye length or Debye shielding distance 
is given by 

2 1/2 

Dp « (e^kT/n^e ) ' meters (1*7) 

where k is Boltzmann's constant, T g is the electron temperature. D g is of 
the order of 1 cm or less below altitude 1000 km and rises to about 6 cm at 
altitude 2000 km. Arecibo Ionospheric Observatory (see section 1.9) operates 
the radar at wavelength of 69.8 cm. This is the radar facility which will be 
used in this work. 

Information about the temperature and number density of electrons, 0 + , 
and H + can be deduced from incoherent -scatter returns. This is done by fit- 
ting theoretical generated spectra with experimental spectra. This will be 
discussed in more detail later. 

1.9 Statement of the Problem 

In this work both experimental observations and theoretical analysis of 
the thermal structure of the topside ionosphere will be studied. The experi- 
mental observations are in the form of incoherent-scatter returns from 
Arecibo, Puerto Rico which is located at a magnetic latitude of 30°N. The 
magnetic dip angle, I, is 50° at Arecibo. The antenna at Arecibo is a spheri 
cal disc with a radius of curvature of 245 meters and a diameter (viewed 
from above) of 305 meters. Details of the system have been described by 
Gordon [1964]. Parameters such as temperature, density, and composition of 
the charged particles can be deduced from incoherent-scatter returns by use 
of incoherent-scatter theory. Discussion of the theory and the analysis of 
the experimental data will be given later. 


13 


Theoretical analysis of the thermal structure of the topside ionosphere 
involves solving the appropriate energy balance equations for the charged 
particles in the topside ionosphere. In this work the major ions are assumed 

to be the positive charged oxygen ion, 0 + , and the positive charged hydrogen 

+ + + 
ion, H . Together the number densities of 0 and H add up to be the elec- 
tron number density or 

n = n(0 + ) + n(H + ) 

& 

where n is the electron number density, and rc(0 + ) and n(H + ) are the number 
densities of 0 + and H + . 

The electron, 0 + and H + energy balance equations are coupled nonlinear 
differential equations. Analytic solutions are impossible to obtain, so nu- 
merical methods are used to solve these equations. By using different atmos- 
pheric models theoretical profiles for the electron, 0 + and H + temperatures 
can be obtained and compared with those of corresponding atmospheric condi- 
tions. 

Experimental analysis of incoherent -scatter radar returns can be per- 
formed to obtain 0 + , H + , and electron number densities and temperatures. 

This involves numerical methods of fitting experimental autocorrelation func- 
tions to theoretical autocorrelation functions. 

The objective of this work is the investigate of thermal processes of 
electrons, 0 + and H + ions in the topside ionosphere. Through theoretical and 
experimental deduction of temperature profiles the importance of various 
heat transfer processes can be verified. Thus, comparison between experi- 
mental and theoretical results will indicate further refinements in the 
theory of heat transfer between charged particles and other charged particles 
or neutral particles. 
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2. IONS IN THE TOPSIDE IONOSPHERE 

2 . 1 Introduction 

Parameters which can easily be measured by existing experimental tech- 
niques are usually used in describing physical processes in the ionosphere. 

The observed values of these parameters are compared with the values calcu- 
lated from theoretical models. These parameters include values of electron 
temperatures, ion temperatures, ion densities, electron density, and neutral 
densities, and velocities of neutrals and charged particles. 

One of the most important parameters is the ion density distribution. 

F region ion densities have been measured or deduced from a number of experi- 
mental instruments, including incoherent scatter radars, ion mass spectro- 
meters, Langmuir probes and retarding potential analyzers aboard rockets and 
satellites, and radio wave sounders or ionosondes. Recent experiments using 
these experimental techniques [Hofjhan, 1967; Hoffman et at., 1974; Moororoft, 
1964; Knudsen and Sharp, 1966; Hagen and Hsu, 1974] have indicated that 0* 
and H* are major constituents in the topside ionosphere. In Figure 2.1 
[Hagen and Hsu, 1974], one notices that 0* is the major ion at F ^ maximum and 
H* is the major ion at altitudes above 1000 km. This is the case in all typ- 
ical profiles of ion composition in the topside ionosphere. He* is usually a 
minor ion near the altitude where the transition between 0* and H* as major 
ions takes place and is less than 10% the total ion concentration at all 
other altitudes. 

In the topside ionosphere the neutral particle density becomes small so 
ion-neutral collisions lose their importance in the ion momentum equations 
at increasing altitudes above the F ^ maximum. Therefore, the ions move along 
the magnetic field lines when there are no perpendicular electric fields 


ALTITUDE (km) 




CONCENTRATION (cm'*) 

Figure 2.1 Daytime profiles of n e , n(0 + ) , «(H + ), and «(lle + ) for the 
F region over Arecibo, Puerto Rico [after Hagen and Hsu , 
1974]. 
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present. The action of the charge separation electric field inhibits the 
movement of the ions even though ion-neutral collisions are no longer impor- 
tant. Ion-ion collisions are important since they also help in keeping the 
ions from moving quit 'y out of the earth's atmosphere. 

2 . 2 Hydrogen lone 

Oxygen and hydrogen atoms have almost the same ionization potential. The 
oxygen ionization potential is slightly greater than the hydrogen ionization 
potential. This implies that charge transfer between oxygen atoms (0) and 
hydrogen ions (H + ) can proceed almost as rapidly as charge transfer between 
hydrogen atoms (H) and oxygen ions (0*) . This reaction can be expressed as 

H + ♦ 0 £ 0* + H + LE (2.1) 

where A E is small compared to thermospheric energies in the F region. In 
Table 2.1 are given some time constants for 0 + - H charge exchange. This 
table gives data for several atmospheric models. The time constant given in 
Table 2.1 is the time constant for changes in n(H + ) when production is zero 
and when flux of H + is zero [Banke and Kookarte t 1973] . 

t(H + ) • (fCj, n(0)(8^/ir mfH)) 1 *)" 1 (2.2) 

_1 5 2 

where Q„ = (1.5 ± 0.5) x 10 cm is the charge -exchange cross section for 

E 

(0 + ,H) collisions, n(0) is atomic oxygen density, T. is the ion temperature 
and m(H) is the mass of hydrogen. 

The continuity equation governing the H + number density is 

3n(H + )/3t + V*(n(H + )?(H*)) * p(H + ) - t(H*) (2.3) 

where ? is transport velocity 

p(H + ) « I (H)n(H) «■ Q e n(H)n(0 + )(8fcr n /mn(H)) ,s 


(2.4) 
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TABLE 2.1 

Time constants' 3 for 0*-H charge exchange (after Banks and KoQkrwte, 1573). 


Atmospheric model thermospheric temperatures 

Altitude 


IKmj 

750 K 

1000 K 

1250 K 

1500 K 

1750 K 

2000 K 

200 

2.3(-l) 

2 . 3 (—1 ) 

1.9(-1) 

1 ,6(— 1) 

l.K-l) 

8.2 ( — 2) 

250 

8 . 0 (—1 ) 

6.2(-l) 

4 . 4 (—1 ) 

3 . 5 (—1 ) 

2.4 (—1) 

1.6(-1) 

300 

2.6(0) 

1.5(0) 

9.4 (-1) 

6.6(— 1) 

4 . 3 (—1 ) 

2.8(— 1) 

350 

9.1(0) 

3.6(0) 

1.9(0) 

1.2(0) 

7.3 ( — 1 ) 

4.5(-l) 

400 

2.5(1) 

8.3(0) 

3.7(0) 

2.2(0) 

1.2(0) 

7.1(-1) 

450 

7.6(1) 

1.9(1) 

7.3(0) 

3.8(0) 

1.9(0) 

1.1(0) 

500 

2.2(2) 

4.3(1) 

1.4(1) 

6.5(0) 

3.1(0) 

1.7(0) 


a Values are in sec. 


b 2 . 3 (—1 ) * 2.3 x 10' 1 . 
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is the production of H* 

*(H + ) - 1 7* n(0)»CH + ) (SWCH^J/wnCH)) 1 * (2.S) 

is the loss of H + . I(H)n(H) is production from photoionization of H. 2* n is 
the neutral gas temperature. T(H + ) is the hydrogen ion temperature. The 
e - H + recombination rate is negligible compared to the charge transfer rate. 
When the conditions for chemical equilibrium (<?(H + ) ■ 0 and 3n(H + )/3t * 0) 
are used, the H* density can be derived from the continuity equation. 


n(H*) 9 nfin 

\ " 

h 

1 , Iff) 

„Co*) 

m*) 


« fi 7»(0 + )(8fcT n /irm(H)) ,J 


Using 1(H) * (8.8 ± 3.3) x 10“ 8 sec"* and T n ■ 1000 K one finds that photo- 

ionization is not important when n(0 ) is greater than 10 cm . Since n ( 0 ) 

4 -3 

is typically greater than 10 cm in the topside ionosphere (2.6) can be 
reduced to 

^ * rsSJ- (V 3 ''"*”'’ < 2 - 7 > 

This implies the hydrogen to oxygen ion density ratio increases rapidly with 
altitude since atomic oxygen decreases more rapidly with altitude than atomic 
hydrogen does. In Table 2.2 are given some n(H + )/n(0 + ) for several atmos- 
pheric models. Equation (2.7) shows that the hydrogen ion density for chem- 
ical equilibrium does not depend on the charge exchange rate even though 
charge exchange is the dominant source of protons in the topside ionosphere. 
n (0) and >:(H) are the important parameters for determining n(H*). 

Tables 2.1 and 2.2 show that t(H + ) and «(H*)/«(0 + ) increase with altitude. 
Eventually r(H + ) and n(H*)/n(0*) become large enough so that ionic diffusion 
becomes important and controls the distribution of number density. The 
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TABLE 2.2 

Hydrogen to oxygen ion density ratios: 
chemical equilibrium (after Fanka and Kockavt8 t 1973). 


Altitude 


Thermospheric temperatures 



(km) ~ 

750 K 

1000 K 

1250 K 

1500 K 

1750 K 

2000 K 

200 

4.6 (-4) 

6.9 (—5) 

3. 3(-5) 

2. 8 (-5) 

2.8(-5) 

2. 9 (-5) 

250 

1 . 4 (—3) 

1 • 4 (—4) 

4.5H5) 

3.1 (-5) 

2 . 8 (-S ) 

2.8(-5) 

300 

4. 2 (-3) 

3 . 0 (—4 ) 

7.0(-5) 

3.8(— 5) 

3.0 (-5) 

2.8(-5) 

350 

1.2 (-3) 

6.7M) 

1.2H) 

5.1 (-5) 

3. 4 (-5) 

2.9(-5) 

400 

3. 5 (-2) 

1.4 (-3) 

2.1H) 

7.3(-5) 

4.1 (-5) 

3. 2 (-5) 

450 

8. 8 (—2) 

3.1 (-3) 

3.’H) 

i.l(-4) 

5. 3 (-5) 

3.7(— 5) 

500 

2 . 7 (—1 ) 

6.8(-3) 

6.8(4) 

1.7M) 

7. 2 (-5) 

4.4(-5) 


a 4.6(— 4) = 4.6 x 10 
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transition from chemical to diffusive equilibrium does not greatly affect 
the H* distribution because usually this occurs when «(<>*) » n(H*). When 
n(0*) » n(H*) the chemical and diffusive equilibrium solutions for the hydro 
gen ion density are nearly the same. 

2.3 Helium Iona 

The source of helium ions in the topside ionosphere is the photoioniza- 
tion of atonic helium by solar radiation with wavelength, X, less than S04 A. 
This reaction is expressed as 

J( He): H e(*S) ♦ hv(X < 504 A) + He*( 2 S) ♦ e (2.8) 

The important loss reactions involve charge transfer and dissociative 


charge transfer with 0 2 and N 2 given by 

(Y 10 ): He + ♦ 0 2 -*■ He + 0 ♦ 0* ♦ 5.85 eV (2.9) 

(Y 10 *): He* ♦ 0 2 -► He ♦ 0 2 * ♦ 12.25 eV (2.10) 

(Y 6 ) : He* ♦ N 2 He ♦ N + N* ♦ 0.28 eV (2.11) 

(Y 6 *): He* ♦ N 2 -► He ♦ N 2 * (2.12) 

See Appendix B for rate coefficients given above. 

Thus the continuity of He* number density is 

3n(He*)/3t ■ V* («(He*)?(He*)) ■ 


I (He) « (He) - n(He*)((Y 10 ♦ Y 10 *)«(0 2 ) ♦ (y 6 + Y 6 *)n(N 2 )] (2.13) 

If the conditions for chemical equilibrium are applied one gets 


w(He ) 


I(He)w(He) 


( 2 . 14 ) 


(Y 10* * (V 6 * 

This equation implies that n(He*) increases with altitude since n(N 2 ) and 
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n(0 2 ) decrease much more rapidly with altitude than n(He) does. At high al- 
titudes diffusion cannot be ignored and so therefore the n(He*) distribution 
will decrease with altitude. 

He* never becomes the predominate species in the topside ionosphere. 

Its maximum density is 5 x 10 3 cm' 3 [Hoffman, 1967]. Haneon at al. ]1969] 
showed that even at 1068 km that only 0.14% He* existed in the topside iono- 
sphere. 

2.4 (keygen lone 

Oxygen ion, 0*, is the principal ion of the F 2 region and the predominate 
neutral gas is 0. The F^ region is formed through the competition between 
ion transport and chemical process. Oxygen is ionized through solar ultra- 
violet radiation (X < 910 A). 0* is lost through slow ion-molecule reactions 
with N 2 or 0, to form NO*, N 2 *, or 0 2 *. 

The altitude profile of the 0* ion density is controlled both by ion 
transport and ionic reactions. Ion transport is important at higher altitudes 
because in the absence of ion transport the steady state 0* ion density would 
continue to rise with altitude due to the decreasing ion recombination coef- 
ficient. However, at high altitudes chemical processes cannot compete with 
ion transport. Thus, a peak in the 0* ion density is reached. 

Short wavelength solar ultraviolet radiation (X < 910 A) ionizes atomic 
oxygen at the rate 

<7(0*) - 7(0)«(0) 

where 7(0) ■ (3.2 t 1.3) x 10"^ sec -1 and n(0) is the density of atomic oxy- 
gen. The importance of the final state of 0* is in this case significant due 

to the rapid reaction of 0*(^P) with N 2 in comparison with 0*(*S) with N 2> 

♦ 4 +2 

0 ( P) and 0 ( P) can be neglected due to relatively small cross sections. 
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f 4 

Thus, only 0 ( S) reactions and transport are important . 

The important ionic reactions which 0* is lost are 

Yj * 0* ♦ N 2 NO* ♦ N 
Y 2 i 0* ♦ 0 2 -► 0 2 * ♦ 0 

where Yj ■ (1.1 t 0.3) x 10" 12 cm 3 sec" 1 and y 2 ■ (2 ± 0.2) x 10" 11 cm 3 sec’ 1 . 

The 0* ion density is extremely variable in its naxinua density and peak 
altitude. Changes in the neutral atmosphere such as temperature and density 
changes and neutral winds can affect t:»e magnitude and height of the peak. 
The neutral winds raise or lower the peak altitude of the F 2 layer if the 
drift is upward or downward. Also north to south winds tend to aove ioniza- 
tion up magnetic field lines. For increasing exospheric temperatures the 
peak density increases and the peak altitude shifts upward [Hehbeth and 
Saiwn, 19o0]. 

Outward and inward flows of 0* also affect the density profile. Shortly 
after dawn outward flow of 0* from the ? 2 region to high altitude magnetic 
field tubes replace ionization which was lost during the night as a result 
of recombination. In the evening and during the night inward flow of 0* 
(created through H* to 0* charge exchange) maintains large densities in the 
region. 

Another factor which affects the F 2 peak is the magnetic dip angle I. 

The vertical component of ion diffusion is reduced for smaller values of I. 

The effects of electric fields on F 2 region electron densities has been 
used to explain diurnal variations [ Stubbe and Chandra t 1970]. But the 
magnitude of these fields (5 mV M’ 1 ) appear to be much larger than experimen- 
tal values [Cauffman and Gumett, 1971], Thus, the role of electric fields 
in the F 2 region is generally neglected. 
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2.5 Ion Diffusive States 

The discussion in the previous sections dealt with chemical equili- 
brium states for H + and He + which exist at low altitudes in the F region. 

Above 500-700 km this is not the case. Diffusion and convection become the 
dominant processes. Mechanical forces such as gravity, electric and magnetic 
fields, pressure gradients and centrifugal forces also become important. The 
Fj region can act as either a source or sink of ionization for the topside 
plasma since diffusion strongly couples the Pj region to the topside iono- 
sphere. 

Magnetic fields are important since charged particles can move only along 
field lines if there are no electric fields. Therefore, when discussing mo- 
tions of charged particles in the topside ionosphere one speaks of motions 
along the magnetic field lines. 

In the topside ionosphere there are three diffusive states which are 
important in determining the ion density distributions. These states are 
diffusive equilibrium, inward plasma flow, and outward plasma flow (see Fig. 
2.2). The first state, diffusive equilibrium, corresponds to a state where 
there is no net transport of ions along the magnetic field tubes. This state 
rarely occurs in the earth's atmosphere at mid and high latitudes. Only at 
low latitudes does the topside ionosphere have a chance to reach a state of 
diffusive equilibrium. Still, this state serves as an important basis for 
the understanding of inward and outward plasma flow. 

Inward diffusive flow results when there is an excessive plasma pressure 
at high altitudes in the magnetic field tube. 11 and He diffuse downward to 
regions where ion-neutral reactions convert these into 0* and thereby 
-nhancing the ion density distribution in the region. Calculations of 


(a) 
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Figure 2.2 The three different diffusive states 
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inward diffusive flow shown in Fig, 2.3 indicate that the ion density distri- 
bution are unchanged from diffusive equilibrium at low altitudes (below 500 
km.) This is discussed in detail later. Outward diffusive flow of light 
ions results from low plasma pressure at high altitudes in the magnetic field 
tube. If there are no magnetospheric electric fields then the upward flow 
plasma replenishes the field tube and eventually diffusive equilibrium comes 
about. During magnetic storms magnetospheric convection and the plasmapause, 
a sharp latitudinal decrease in plasma density which normally occurs near 65° 
magnetic latitude, can move to magnetic shells as low as L ^2.5 [Ryoroft 
and Burnell, 1970). Following the storm the field tube is essentially empty 
because field tube convection moves to higher L shells [Banks et al. , 1971] 
and plasma flow must eventually fill the field tube. 

2.6 Continuity and Momentum Equations 


The basic equations for determining the ion density distribution are the 
continuity and momentum equations. The ion momentum equation for flow par- 
allel to the magnetic field line is 


9u . Bu . 

JT * u i 51T + i7 | 


l 


Bp. Z.e 
• t r 


n.m , Bs 

t* t* 


m , 




where u ^ is the ion bulk transport velocity, u ^ is the ion bulk transport 
velocities of other ions, p| | is the acceleration of gravity (positive in- 
wards) parallel to B (the magnetic field), a is a coordinate along B , A is 
the area of the field tube at point a, Z.e is the atomic charge of the ion, 
«£ is the ion density and is the ion mass, E j j is the electric field par- 
allel to B and \>^ is the momentum transfer ion diffusion collision frequen- 
cy, p. and p are the ion and electron gas pressures (p. ■ n.kT.), the sum- 
motion is over all neutral and other charged particles. 


ALT ITUDEOtml 


Outflow 
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ION CONCENT RAT I ON (cm"') 


Profiles of n(H + ) and n(0’ , ‘) for a 1000 K model 
thermosphere and different escape fluxes. The 
m(0 + ) profile applies to the largest outflow. 
Also, «(0 + ) 1 400 km * 2 x 10$ crn-o, 2^| 400 km = 
2500 K, and th^ heat flux at 400 km is 
5 * 10^ eV cm"“ s“^ [after Banks and Koakavts 3 
]• 
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The continuity equation for lone in general is 


(Bn^/Bt) ♦ V-(n^) * 


(2.16) 


where n. is the number density of the ion, p. is the production rate of ions, 

t v 

is the loss rale of ions. 

% 

The ion continuity equation parallel to the magnetic field is 


(3^/3t) + j 


1 A) 


p t ’ l i 


(2.17) 


where s is length along the field line and A is the area of the field tube 
at point 8. 

Since the ion mass is much greater than the electron mass, the electron 
momentum equation can be expressed simply as 


1 

n s’ 
e 


♦ &E i i ■ 0 


(2.18) 


Using (2,18), (2.17) can be expressed as 


3u . 3 u . 

i + u % 


Bp. Z. Bp 

‘ i . t t c 


W * u i W + *1 I + 55“ + fTmT 55” " ? v tfc (u i ■ V 


(2.19) 


where charge neutrality requires that 


n ■ E 2 .n . 
e .it 
t' 


( 2 . 20 ) 


Equations (2.17) and (2.19) are nonlinear differential equations and are 
difficult to solve analytically without making some initial approximations. 
Fortunately, the flow conditions of diffusive equilibrium and low speed flow 
permit explicit analytic results. At low latitudes these two conditions are 
enough to describe typical conditions in the topside ionosphere [Nagy et al. t 
1968]. 


■L 
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2.7 Diffusive Equilibrium 

The complex problem of solving for ion density distributions is simpli- 
fied in the case of diffusive equilibrium. When considering an ionosphere 
made up of several different ion species the assumption, • 0, for all ion 
species makes explicit results possible. Ignoring centrifugal force and 

letting u , ■ 0 equation (2.19) for steady state conditions can be expressed 

t 

as 


1 3p j 

n 3s m i 


- 2^| 


The electron momentum equation becomes 
1 3p a 

+ eEt | • 0 

% II 


( 2 . 21 ) 


(2.22) 


Charge separation electric fields in the electron momentum equation are large 
enough so that inertial, gravitational and frictional effects can be ignored. 

The effect of gravitation in equations (2.21) and (2.22) is that the 
light ions are displaced to regions of the weakest gravity. Parallel elec- 
tric fields are also important in keeping the ions in layers at higher alti- 
tudes. The gradients of density with respect to a can be removed from equa- 
tions (2. 21) -(2. 22) to give the parallel electric field. 
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For an example take the case where n(0 + ) >> n(H + ) and n(0*) >> n(He*) 
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Equation (2.24) shows the parallel electric field decreases the effect of 
gravity on ion distribution because if | is substituted into (2.16) the 
coefficient of the gravity term is less than one. 

Now, since eE \ | is known the ion density distributions can be calculated 
from the equations below. 
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For an example take the case where j is a minor species in a plasma dominated 

by a singly charged ion (H + in 0 + plasma). Also take 3!T./3s ■ 3 T /3s = 0 

a 

then 
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This equation shows that n(H ) can decrease or increase with altitude depend- 


ing on the sign of the term in the brackets. In the topside ionosphere T./T, 

t e 

will be between 0.25 and 1.00; so >?(H + ) will increase with altitude. 
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Equations (2.25) and (2.26) can be simplified with the use of the mean ionic 
mass. 
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Also, the principal ions of the topside ionosphere are singly charged; so the 
expressions for the electron and ion density are 
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(2.32) 


If equations (2.31) and (2.32) are integrated along the magnetic field 
tube from a to a the coupled equations for the electron and ion density 
distributions become 
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Equations (2.33) and 2.34) must be solved numerically to obtain electron and 
ion density distributions for arbitrary electron and ion temperature varia- 
tions. An example of these calculations is shown in Figures 2. 4-2. 7. These 
profiles were obtained from (2.34) by numerically integrating upward from 500 
km using different boundary conditions at 500 km. Equation (2.7) was used 
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as a boundary condition for n(H + ) in one case. In all cases n(0 + ) was as- 
turned to be 1 x 10 cm" at 500 km. Also, different temperature models were 
used for ion and electron temperature profiles. They are shown along with 
n(H + ) and n(0 + ) on each figure. The boundary condition for H + in Figures 
2.4 and 2.5 were applied at 200 km using equation (2.7). This was due to the 
large H* number densities which resulted from using equation (2.7) at 500 km. 
The H + number densities were too large to satisfy the assumptions used in 
equations (2.31) and (2.32). The 0* boundary condition was still applied 
at 500 km. Below 500 km an experimental profile for the 0* number density 
was used in all the cases. 

Figures 2. 4-2. 7 show that H + dominates the ion distribution at high al- 
titudes. The level where the composition transition between H + and 0 + takes 
place is at slightly higher altitudes for atmospheres having larger thermo- 
pause temperatures. 

In Figure 2.6 H + and 0 + number density profiles are shown for the cases 
where H + at 500 km is 1 x 10 3 cm’ 3 and that found from equation (2.7). The 
0 + profile decreases more rapidly in the first case. The calculations for 
the case of n(H + ) ■ 1 x 10 3 cm" 3 at 500 km are not shown in Figures 2.5 and 
2.8 because the profiles are almost exactly the same and the temperature pro- 
files are all about the same from 500 to 1200 km. 

In Figure 2.8 the dependence of ion composition on electron temperature 
is shown. The composition transition level moves to higher altitudes as the 
heat flux and T a + T. increase. This is due mainly to a slower decrease in 

6 t* 

0* densities as T * T. increases. 
e % 

2.8 Low Speed Plasm Flew 

Besides diffusive equilibrium the atmosphere of the low and mid latitudes 
can exist in a state of low speed plasma flow. Indeed, low speed plasma flow 
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Figure 2.4 Profiles of n(0*) and n(H + ) for a 500 K model thermosphere. The ion 
temperature profile is also shown. 
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is important because the coupling between the magnetosphere and the P 2 region 
results in almost continual ion transport. 

The problem of ion flow in the topside ionosphere should not he confused 
with the ambipolar diffusion of 0 + in the P 2 region* In the topside iono- 
sphere several ion species are involved in transport. Also, the scattering 
cross section for ambipolar diffusion is very small compared to the Coulomb 
scattering cross section of thermal protons (dth oxygen ions [Hanson and 
Ortenburgev, 1961], 

The low speed flow momentum follows from (2.19) 


3tt. , tyj 2. 
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where the inertial term is dropped since this case applies to low speed flow. 
(2.30) is valid as long as u « (fcT./m.) i.e. subsonic flow. 

Due to the nonlinear coupling between the ion equations of motion, 
analytic solutions for the ion densities have not been obtained. But in the 
regions below 1000 km H* and He 4 are usually minor constituents. Thus, the 
electron pressure gradient term (in equation (2. 35)) becomes linear (in «(0*).) 
Using the assumptions: 1) a minor ion species, j, diffuses through a sta- 
tionary background ion gas of density n. with n . « n 2) steady state is 

t' J t- 

used in the equation of motion for electrons and major ions; 3) the tempera- 
tures of the two ion species are the same, one can derive analytic v?xpressions 
for the ion density and ion flux [flanks and Kooharte, 1973], 

The equations of motion for electrons and ions are 
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where u . * 0 and v . . is defined such that njn.v . , = njn ,v . . . 
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Eliminating eEi i from (2.38) with (2.36) one finds that 
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where = (W^O/v^. 

1 3n 

In (2.39) — -b- 2. is unknown, but in most real cases the distribution of 

ft o8 

e + 

of electron is assumed to be equal to the 0 density distribution. Also if 

it is assumed that the momentum transferred to 0 + from the minor constituents 

does not affect 0 + distributions, i.e. m.v..w. = 0, then (2.36) + (2.37) gives 
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which gives 
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Substituting this into (2.39) and using the definitions 
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Equation (2.44) in conjunction with the j‘th species continuity equation 
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can be used to find a second-order differential equation in n . which is 
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(2.46) 


where p. is the ion production rate, t. is the loss rate, *(e), y(s), and 
Z(s) are functions derived from quantities in equations (2.45) and (2.44). 

To solve (2.46) numerical methods are required but analytic solutions 
can be found when the production, loss, thermal gradients and gravity varia- 
tions are ignored. The two linear independent solutions are given as 
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where a represents the lower reference point. 


The first solution (2.47) is the solution for the condition «. ■ 0, non- 

diffusing minor constituents. This solution gives a minor ion density which 

increases with altitude if H. < H T./T . At high altitudes the condition 

0 p t. e 

n. « n. is violated by this solution. 

J t. 7 

The second solution (2.48) gives a minor constituent which decreases in 
density with altitude. The diffusion velocity associated with this solution 
can be derived from (2.44) as 
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The case of H + flowing through 0 + has a limiting flux of 
- n.u. - 1.1 x 10 4 [n(H)/«(0)] T ** T. 

t J J IT T* 


(2.50) 
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GeiaUr [1967] has found an analytic expression for the limiting flux 
of H* moving through 0 + for thermal equilibrium. Using a different approach 
his results predict H* limiting fluxes which are approximately the same as 
those presented in Table 2.3. 

Typical limiting fluxes for different thermospheric models are shown in 
Table 2.3. Values in Table 2.3 are based on thermal equilibrium between 
electrons and ions and thus the limiting fluxes are underestimated. Numerical 
solutions give somewhat larger fluxes. 

The previous discussion dealt with approximate solutions to the ion 
momentum and ion continuity equations. H + and He* are not minor constituents 
at high altitudes and therefore the solutions for the density distributions 
of H* and He* are not generally analytic. Equations (2.39) and (2.45) must 
be 


solved by numerical methods. 
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3. THERMAL PROCESSES IN THE TOPSIDE IONOSPHERE 
3,1 Intpoduotion 

In this chapter all the various heat transfer processes in the topside 
ionosphere will be discussed. This includes heat transfer through electron- 
neutral elastic and nonelastic collisions, electron-ion Coulomb collisions, 
ion -ion Coulomb collisions, and ion-neutral elastic and nonelastic colli- 
sions. Also discussed is heat transfer through thermal diffusion, convection, 
and thermal conduction. Particular attention will be paid to the role of 
thermal conduction in determining electron and ion temperatures. 

Electrons created in the photoionization of neutral gas particles was 
recognized as a source of electron gas heating in 1946 [Dmkarev, 1946]. 

When a solar photon ionized a neutral gas particle the excess energy, the 
energy difference between the incident photon and final ion is given to the 
photoelectron. In fact, most of the energy goes to the photoelectron because 
the mass of the ion is much greater than the electron mass. This photoelec- 
tron loses its energy as it travels through the atmosphere. At energies near 
15 eV the photoelectron exchanges its energy with neutral gases through non- 
elastic excitation collisions. At energies near 1.5 eV the photoelectron 
transfers its energy to the ambient electron gas. 

The important sources of energy for the atmospheric gases all originate 
from ultraviolet radiation from the sun (see Figure 3.1). The ambient elec- 
tron gas receives most of its energy through Coulomb collisions with the 
hotter photoelectrons. The thermal electrons, on the other hand, are an 
important source of heating of the ion gases. The neutral gas receives its 
energy through superelastic collisions with excited neutrals produced when 
ultraviolet radiation from the sun is absorbed. 
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Figure 3.1 Sources of energy for particles in the F region. 
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The importance of heating of the ambient electron gas by photoelectrons 
was first discovered by Hanson [1963] who found that the electron tempera- 
ture is substantially greater than the neutral temperature in the F region. 
This is due to the heating of ambient electrons by photoelectrons. Experi- 
mental data from a Langmuir probe on a Javelin rocket shows that the electron 
temperature is greater than the neutral temperature (see Figure 3.2) 

[Hanson et at., 1969]. Incoherent scatter radar measurements [Evans, 1962; 
Hagen and Hsu, 1974], other rocket measurements [Spencer et at., 1965, Nagy 
et at., 1963], and satellite measurements [Smith, 1968, MdClx&e et al. , 1973] 
show conclusively that the electrons are heated to higher temperatures than 
neutrals which is due to photoelectron heating. 

The importance of electron heat conduction was also deduced by Hanson 
[1963] who found that heat conduction is responsible in keeping the electron 
temperature constant even at high altitudes. Photoelectrons produced at 
lower altitudes were found to contribute to non-local heating at higher al- 
titudes [Geisler and Bouhill, 1065), which in turn is conducted back down to 
lower altitudes along the lines of magnetic force. 

Photoelectrons which transfer most of their heat to the ambient electron 
gas do not heat the ion gases because of the mass of ions is much greater 
than the mass of electrons. The thermal electrons are sources of energy for 
the ion gases. Below 300 km neutral gases rapidly cool the ion gases and so 
T. « T [Hanson, 1963] . At high altitudes the energy given to ions by elec- 
trons cannot be dissipated by the ions to the neutrals. In fact, the ion 
temperature approaches the electron temperature at high altitudes [Banks, 
1967a] . 

Thermal conduction for the ion gases is important in keeping the ion 
temperature substantially below the electron temperature [Banks, 1967b]. 


2000*K 3000*K 4000*K 

ELECTRON TEMPERATURE CK) 

Figure 3.2 An electron temperature profile from a Javelin 
rocket experiment [after Hctneon et al. t 1969], 
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Also important is the effect of heat transfer between individual ion species. 
This heat transfer is weak enough to permit significant temperature differ- 
ences between 0* and H* at altitudes near 400 km [Banks t 1967a]. 

In the following sections the important heat transfer processes of the 
ionosphere will be discussed in more detail. 

3.2 Energy Balance Equations 

The energy balance equation for a plasma composed of electrons, several 
ion species, and neutral gases can be derived from Boltzmann's equation as 

n DE /Dt - 3p /3t ♦ V*£ • M + nl* i (3.1) 

where n p is the gas number density. p p » n p kr p is the gas partial pressure, 
T p the gas temperature, u p the rth species transport velocity, K p the rth 
species heat flux vector, F p the external forces. M p is the kinetic energy 
moment of the collision term 


W r " 7 J *°r 

(3.2) 

E r ’ * jVr 2 

(3.3) 

DEjDt - 3 Ejit * U *VE 
¥ ¥ ¥ ¥ 

(3.4) 


The term (3/ /H) is the change in the velocity distribution function 

per unit time for a collision, M p is the mass of the rth particle, and k is 

Boltzmann's constant. The term M p represents the formation of energy per 

unit volume minus the energy loss rates per unit volume due to interactions 

with other eases. M includes both transfer energy due to unequal mean flow 

r 

velocities and unequal species temperatures. This is expressed in the gen- 
eral form M p • l M , where 
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% ■ [avyov + VV*[K ‘ 

+ K“r + m P “p)$p m K )g (i,/a) ] (3 * s) 

for elastic collisions, with the nonentun transfer collision frequency, 
y - |u r - u p | , o ■ 2k(3’ r /m r + T^/m^) , and f(v/ a) and g (y/o) are functions 
determined by the velocity dependence of the interaction cross section p is 
an index indicating a neutral or charged particle. 

The effect of an external magnetic field is present in the flow vector 
which becomes increasingly anisotropic as the ratio of the collision 
frequency to gyrofrequency becomes small. Then most of the heat flows par- 
allel to the magnetic field. Another effect of the external magnetic field 
is the confinement of charged particle flow along the magnetic field lines. 

The term M can be reduced to specific forms depending on the type of 
interactions used. For example, if one considers the hard sphere interaction 
then [Banks and Kockarts, 1973] 
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erf is the error function, and o is the collision cross section, 
interaction, one has [Banks and Kookarta , 1973] 
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where 
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u is the reduced ness given by 
P - mmj (m r ♦ m g ) 

t Z g are the atonic charge of the rth and ath particle, and In A is the fac- 
tor which depends on the minimum scattering angle (it will be discussed later). 
For inverse velocity dependent cross section 

fiv/a) - g{vl o) - 1 (3.12) 

In this case the energy transfer term is 
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Using this result equation (3.1) can be reduced to 
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Equation (3.14) can be further reduced if one uses the fact that the 
heat flow vector is essentially along the magnetic field line. A cylindrical 
coordinate system whose axis is parallel to t and where changes in the inten- 
sity of 5 are taken into account gives 
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where a is a coordinate parallel to t, A is the area of magnetic field tube 

_ i t 

at point a, A * D , and K is the rth species thermal conductivity. 
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3.3 Electron Energy 

The source of energy for electrons is the energy from photoelectrons. 

Since photoelectrons are created from photoionization of the neutral atmos- 
phere, one must know the ultraviolet flux spectrum, the composition of the 
neutral atmosphere, and relevant absorption and ionization cross sections in 
order to determine the photoelectron energy spectrum. Photoelectrons are 
transported and can deposit their energy elsewhere in the atmosphere. Thus, 
the photoelectron flux is needed in order to solve the problem of photoelec- 
tron heating of the ambient electron gas. The solution of the photoelectron 
transport is a difficult problem and will be discussed later. 

Usually when dealing with electrons one divides the energy distribution 
into two components. The first component deals with the ambient electron 
energy distribution. This component is assumed to have a Maxwellian distri- 
bution with a characteristic temperature T g , This component also covers the 
low energy portion of the energy distribution (0.01 - 0.3 eV). The second 
component of the energy distribution is the so-called "tail" of the electron 
energy distribution. This tail contains photoelectrons with energies greater 
than 1.5 eV. 

The ambient electrons gain energy from electron-electron Coulomb collisions 
with a rate given by [Schrnk and Hays, 1971] 

g| - - 2v7 w n g e 4 /^) 5 * In Cj (3.16) 

•j| ■ - 7.7 x 10' 6 n e /J* In C j (3.17) 
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2 4 
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(3.19) 


with In y * 0.5771 > Euler's constant, h ■ Planck's constant and • plasma 
frequency. 

The photoelectron energy loss per unit length is 
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Equation (3.20) can be used to calculate local photoelectron heating rates. 

The electron energy balance equation must include photoelectron heating 
and cooling through elastic and inelastic collisions with neutral particles. 
Inelastic collisions are especially important because cooling through elastic 
collisions is generally slow. This will be discussed in detail later. The 
electron energy balance equation is 

*p, 
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where Q is heat production due to photoelectrons and L is heat loss due to 

p 8 

inelastic and elastic collisions with neutral particles, and Coulomb colli- 

«+. + 2 

sions with ions. jj m n v (u - u ) is the rate of energy production due 
8 88 88 8 8 

to elastic thermal electron collision with ions and neutrals, i.e. Joule 
dissipation. 

3.4 Production of Photoelectrone and Local Heating 

The production of photoelectrons is a very complex problem due to varia- 
tion of the incident flux of solar photons. There have been a number of 
measurements of the solar EUV flux by rocket and sate llite -borne instruments 
with varying degrees of resolution and accuracy [Hall et al. , 1963, 1965, 
1969; Hinteregger, 1965; Hall and Hintereggev, 1970; Hinteregger, 1970], 
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When considering the ionizing effect of solar ultraviolet photons of 
flux T m at the outer boundary of the earth's atmosphere, one has 

T J e " T (3.22) 

where is the ^-th neutral particle number density, is the production 
rate of ionization, of. is the ionization cross section and t is the optical 

T* 

depth given by 



(3.23) 


where c r is the absorption cross section, and the integration is over the 
optical ray path. 

Equation (3.22) gives the ionization rate for a particular neutral spe- 
cies. To find the total ionization rate one must sum over all the neutral 
particles. Thus 
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The concept of the "heating efficiency" [Hanson and Cohen, 1968] can be 
be used to calculate the rate of heat input given locally to the thermal elec 
tron gas. Heating efficiency Is defined as the energy input to the thermal 
electron gas per electron-ion pair produced by photoionization. The photo- 
electron-thermal electron heating rate is given by 

P (s) « c(a) q » e(s) l n.(<A I ) e T ('.24) 

V 

where e(s) is the heating efficiency. 

As an example of equation (3.24) take the case of an isothermal hydro- 
static neutral atmosphere of temperature T , and total density given by 


Em . * n(a) exp(-(*-tf)/P) 
* *7 


( 3 . 25 ) 
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where H * kT^fmg is the scale height, n(a) is density at altitude a, and a 

is the altitude. If o is constant then the rate P fa) is 

a a 

P t (a) * e(a) n(a)c^ I m 

exp(-(a-a)/ff - n{a)Ho a exp (-(a-a)/P)) (3.26) 

This shows that the rate of heat input decreases exponentially with altitude 
if e( 2 ) is constant. The heating efficiency, e(z), plays an important role 
in the solution of the electron temperature. The altitude dependence of e( 2 ) 
should show the effectiveness of neutral particles in taking energy from the 
photoelectrons by inelastic collisions. At altitudes above 250-300 km at 
moderate to high magnetic latitudes the heating efficiency loses its physical 
significance because photoelectrons are transported and deposit their energy 
non-locally. Above the altitude (escape altitude) where transport is assumed 
to be non-negligible, the non-local contribution is added to the heating rate 
[Geisler and Bowhill, 1965; Nagy et at., 1969; Swartz and Nisbet, 1972]. 

The theory of non-local heating of the topside ionosphere by photoelec- 
trons as developed by Geisler and Bowhill [1965] is based on a number of 
assumptions. At altitudes above about 300 km collisions between photoelec- 
trons and neutral constituents are ignored due to low neutral densities. The 
photoelectron pitch angle a is assumed to remain constant as the photoelectron 
is slowed to thermal energy and the photoelectron flux from the conjugate 
ionosphere is ignored. The rate of non-local heating is given by the product 
of the photoelectron flux and the absolute value of dE/ds, where dE/ds is the 
rate of energy loss per unit lengths along the magnetic field [see Eq. 3.20 
In C l = 15] 

^ - Kn/E t (3.27) 

-12 2 2 

where K = 1.95 x 10 eV cm , E is the energy of the photoelectron, n is 


che ambient electron density. 


This method of calculating the non-local photoelectron heating rate 
serves as a basis for further calculations by Mantas [1973]. In his method 
the electron gas heating rate is also calculated from the product of dE/ds 
and the photoelectron flux as 




(dE/ds) * r (E,s)dE 


(3.28) 


where ♦ T (F,a) is the photoelectron flux. But the quantity (dE/ds) depends on 
the energy of the photoelectrons in a more complex way than that given in 
(3.27) since collisions of photoelectrons with neutrals are not ignored. Also 
♦ T (E t s) depends on the photoelectron energy. 
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where v, . is the single collision frequency, v is the photoelectron velo- 

KJ $ 

city, and J7 . is the mean energy lost per collision of the j type. 
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(3.30) 


where n ^ is the density of the fc-th particle, g is the relative velocity of 
of the colliding particles, g ■ |pj - v,|, and o^(c7) is the velocity depend- 
ent cross section for the event j. 

The photoelectron flux energy spectrum has been calculated by Mantas 
[1973]. In his calculations pitch angle distributions and the effect of the 
photoelectron flux from the conjugate ionosphere are taken into account. 

The lower limit of the heat input integral (3.28) is taken to be leV since 
electron temperatures at altitudes below 300 km have energies less than this 
value. In Figure 3.3 profiles of the electron heating rate are shown for 
solar zenith angle x * 0° and 90°. 
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3.5 Other Energy Sources for Bleotrone 

Although photoelectron heating is the major source of heating for the 
ambient thermal electrons, one cannot ignore other possible sources of heat 
for electrons. Electrons can gain heat through interactions with electric 
fields or Joule heating. Electrons can also be heated by hot neutrals. 

Charge particles drift through the neutral gas when there is an electric 
field and a magnetic field applied to a region in the atmosphere. This 
drifting causes collisions between neutrals and charged particles which act 
to limit the final transport velocity of the charged particles. Now since 
the final transport velocity is lowered due to the neutral and charged parti- 
cle collisions, motion is converted into random thermal energy. This energy 
is important in the auroral atmosphere [ Walker , 1971] . The rate of energy 
transfer for electrons is small compared with that for ions. The electrons 
are not heated significantly by the electric fields because the mass factor 
is much greater for ions. This heating is only important in the auroral 
ionosphere because the field-aligned currents are a result of auroral electron 
precipitation. The heating can be ignored at lower latitudes where electron 
precipitation is not significant. 

Joule heating can also be caused by collisions due to differing transport 
velocities of electrons and ions. Thus, electric currents flow and the elec- 
tron heating rate follows from the M equation (3.5) and the Coulomb interac- 
tion equation (3.10). The Joule heating rate is given by 

J e “ m e * e (3.31) 

where g(v/ o) = 3^/4 (a/y) ^ erf(y/a) - 3/2(a/v)^ exp(-v^/o^) (3.32) 

y ■ lu - w. I , = 2 kTjm^ 

1 c t 1 e e 


(3.33) 
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and v Is mean electron-ion collision frequency. 

Electron heating by hot neutrals is possible because of the presence of 
significant numbers of hot neutrals. Energetic oxygen atoms with character- 
istic temperatures as high as 5400 K have been found to have a significant 
number density above 1000 km under quiet solar activity condition and above 
1600 km under active conditions [Rohrbnugh and Nisbet , 1973]. Thermal elec- 
trons can be heated by these hot neutrals through elastic collisions. Cal- 
culations shown below indicate that this source of heat for electrons is much 
smaller than heating by photoelectrons. 

If, n(0) “ 10 4 cm" 3 (3.34) 

<E> Q “ 0.7 eV “ 5400 K (3.35) 

where «(0) is the concentration of the energetic oxygen atoms and <F> Q is the 
average energy of these 0 atoms. In equation (3.35), a Maxwellian velocity 
distribution is assumed to exist, since the mean free path is very long and 
these energetic oxygen atoms have a collision frequency which allows a random 
distribution around the average squared velocity. The energy transfer rate 
for elastic collisions between electrons and 0 atoms is given by [Banks, 
1966a] . 

L = -3.74 X 10" 18 n e «(0) TgVg- T) (3.36) 

A X 

With n =1x10 cm , T = 3000 K, and T » 5400 K, one has 

& & 

L = 4.9 x 10” 3 ev cm" 3 s~*. (3.37) 

eo 

Now, for the photoelectron heating rate at 1000 km one has [Mantas, 1973] 

X « 0‘ Q - 7 x 10" 1 ev cm" 3 s" 1 
X = 90“ Q * 1.5 ev cm" 3 s~* 


(3.38) 
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X is solar zenith angle, Q is photoelectron heating rate. Comparison of L 
with Q indicate that heating through energy transfer between electrons and 
energetic oxygen atoms can be neglected during the day. 

Another source of electron heating is energy gained through de-excitation 
of the vibrational levels of nitrogen which have characteristic temperatures 
of about 3100 K [Walker, 1968]. This nitrogen vibrational energy is pro- 
duced by the collisions with photoelectrons due to the large cross sections 
for vibrational excitation of molecular nitrogen by electrons. This source 
of electron heating is principally effective in the E region (100-140 km) . 

It will be neglected in F-region electron heating. 

3.6 Energy losses of Electrons for Elastic Collisions 

Elastic collisional energy transfer between mixed gases having separate 
Maxwellian velocity distributions is one of the physical processes by which 
electrons lose energy. A deeper physical insight into this process can be 
achieved by considering a simple model of energy transfer. 

First, the average energy loss per collision of a single particle of 
mass and kinetic energy traveling through a gas composed of particles 
of mass n?2 and average energy F 2 is [Huxley and Crompton, 1962] 

Aej * + m 2 ) ] (cj - e 2 ) (3.39) 

The rate at which the single particle loses energy per unit time is given 
by the single particle collision frequency 

v 12 " n 2 g q D (3 ‘ 40) 

where n 2 is the ambient gas number density, g is the relative velocity given 
previously, q ^ is the velocity dependent momentum transfer cross section, 

* 2tt f cr(j7,0)(l - cos©) sinQ dQ 


(3.41j 
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and o(£,0) is the differential scattering cross section. 

The average rate of energy transfer of a single particle follows from 
(3.39) and (3.40) as 


W -*1*2 . 

M " (mj + mp (e l * e 2 )v 12 


(3.42) 


Now, the total average energy transfer exchange rate follows from com- 
bining the average rate of energy transfer of a single particle of Maxwellian 
gas mixed with the original gas as 


dU 

1 . 

W 


J"l m 2 


(m 1 + m 2 ) 


j njlc, - t 2 )v 12 


(3.43) 


°l ' 5"lVl 2 ' Vl 


(3.44) 


where is now the average energy of particles with mass and v^ 2 is the 
average collision frequency which takes into account the many different rela- 
tive velocities between the various gas particles. 

The result (3.43) shows that the energy exchange rate can be decomposed 
into three factors: mass ratio, a difference in average particle energy, and 
an energy transfer collision frequency. The first two quantities are inde- 
pendent of the mode of interaction between the two gas species. The colli- 
sion frequency contains the factors which depend on interparticle forces. 

The exact result for the total average energy exchange rate is given by 
[Desloge, 1962] 


dU x 

dt 


- r 2 ) 


Orij ♦ m 2 ) 2 (2 1 rk) 3/2 (m 2 r 1 + m 1 2* 2 ) 5/2 


r 9 5 q D (g) exp (-Kg 2 )dg 

* ft 


-4irn 1 n 2 ’ 


(3.45) 
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where 

u i ■ j TVi 2 f \ ** V 1 

(3.46) 


K - (2 kT l /m 1 + 2 k!T 2 /m 2 ) _1 

(3.47) 


q D (g) ■ 2s| 0(^,9) (1 - cos0) sin0 dQ 

(3.48) 


U - gas total kinetic energy 

« - particle number density 

m - particle mass 

T - Maxwellian temperature 

k - Boltzmann's constant 

g - relative velocity between particles 

q p(g) - velocity dependent momentum transfer cross section 

v - particle velocity in laboratory system 

d?v - velocity space volume element 

@ - center of mass scattering angle 

ofa.Q)- differential scattering cross section 

f - velocity distribution function 


If one lets 


and 

where 


Q d ■ X s | g S q D (g) exp(-Kg 2 )dg 
7 12 " J n 2 


h f 7&\ 



(3.49) 

(3.50) 

(3.51) 


or if /j and f 2 are Maxwellian velocity distributions 
g - (Sfe/ir)* 4 (fj/mj ♦ T 2 /m 2 ) H 


(3.52) 
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— 3 — S 

and *1 " T ^ E 2 " 7 ^*2 

then 

dU. m.m 9 _ 

W " " 3n l ' .2 * (r i " r 2 )v 12 (3,S3) 

(ntj ♦ m 2J 

There is an arbitrary numerical factor of 4/3 which has been noted by Nioolet 
[1953] in an analysis of electron collision frequencies based upon an analysis 
of collision integrals and diffusion coefficients derived by the velocity dis- 
tribution method. Since m & « m 2 , for electrons becomes 

% 9 ( w’ ) T 1,5 *0°° ex pc- m « t,2 / 2w, # ) <to # < 3 - 54 ) 

6 '0 

To apply this equation to electron-neutral collisions one must have 
analytical expressions for the momentum transfer cross sections. But unfor- 
tunately these results don't exist and one must use laboratory results. A 
description of the current methods used to measure elastic electron-neutral 
momentum transfer cross sections can be found in McDaniel [1964], and Hasted 
[1964]. Empirical expressions are found for q ^ for the atmospheric gases and 
are presented in the following sections. 

3.6.1 Molecular nitrogen . "Pack and Phelps [1961] have measured the 
drift velocity of electrons under the influence of a constant electric field. 
Their data for the electron momentum transfer cross section for electron 
energies between 0.02eV and O.leV can be represented by the equation 

q D - 1.88 x 10“ 15 fik cm 2 (3.55) 

where E is the electron energy measured in eV. For energies above O.leV 
there is a correction which has been found by Englehardt et al. [1964] who 
calculated q ^ from measurements of electron mobilities and diffusion coeffi- 
cients. The data gives 
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q p - (1.83 - . 73 ^)^ x 10* 15 cm 2 (3.56) 

Substituting this equation into the Q D equation gives 

Q d W 2 ) - (2.82 - 9.44 x 10- 3 T*)T* x 10“ 17 cm 2 (3.57) 

which is valid over the range 100 K <, T <_ 4500 K. 

3.6.2 Molecular oxygen. Hie Boltzmann equation can be used to evalu- 
ate the drift velocity, the diffusion coefficient, and the electron mobility 
in molecular oxygen if the cross sections are known. Measurements of elec- 
tron mobility and diffusion coefficients can be compared with predicted 
values. Therefore, the momentum transfer cross sections can be adjusted so 
the predicted and measured values of electron mobility and diffusion coeffi- 
cients agree within 20 percent. Hake and Phelps [1967] have done this analy 
sis and their results are shown in Figure 3.4. As a first order approxima- 
tion to this data one has 

q D - (2.2 + 5.1^) x 10‘ 16 cm 2 (3.58) 

which applies over the range 0.02 <_ E <. 1.0 eV. 

Substituting q^ into the Q D equation gives 

$p(0 2 ) - (2.2 + 7.86 x 10" 2 r **) x 10" 16 cm 2 (3.59) 

which applies over the range 150 K <_ T <_ 5000 K. 

3.6.3 Atomic oxygen. The experimental results for atomic oxygen give 
only total scattering cross sections due to the chemical activity of oxygen 
in any closed container. The momentum transfer scattering cross section can 
be calculated theoretically from quantum theory partial wave phase shifts as 

q n » — <r ? (4 ♦ 1) sin n 

0 V 


(3.60) 
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where ■ {2m 0 V 0 /h) is wave number , h is Planck's constant , % is the angu- 
lar momentum quantum number, and is the partial wave phase shift due to 
the action of the scattering potential. The phase shifts can be calculated 
from measurements of photodetachment cross sections [Cooper and Martin, 1962] . 
If one accepts that the values for the scattering phase shifts are good for 
for determining the momentum transfer cross section for atomic oxygen then 
[Banks, 1966a] 

q D « (3.4 ± 1.0) x 10 -16 cm 2 (3.61) 

or Q d (0) » (3.4 ± 1.0) x 10’ 16 cm 2 (3.62) 

independent of the electron temperature for T < 4000 K. 

3.6.4 Atomic hydrogen. Like atomic oxygen, there are no experimental 
measurements for the electron-hydrogen momentum transfer cross section. 
Theoretical calculations of scattering phase shifts for electrons in hydrogen 
[Smith et al., 1962] can be used to calculate both total cross sections and 
momentum transfer cross sections via equation (3.60). These predicted total 
cross sections agree with measurements of total cross sections oy Braohman 

et at, [1958] and so errors in the calculation of the momentum transfer 
cross sections should not be too large. Using equation (3.60) and the results 
of Smith et al. [1962] and fitting a suitable analytic expression for the 
energy dependence gives 

q D « (54.7 - 28.7 E) x 10~ 16 cm 2 (3.63) 

and Q d (H) ■ (54.7 - 7.42 x lO' 3 ^) x 10‘ 16 cm 2 (3.64) 

over the temperature range 150 K T g <_ 5000 K. 

3.6.5 Helium . The experimental data for electron-helium momentum 
transfer cross sections give a constant cross section of S.6 x 10’* 6 cm 2 
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[Paok and Phelps, 1961] over the energy range 0.0 to 0.4 eV. Thus, 


$ D (He) - 5.6 x 10" 16 cm 2 


(3.65) 


which is independent of the electron temperature over the range 0 <_ <_ 

6000 K. The error is within 10% because of the good agreement between ex- 
perimental data. 

The values of for the different atmospheric gases are shown in Figure 
3.5 which are calculated from the analytic formulas given previously. The 
errors involved in using these expressions can be estimated by considering 
how close they agree with experimental data and by the accuracy of the ex- 
perimental data. For molecular nitrogen and molecular oxygen the uncertainty 
in Q d is approximately 20 percent or less [Englehardt et al., 1964; Hake and 
Phelpe , 1967]. He) is accurate to within 10 percent, but for atomic hy- 
drogen and atomic oxygen 3^ should be less than 25 percent and 30 percent, 
respectively based on quantum calculations. 

3.6.6 Changed particles. The momentum transfer cross section for 
charged particles can be calculated from equation (3.48) and the Rutherford 
differential scattering cross section 

2 

cosec 4 (6/2) (3.66) 

where p is the two particle reduced mass, 2 j ^ are the respective atomic 
charges, e is the electron charge and 9 is the center of mass scattering 
angle. Integrating o(e,^) over e, one gets 


®( 9 . 0 ) 


* 1 * 2 ® 

J 


■ 


8ir 


w 

2 ug J 


In 


COS 0 _ 


( 3 . 67 ) 


m 


where 9 is minimum scattering angle determined from the parameters of the 


m 


charged particles. The relation between the impact parameter, b, the 
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scattering angle, 6, and the relative velocity, g, is given by 
(1 - cos 6) ? ■ 

1 ♦ (&/2> 0 r 


(3.68) 


2 2 

where = 2^0 /\ig , other quantities defined previously. The minimum 
scattering angle 6^ can be calculated from this equation when b is a maximum, 
by Calculation of the maximum impact parameter, by follows from consider- 
ing the shielding of charge due to electrostatic polarisation effects [Cohen 
et al ai 1950], The Debye length, D , which represents the maximum distance 
over which microscopic density fluctuations are correlated by random electric 
fields, can also be used as the effective range of the Coulomb interaction 
and so can be used in equation (3.68) as the maximum impact parameter. Since 


<D 


16ir 


2 ' 

2 

r 2 _ 1 

w 

In 

W D e 

1 

cm 

CM 

1 

2 

LV2* j 


(3.69) 


which is the momentum transfer cross section between two arbitrary charged 
particles. 

The argument of the logarithm in equation (3.69) can be expressed in 
terms of the center of mass energy of the colliding particles, E, as 


ug\ 2 ED e 

1 = 

z^ Z 2 e z \ z 2 e 


(3.70) 


Calculations of In A for various particle energies and Debye lengths are 
shown in Table 3.1 where 2 j = 2 2 = 1. For the earth's atmosphere particles 
energies and Debye lengths are such that the approximation In A = 15 can be 
used with an uncertainty of less than 10%. 

Using (3.70), (3.69) and (3.49) for two gases with different Maxwellian 


Values of 
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temperatures one gets 


Q D mj * 


In A 


(3.71) 


(jymj ♦ T 2 /m 2 r 

which applied to charged particles of arbitrary mass and charge. 

3.6.7 Electron collision frequencies and energy transfer rates . The 
momentum transfer cross section for electron-neutral elastic collisions 
derived in the previous sections can be used to calculate electron energy 
transfer collision frequencies and energy transfer rates via equations (3.50) 
and (3.53). The final results are given in Tables 3.2 and 3.3. 

For charged particle the electron energy transfer collision frequency 
and energy transfer rate has been calculated by Banks [1966a]. With 
In A * 15, the collision frequency is 


v = (54 ± 5) n./JF ' ^ sec~* 
e v % e 


(3.72) 


and the energy transfer rate is 


dU. 


-^r=~ (7.7 ± 0.8) x 10" 6 w^» 


e i 


Ve-h'* „ -3 -1 

jjj ev cm sec 


(3.73) 


A . T 
% e 


where A. is the ioA atomic mass in amu. 


Note in equation (3.73) that for a fixed I T. there occurs a maximum 


energy transfer rate at T =3 T.. If the temperature dependence of the term 

Q %r 


In A where included in (3.73) one has instead T g = 3.5 T.. Equation (3.73) 


is important since it expresses energy transfer between charged particles in 
a single equation. 

3.7 Energy Losses of Electrons Through Inelastic Collisions 

Besides elastic collisions, inelastic collisions are also important for 
electron cooling. Electron*: can excite neutral particles and thereby trans- 
fer some of their energy to the neutral particles. Rotational, vibration, 
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TABLE 3.2 


Electron energy transfer collision frequencies for elastic 
collisions with atmospheric neutral particles . 


Neutral 

species 

Collision frequency 

N 2 

(2.82 - 9.44 x 10 -3 tSt % x 10 

& 6 

°2 

(2.2 ♦ 7.86 x 10" 2 T S x 10 -16 

0 

(3.4 ± 1.0) x 10” 16 

H 

(54.7 - 7.42 x 10' 3 T ) x 10 -16 

H 

5.6 x 10' 16 


-17 
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Neutral 

Species 

N„ 


H 


H 


TABLE 3.3 

Electron energy transfer rates for elastic collisions 
with atmospheric neutral particles . 

Energy transfer ^ e . ,, -1 -3* 

rat. - ~3t CeVs ] 

nCN,)n (T - T) 1.77 x 10' 19 (1 - 3.35 x 10" S T 

fc v e W Q 6 

n(0~)n (T - T ) 1.21 x 10 _18 (1 + 3.57 x 10“ 2 tSt* 

C & 71 & & 

n (°) n e^ e - T r? 3 * 74 x 10 ’ 18 T e H 

n(H)n e (T e - TJ 9.63 x 10' 16 (i - 1.35 x 10 -4 T 

”<V n « (J '* ' T n> 2 ‘ 46 * 10 ‘ 17 *1* 
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electronic, and fine-structure excitations by thermal electrons are possible 
for the atmospheric neutral particles N., 0-, 0, H, and He. 

Electronic excitation of N 2 is not a significant energy loss process 
for electrons in the ionosphere due to the relatively large excitation poten- 
tial of the first excited state (6.2 eV) . The electronic excitation poten- 
tials of H and He are larger than for N 2 so energy loss through electronic 

excitations of N 2 , H, and He can be neglected in the ionosphere. For 0 the 

1 3 

first excited state, 0( D) lies 1.97 eV above the ground state 0( P) . Thus, 

for electron temperatures greater than about 23000 K electron impact excita- 
tion of the state can be important. Since electron temperatures in the 
F region are much less than 23000 K, electronic excitation of 0 is also not 
important. 

The excitation of the second excited state of atomic oxygen 0(*S) is 
not a significant energy loss process for electrons since the excitation 
potential of 4.2 eV is too large for ionospheric electrons. 

The electronic excitation potential for 0 2 (a*A ) lies only 0.96 eV above 
the ground state. But due to very small excitation cross section leading to 
O^a^Ag) [Schulz and Dowell, 1962] this is not a significant sink of electron 
thermal energy. 

Electron energy loss rates from electron-impact-induced transitions of 
the fine-structure levels of atomic oxygen are important because the ^P^ and 
levels are 0.028 and 0.020 eV above the ground state 3 P 2 level and must 
be included in the electron energy balance equation. This cooling is brought 
about by the collision process 

e + 0( 3 Pj) ■+ e + 0( 3 P Jt ) (3.74) 

The distribution of levels J follows from assuming that the reaction 
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(3.74) proceeds rapidly. The distribution is given by [Dalgamo and Deggee, 
1968] 

n( 2J + l) exp {-EJkr ) 

- - 1 - 2 (3.75) 

JmQ W ♦ 1) «pHy*r„) 

where n is the number density of atomic oxygen , Ej is the energy of the state 
with total angular momentum J and T n is the neutral particle gas temperature. 
The time rate of change of the electron energy E through excitation and de- 
excitation of the fine structure levels of atomic oxygen is given by 

2 

§ - v do n J J'V < 3 - 76 ’ 

where v is the electron velocity, and Q(J, J' , E ) is cross section for col- 
lision (3.74). The cross section Q(J, J' , E ) can be calculated from theoreti- 
cally derived collision strengths. Breig and Lin [1966] were the first to 
derive these collision strengths. More recently Tambe and Henry [1976] and 
LeDoumeuf and Nesbet [1976] have recalculated these collision strengths. 

Cross sections calculated from these more recent data are shown in Table 3.4. 
The mean rate of change of E can be calculated from [Dalgamo and Degges t 
1968] 

I - 2 ' (' W e )' 3/2 (2 1*,/ 12 X V 


exp(-ff/fcT e ) (Ej - E J% ) Q{J , J\ E) 


dE 


(3.77) 


where the electrons were assumed to have a Maxwellian velocity distribution 
with characteristic temperature T . The mean cooling rates are shown in Fig- 
ure 3.6. To a first approximation the cooling rates of Figure 3.6 can be 
represented by 



Total cross sections for de-excitation of the 
spin-multiplet levels of 01 by electron impact 


Transition ( 



Q(J+tT') = 


ENERGY LOSS RATE (eV cm 3 sec’ 1 ) 
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-11 -1 
■ -3.4 x 10 n g n(0) [5*3 exp ( -ErfkTj ♦ exp(-E 3 /fcf M )] 

x T g 1/2 {(.81 + 8.6 x 10 -5 T g ) F 2 [exp(-E 2 /W e ) - exp(-£ 2 /fcr n ) ] 

+ (.32 + 8.7 x 10" 5 T g ) E z [exp(-E 5 /kT e ) • exp(-E 3 /kr n )] 

+ (.32 + 8.9 x 10' 5 T g ) E^expt-E^kTJ - exp^/kT^)] } 

(3.78) 

where ^ « 0.00833 eV, E 2 - 0.01965 eV f ffj « 0.02798 eV, and k = 8.62 x 10" 5 
eV X" 1 . The integrals in equation (3.77) are approximately equal to coefficients 
of the terms [exp(-E/kf ) - exp(-F/kT )] in equation (3.78). The cooling rate is 
proportional to (T g - T ) as T g approaches 7^. 

Election cooling by vibrational excitation of 0 2 has been calculated by 
Lane and Dalgamo [1969]. Their results are shown in Figure 3.7. Although 
the results show that the heat loss by vibrational excitation of 0 2 is not a 
dominant heat loss process, it may contribute up to 20% of the total hea* loss 
in the E region. So, it must be included in any accurate quantitative study 
of the thermal balance in the E region. The electron cooling rate by vibra- 
tional excitation of 0 2 can be represented by 

10~ 16 «(0 2 ) n e {T e - T n ) T e ltl eV cm -3 s _1 (3.79) 

The vibrational excitation of N 2 is a more efficient electron cooling process 
than rotational excitation of N, for electron temperatures above 1500 K 
[Dalgamo and Fenny, 1965], In Table 3.5 cooling rates due to vibrational 
excitation of N 2 from Fees et at. [1967] are presented. 

Electron cooling through rotational excitation of 0 2 has been calculated 
by Dalgamo et at. [1968], They used the expression 






TABLE 3.5 

Cooling rates through vibrational excitation 
of Nj (after Dalgamo at at, 1967). 



200 3 x 10' 19 2500 1 x 10* 11 
300 7 x IQ* 17 3000 5 x 10" 11 
500 7 x 10‘ 1S 3500 1 x 10 -1 ° 
1000 2 x 10‘ 13 4000 3 x 10 -10 
1500 8 x 10" 13 4500 2 x 10“ 9 
2000 2 x 10" 12 5000 7 x 10’ 9 
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eV cm* 3 sec' 1 
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te 



- T nftn - I in* 14 (r# " r,l) 

- -v (0 7 x 10 — 7772“ 


(3.80) 


where n(0 2 ) is the number density of 0g end the other quantities have been 
given previously. 

Rotational excitation of Ng has been studied by Dalgamo and Moffett 
[1963) . The computed cross sections are in satisfactory agreement with the 
analysis of the experimental swarm experiments of Engelhardt et al. [1964]. 
They can be represented by the cooling rate 


-ft - -voy 2 x 10* 14 - T n )/T 0 l/2 ev cm* 3 s’ 1 (3.81) 

with a possible uncertainty of a factor of 2. 

3.8 Thermal Conduction 

Particle collisions determine the thermal conductivity. In ionized 
gases charge separation electric fields have the effect of limiting the flow 
of electrons. Thus, electric fields tend to reduce the effective thermal 
conductivity. 

Owing to the dependence of the thermal conductivity of a gas of charged 
upon the inverse square root of the mass of the particle, the thermal con- 
ductivity of the electron gas is much larger than that for the ion gas 
[ Chapman and Cowling , 1970; Spitzer, 1962], Nevertheless, the ion thermal 
conductivity cannot be ignored comp' itely since at high altitudes the rate 
of energy transfer from electrons to ions is slow. Thus, ion conduction will 
be important in determining ion temperature profiles. 

3.8.1 Electron thermal conduction. The importance of the electron 
thermal conductivity in the electron energy balance equation was studied by 
Geieler and Bowhill [1965]. In their work three situations that probably 
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could be found within the ionosphere were investigated. The first case is 
that of an overheated electron gas at all altitudes. In this case the effect 
of thermal conduction is to conduct downward the excess heat through a strong 
electron temperature gradient to the heat sink in the lower thermosphere. 

The second case is that there is heat sink at high altitudes where heating of 
electrons is slow. In this case even a small temperature gradient will re- 
sult in a heat flux that is sufficient to exceed the capacity of the heat 
sink and so the electron gas is not cooled by the heat sink. The third case 
is t..at the electron gas is only overheated in the region between 200 and 
250 km. This case leads to a temperature maximum around 200 km and the tem- 
perature is controlled by neutral cooling rather than by thermal conduction 
at the maximum. 

The calculation of the electron thermal conductivity is carried out by 
considering the mean free path within a weakly ionized gas. This permits the 
separation of the collision effects of neutral and charged particles colli- 
sions. The effective thermal conductivity can then be found from the indi- 
vidual conductivities of neutrals and electrons. This has been executed by 
Bartke [1966c] who arrived at 


7.7 x 10 5 T g S/2 


P _ 

1 . 3.22 x 10 4 J- l « n S B 
e 


,, -1 -1 „-l 

eV cm sec K 


( 3 . 82 ) 


where K is the effective electron thermal conductivity, the summation is 
e 

over all neural gas species present, Q ^ is the average momentum transfer 
cross section. For low neutral particle number densities this expression 
(3.82) reduces to 


K • 7.7 x 10 S T S/2 eV cm" 1 sec" 1 K" 1 
e e 


( 3 . 83 ) 
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which is the equation derived by Spitae r [1962] for a fully ionized gas. 

The effect of the neutral atmosphere on the electron thermal conductivity 

for several values of T is shown in Figure 3.8. Hie transition between the 

e 

fully ionized case and the partially ionized case of (3.82) occurs when the 
curves depart from the horizontal. 

The effect of the earth's magnetic field on the thermal conductivity is 

important since the thermal conductivity is anisotropic with different values 

parallel and perpendicular to the field lines [Chapman and Cowling , 1970]. 

The value of the thermal conductivity is unaffected when temperature gradients 

are parallel to the magnetic field. The perpendicular component is reduced 

2 2 

by the factor v /(v + Q ) [Chapman and Cowling , 1970] where v is the mean 

B B B B 

collision frequency and ft is the electron cyclotron frequency or gyrofrequency. 

I? 

2 2 -9-1 

In the earth's atmosphere v /(v + n ) « 1 since v * 10 n sec and 

& b b b n 

JJ * 10^ sec. Thus, only electron temperature gradients parallel to the field 

B 

lines will be effective in transporting energy through the electron gas. 

At high altitudes the electron temperature profile can be determined 
from thermal conduction alone. The heating and cooling terms in equation 
(3.14) can be neglected. If steady state is assumed and the effect of exter- 
nal forces is ignored, then the electron temperature can be calculated from 


a * T 0 
| - K -J- = 0 
32 e 3 2 


Integrating one gets 
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where the subscript a refers to some reference altitude. Integrating again 


C ' 


42 - 


gives 
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Figure 3.8 The effect of the neutral atmosphere on the 
electron thermal conductivity [after Banka 
and Koakarte, 1973]. 
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showing that the electron temperature must continually rise with altitude. 


The presence of the heat flux can be inserted into (3.85) giving 


T e * 0 ?/2 + 7 *\\ (8_a) : 


(3.86) 


where $|| is the heat flux at some reference altitude. 

3.8.2 Ion themal conduction. Although ion thermal conductivity is not 
important at lower altitudes it can be effective in keeping the ion tempera- 
ture significantly below the electron temperature at high altitudes [Banks, 
1967b). Thermal energy flows along the lines of magnetic field when there is 
a gradient in the ion temperature. Generally, heat flux for a single ion is 


given by 


77 i 


(3.87) 


but diffusion and thermal diffusion can alter this expression when the ion 
temperature is small and when ion densities are large. The importance of 
diffusion heating or cooling due to ion transport can be calculated by com- 


paring the term 


37 * • 

(f n e k) 


(3.88) 


where u is the ion diffusion velocity, s is the vertical coordinate, and n e 
is the electron concentration. In Table 3.6 are shown some calculations of 
diffusive heating and thermal flux from conduction where the ionization fluxes 


have been taken from Table 2.3. 


TABLE 3.6 

Comparison between diffusive heating and thermal flux from conduction for ions with 
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The calculation of the ion thermal conductivity for a single ion gas 
of atomic mss has been performed by Chapman and Cowling [1970] and is 

- 4.6 x 10 4 (2\ 5 /^) 1/2 eV cm -1 sec’ 1 K" 1 (3.89) 

where the second approximation correction has been used. Below 500 km the 
ion conductivity is damped because of ion-neutral collisions, and is divided 
by [1 ♦ 3.22 x 10 4 (T i 2 /n i ) Z Q D ). 

When there is a mixture of ions in the ionosphere and a single ion tem- 
perature is assumed, then the thermal conductivity must be adjusted to take 
into account the different values of the thermal conductivity for each indi- 
vidual ion. The exact expression is difficult to obtain but a density 
weighted thermal conductivity can be used as [Banks, 1967a] 

K. - 1.2 x 10 4 [n(0 + ) ♦ 2n(H + ) ♦ 4n(He + )]/n (3.90) 

w & 

which is accurate to within several percent in binary mixtures. (n(0 + ) is the 
number density of 0 + ions, n(H + ) is the number density of H + ions, and n(He + ) 
is the number density of He 4 ions). When there is a mixture of ions and sepa- 
rate ion temperatures are calculated, then equation (3.89) is used for thermal 
conduction of each individual ion. The solution for the separate ion tempera- 
tures is discussed later in Chapter 6. The effects of ion-neutral .ollisions 
which is included in the electron thermal conductivity can be neglected abov~ 
300 km. 

3.9 Energy Sources for Ions 

The major source of heating of ions is energy transfer through Coulomb 
collisions with thermal electrons. The energy transfer rate has been calcu- 
lated previously and is given in equation (3.73). 

In addition to heating by the electron gas, ions can be heated through 
dissipation of wave energy, or through electric fields, chemical reactions, 
or inelastic collisions. But experimental observations indicate that the 
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electron gas is the primary source of ion heating altitudes above 300 km. 
Incoherent scatter measurements of Foana [1969] and Began and Bsu [1974] show 
that the ion temperatures are always equal to or less than the electron tem- 
perature. This implies that the electron gas is the major ion heat source. 

Different ion species can have separate temperatures . H* ions are heated 
to higher temperatures than 0* ions due to the larger collision frequency of 
electron with H* ions. Thus, H* ions can transfer energy to the cooler 0* 
ions through Coulomb collisions, and H* ions are then an energy source for 0* 
ions. 

Below 300 km, the effects of neutral winds, electric fields, and elec- 
tric currents can cause ion heating through conversion of bulk transport into 
random thermal motions. This is ion Joule heating. Also, since heat produc- 
tion Mtd energy loss are governed by the same collisional processes, the ion 
temperature can be found from (3.14) when steady-state conditions apply as 

T. - T + m (u . - uj 2 /3k (3.91) 

t n n t n 

Generally, electric fields and neutral winds are not large enough to have a 
significant effect in raising the ion temperature above the neutral tempera- 
ture. 

3.10 Energy Leases for Ions 

There exist for ions three methods of energy loss: ion-neutral elastic 

collisions, ion-ion Coulomb collisions, and ion-neutral resonant charge ex- 
change. The process of ion-neutral resonant charge exchange is fundamentally 
elastic in nature but it cannot be described by the equation of collisional 
energy transfer. This is discussed later in Section 3.10.2. The analysis 
of laboratory and theoretical cross section data are used to obtain the res- 
onance ion-neutral energy transfer rates. In the following subsections the 
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collision frequency and energy transfer equations are discussed for elastic 
ion-neutral and ion-ion collisions. 

3.10.1 Ion-neutral oolliHons. The energy transfer rate for two gases 

with Maxwellian velocity distributions characterized by the temperatures 2^ 

and ?2 11114 particle masses Wj and m 2 is given by equation (3.53). Accurate 

expressions for both energy transfer rates and collision frequencies can 

easily be derived from equation (3.53) once the velocity dependent momentum 

transfer cross section is known. At low temperatures the induced dipole 

2 2 

attraction with a potential of Q • - a e fir is the most important ion- 

j 

neutral interaction (a is the neutral atom polarizability, and r is the ra- 
dial separation) . Using this relation it is found that the average momentum 

- 1/2 

transfer collision cross section is proportional to (2^ ♦ ' and that 

the collision frequency is independent of the temperature. 

The induced dipole force of attraction is countered by a short -range 
quantum mechanical repulsion at temperatures greater than 300 K. The nature 
of this repulsive force is linked directly to the details of the ion and neu- 
tral orbital electron structures. This creates a problem since there are 
almost no data available for either the cross sections or collision frequen- 
cies for ion-neutral collisions at high temperatures. Thus, ion-neutral 
energy transfer rates are somewhat uncertain. The rates are estimated by 
using the low temperature polarization interaction with a suitable truncated 
repulsive potential to account for short range quantum effects. 

dalgamo et at. [1968] has evaluated the ion-neutral interaction po- 
tential by assuming the polarization potentials at distances greater than a 
small atomic radius and an elastic sphere repulsive potential at distances 
less than this atomic radius to simulate quantum repulsive effects. The 
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ion-neutral collision frequency for singly charged ions is given by 

v^ n - 2.6 x l<f 9 » n (« tf / u <n jl/2 8ec " 1 ( 3 * 92 ) 

-24 3 

where is the neutral gas atomic polarizability in units of 10 cm , and 

u. is the ion-neutral reduced mass in ami. The ion-neutral elastic energy 
tn 

transfer rate is calculated from equation (3.53) as 

dU * ■ n. . i 

3F " * 2 * 6 x 10" 2 v * n * V « v « sec 

{A i 4 V 


(3.93) 


where A is the atomic mass of the indicated particle. In Table 3.7 are 
shown the ion energy loss rates to the atmospheric gases where the polariza- 
bilities are from Banks and Koohxrts [1973]. The differences in the rates 
are due mainly to changes in the mass factors rather than changes in the ato- 
nic polarizability. 

3.10.2 Ion-ion collisions . In studies of ionospheric thermodynamics a 
single ion temperature is often used throughout the ionosphere . However, 
Banks [1967a] has shown that a temperature difference exists between H 4 and 
0* ions and between He 4 and 0* ions. This difference is largest for H 4 and 
0* ions and has a maximum of 200 K in the region between 250 and 650 km. 

The important factor limiting the size of the temperature difference is the 
rapid ion-ion energy exchange rate. The Coulomb interaction determines the 
ion -ion energy transfer frequency which follows from equation (3.71). The 
ion- ion energy transfer rate is calculated from equation (3.53) using equa- 
tion (3.71) as 


3.3 x 

m 


10 


-4 


n l n 2 


Pi - T, 


2 5 

m 


eV cm sec 


V 2 [(W 4 (W 1 


(3.94) 


TABLE 3.7 


Ion energy loss rates for elastic collisions with 
atmospheric gases (after Banka, 1967b). 


Neutral 

Species Ion Energy 



loss rates (10"**eV c»" 3 sec**) 

6.5 n(0 + )n(N 2 )(r € - TJ 

5.8 n(0 + )n(0 2 ) (2\ - TJ 

2.8 N(0 + )n(H e )(!T i - T n ) 

3.1 n(N 2 )n(H + )(T £ - 2* n ) 

2.8 n(0 2 )n(H + )(2\ - iy 

5.5 n(H e )n(H + )(^ - TJ 
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where the subscripts apply to the two ion species and A is the ion mss in 
ami. 

The H + ions lose energy to 0* ions and the 0* ions in turn gain energy 
through 0 + - H + Coulomb collisions. This is because H* ions are heated 16 
times faster by electrons as shown in equation (3.73). 

3.10.3 Resonant ion-neutral oollieione. The importance of charge ex- 
change in the ion energy balance equation has been investigated by Banke 
[1966b]. In process of exchanging charge both of the original particles tend 
to retain its original kinetic energy. Thus, the reaction is fundamentally 
elastic in perserving the total kinetic energy even though the identity of 
the ion has changed. Thus, charge exchange provides energetic ions a rapid 
way in which to transfer energy to neutral particles. For many atmospheric 
ion-neutral pairs the probability of charge exchange collisions is much larg- 
er than that for normal elastic collisions. 

The energy loss rate of an ion gas is derived from considering energy 
balance in ion production and loss. For Maxwellian velocity distributions 
with characteristic temperatures , T. and T the ion-neutral energy loss rate 
is [Banke, 1966b] 

% > - 7 »i* V T i - V < 3 - 95 ’ 

where is the ion-neutral charge exchange collision frequency. In Table 
3,8 are given energy loss rates for ions of the topside ionosphere. 

Besides resonant ion-neutral charge exchange collisions between atmos- 
pheric gases of the same species, there exists an accidentally resonant charge 
exchange between H* and 0. This has been discussed previously in Sec ion 2.2 
with regard to hydrogen-ion production and loss. The energy loss rates for 
oxygen and hydrogen ions due to charge exchange are [Banke, 1967a] 


TABLE 3.8 


Resonance charge exchange energy loss rates. 


Species 


0 + , 0 



1 1 1 

n i "n (T. ♦ tJ* 

t ft 

2.1 

13.0 


3.9 
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dU i 

~d? 


. ■ 3.8 x 10" 1S n(0*)n(H)T 1/2 [r(0 + ) - 

0 -H " 


l a ffl to flCl r'Kf ^(H 4 )] eV cm" 3 sac- 1 
9 n(0 )n(H) n 


du € 

St 


H + -0 


3.4 x IQ' 15 n(H 4 )n(0)r 1/2 (H fc ) 


(3.96) 


x [T( 0*) - ? n t H - M° . ) (T $ /T[ H + )) 1/2 ] eV c*“ 3 sec’ 1 
8 n(H )n(0) n 


(3.97) 

The resonant charge exchange energy loss rates are based on experimental 
results reported in [Banka, 1966b). While the maximum error in He 4 - He 
charge exchange cross sections is less IS percent, the maximum errors associated 
with 0 4 -0, 0 4 -H, H 4 -0, and H 4 -H charge exchange cross sections are less than 
11, 10, 10, and 1 percent respectively. 
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4. INCOHERENT SCATTER THEORY 

4 . 1 Introduction 

Information about the temperatures of the ions and the electrons and 
about the mixture of ions can be obtained from observations of the scattering 
of electromagnetic waves from density fluctuations in the medium. Gordon 
[1958] was the first to predict that the width of the frequency spectrum would 
be determined by a Doppler spread caused by the thermal velocities of the 
electrons. Later observation? oy Boalea [1962] and Pineo et at., [1960] 
showed the frequency spectrum to be much narrower than anticipated because the 
Doppler spread is determined by the thermal velocity of the ions. The theory 
of this phenomenon has been developed by a number of workers \Fejer, 1960, 

1961; Dougherty and Farley, 1960; Renau t 1960; and Salpeter, 1960a, b] and 
all arrived at identical conclusions. 

Fejer [i961] has developed a theory of scattering of radio waves bv an 
ionised gas in thermal equilibrium in the presence of a uniform magnetic field. 
This is a generalized treatment in that the results would be used for s'veral 
types of ions present and for unequal values of the individual ion tempera- 
tures and electron temperature. 

Moororoft [1964] has studied the effects of ionic composition on the 
nonmagnetic spectrum. From his work Hagen and Heu [1974] have deduced the 
structure of the protonosphere above Arccibo. 

In this chapter attention will be paid to the effects of unequal ion tem- 
peratures on the spectrum where a magnetic field is not present which is rele- 
vant for all directions of propagation except those nearly perpendicular to 
the magnetic field lines. The ions, 0* and H*, will be treated and considered 
in detail. The ion, He*, is not considered in this treatment because of its 
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of its relative unimportance with respect to H 4 and 0 4 . Experimental obser- 
vation of He 4 indicate that its number density is less than ten percent of the 
H + number density [Rees et al., 1971; Gleeson and Ax ford, 1967]. 

The effects of ionic composition, electron temperature, and Debye length 
on the nonmagnetic spectrum will also be considered. The H 4 number density 
plus the 0 4 number density is assumed to equal the electron number density. 
Electron temperatures are assumed to be greater than or equal to the ion tem- 
peratures. This has been shown to exist by a number of experimental observa- 
tions of the topside ionosphere [Hagen and Hsu, 1974]. Theoretical expres- 
sions for scattering from the thermal density fluctuations in an ionized gas 
have been taken from a paper by Fejer [1961]. Computations in this work were 
programmed on the H-P 9830, the CDC Cyber 1700, and the IBM 360 at the Univer- 
sity of Illinois. 

4.2 Incoherent Scatter from a Thermal Unmagnetized Plasma 

The influence of ions on the scattering characteristics has been shown 
to be important when the exploring wavelength is very much larger than the 
Debye length [Dougherty and Farley, 1960; Salpeter, 1960a; Hagfors, 1961]. 

The Debye length or shielding distance D is 

D = (e n kT In e"^-^ meters (4.3) 

e o e e 

where is the permittivity of free space, k is Boltzmann's constant, T g is 
the electron temperature, n is the electron density, and e is the charge on 
the electron. The scattering can be thought of as arising from density fluc- 
tuations brought about by longitudinal oscillations in the plasma. Ion- 
acoustic waves and electron-induced waves at the plasma frequency and the 
electron gyrofrequency comprise a major part of these oscillations. In the 
absence of collisions the frequency spectrum of the scattered power from 
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randomly distributed electrons is 

2 


where 


and 


S(w) 




1 - il c t 


exp(-u 2 /nj 2 ) 


1 A ^ 

J»1 2 


jl 2 n /‘ *j ", 


*j ■ V m i ,n j )in 

Cj - P^" 2 fcy’ 2 [-1 ♦ 2(u/fl^) exp(-«j 2 /0^ 2 ) 


X exp y 2 dy * i w 1/2 (u>/SK) exp(-<o 2 /fij 2 )] (4.4) 


(4.2) 

(4.3) 


*i • (e o^/V 2)1/2 


k w - 4w/x 


(4.5) 

(4.6) 


where X is the transmitted wavelength 


Subscript 1 refers to electrons, 2 refers to 0*, and 3 refers to H*. The 
Thomson scattering coefficient of an electron is 

" CV* (4.7) 

where w Q is the permeability of free space, y is the angle between the elec- 
tric field vector of the incident wave and direction of scattering, and the 
other quantities have been defined previously. 

Tn equation (4.2) the first term is called the electronic term and 
second term is called the ionic term. For large values of the parameter 
a - 4sP^/X the spectrum function is very broad and most of the scattered 
power is in the electronic part of the spectrum 1961]. This spectrum 
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has a Gaussian shape. For small values of a the spectrum becomes narrow and 
the ionic part of the spectrum contains most of the scattered power. It is 
this feature which allows radar investigation of the ionosphere using appa- 
ratus of lower sensitivity. In the topside ionosphere a is generally quite 
small. At larger values of a the scattered power can be calculated by use of 
approximations. 

The only term in equation (4.2) which is not analytic is the integral 


«i)/n 

l ' 


7 2 

* exp y dy 


(4.8) 


This can be expressed as the power series 

2 k - 1 


0) 


u» r.2tn2\ t. (fc-1)! -0) 

ff exp( “ /n ) fc£i 

(see Appendix B), or as an asymptotic expansion 

o 2 2 

exp(u /n ) 


(4.9) 


2 2 2 2 3 2 4 
n fl Or tr 

x (1 + --a- ♦ 3 -— *■ ♦ 15 -iiy + 105 A? ♦ ...) 

2u 2« 2u 2w 

Both the power series and the asymptotic expansion used in the calculation 

2 2 

of the power spectrum depends on the magnitude of the quantity, m . When 
2 2 

u /£! is small, the power series is very useful since it converges rapidly. 

2 2 

When id /n becomes large and the power series fails to converge, the asymp- 
totic expansion is used with very little error entering into power spectrum 
calculations. 

Buneman [1962] has shown that the total cross section attributable to 
the ionic component when T » T . is given by 

6 u 


a * o /(I ♦ a 2 ) (2 ♦ a 2 ) 
& 


(4.10) 
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which approaches of 2 as a •» 0. In the topside ionosphere, T is often not 
equal to T.. Generally T/T. is in the range 1 to 4. When this is true the 
effect of varying T/T. on the total echo power can be expressed as 

o - o e /(l ♦ o 2 )(l ♦ TjT i ♦ a 2 ) (4.11) 

The exact variation of the total cross section a as a function of T/T , 

0 % 

has been calculated by Mooraroft [1963] and is shown in Figure 4.1 for the 
case a -► 0. 

4.3 Characteristics of Theoretical Spectra 

4.3.1 Spectra for one type of um . In order to better understand the 
more complicated case of two ions, spectra for one type of ion will be re- 
viewed and examined first. Power spectra are shown in Figure 4.2 
for scatter from ionized gases where the ion present is either oxygen or hy- 
drogen for electron-ion temperature ratios T /T. of 1.0, 2.0, and 3.0. Only 

0 S' 

half of the spectrum is shown since the spectrum is symmetric about the 

transmitter frequency. The scale for the power scattered per unit incident 

power density (power spectral density) is normalized by the power spectral 

density at zero Doppler shift. The frequency scale is normalized by the fac- 
1/2 

tor X(m /8kT ) . This indicates that the position of the peak in the spec- 

1/2 

trum depends on T ' , since the peak position doesn't change significantly 

when the frequency scale is normalized by this factor. Figure 4.2 also shows 

that the width of the spectrum is nearly proportional to the square root of 

1/2 

the ion mass. Thus, if the frequency scale is normalized by m. as is 
shown in Figure 4.3, then the peak of the spectrum appears at almost the same 
Doppler frequency shift for both 0* and H + . 

The width of the spectrum also depends on the ion temperature. This 

is shown in Figure 4.4 for 0 + where T/T, is assumed to be equal to 1. The 

0 \ 





Figure 4.1 Total scattering cross sections as a function of 
T 0 /?i for a -*■ 0 [after Moorovoft , 1963]. 
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Figure 4.3 Spectrum shapes for the ionic component for three ions for 
three values of T e /Ti. These spectra have been plotted 
against a frequency scale normalized with respect to the 
ion mass and electron temperature so that they may be 
compared directly [after Mooravoft, 1964]. It can be seen 
that the shapes arc almost identical in all cases. 










width of the spectrum increeses es increases. This characteristic can 
be used to calculate T ^ from experimental data. 

When only one type of ion is present and T ■ the calculation of 
T^ t T 0 and is straightforward. n g can be found from the power scattered 
per unit voluae o. T. can be calculated from the width of the spectrua (A/) . 
T* is equal to 2\. Generally in the topside ionosphere the ions are hotter 
than the electrons. The quantities a and tf are then functions of the tea* 
perature ratio 2^/2^. However, the height of the peak of the spectrua rela- 
tive to the spectrum value at zero Doppler shift, R ■ max(S(w) }/S(0) , is a 
function only of the temperature ratio. So, if R is measured then can 

be found which in turn can be used in the calculation of T. and It.. Figure 
4.5 shows the dependence of R on 2’ fl /2\ for 0* and 2*(0*) ■ 1000 K. 

4.3.2 Spectra for a mixture of H* and 0* with a single ion temperature. 
The interpretation of the scattering from an ionized gas is considerably 
more complicated when two types of ions are present than when wily one kind 
is present. The discussion in this section is related to a mixture of H* 
and 0*. This is appropriate for the topside atmosphere since densities of 
other ions are very small compared to the combined densities of H* and 0*. 

Moorcroft [1964] has calculated the spectrum for the case where the mass 
of one ion is four times the mass of the other ion. This could apply to mix- 
tures of 0* and He* ions or mixtures of H* and He* ions. Although this anal- 
ysis may be useful where the He* density reaches its maximum value, generally 
the He* density is much too small to have any appreciable effect on the 
incoherent-scatter spectrum. The effect of He* will therefore be ignored 
and the mixture of H* and 0* ions will be studied in detail. 

While the shapes of spectra for two different ions are similar for a 
given temperature ratio, the shape of spectra for a mixture of the two ions 
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li quit* different from those of either ion time. This is evident to Pig* 
ure 4.6 where spectre are shown for various mixtures of 0* and H* for a tea* 
perature ratio of 2.0. This figure also indicates that the peak-to-center 
ratio R is now a function of n(O*)/n 0 as well as of TjT^ which coaplicates 
the calculation of TjT^. The R dependencies oust be represented by a whole 
family of curves rather than a single line as to Figure 4.5. Figure 4.7 
shows £ as a function of n(O*)/n 0 and T^/T^ for mixtures of 0* and H*. 

The appearance of a double hump feature to Figure 4.6 is the result of 
the Doppler frequency separation of the 0* peak and the H* peak. The peak 
of the H* part of the spectrum doesn't show clearly as a peak until n(0 4 )/n # 
is 0.4 or less. Likewise, the 0* peak doesn't become, a peak until n(O*)/n 0 
is 0.6 or greater. Thus, when n(O*)/n 0 is greater than 0.4 and less than 
0.6 there is no peak in the spectrum and the maximum value of the spectrum 
occurs at zero Doppler shift. This makes the determination of TjT^ or 
n(0 *)/n 0 from measurements of R impossible since T /T^ and n(O*)/n 0 can t'- s 
on many values for a single value of R. Thus, an extra parameter is needed 
in order to determine T/T . and n(0 + )/n when n(0 4 )/« is in the range be- 
tween 0.4 and 0.6. When n(Q*)/n is between 0.6 and 1.0 and between 0.0 and 
0.2, R changes rapidly as n(0*)/n g changes. 

In contrast to R, (see Figure 4.8) the width of the spectrum at half 
maximum, changes much more rapidly when n(O*)/n 0 is between 0.2 and 0.6 than 
when n{Q*)/n is between 0.0 and 0.2 and between 0.6 and 1.0. When n(0*)/n 

9 9 

is between 0.8 and 1.0, is almost constant and changes only slightly from 
one value of T/T* to another. This makes the determination of n(0 + )/n and 
T g /T^ from measurements difficult when n(O*)/n 0 is between 0.8 and 1.0., 
with the result that the error in n(0*)n and T/T. would be much larger. 








f 1/2 






107 


One quantity which changes only slightly as the composition is varied 
is afn where a is power scattered per unit volume [Moororoft , 1964], o is 
a function of T which is shown in Figure 4.9. The change in afn Q is 
sufficiently small as the composition varies that it can be either ignored 
or obtained by a simple interpolation between the curves of Figure 4.9. 

Since very little error is introduced by ignoring the dependence of a on 
composition, a is assumed hereafter, to be a function only of and T /T*. 

Generally, the determination of the temperatures and densities of the 
electrons and ions requires at least the measurement of four parameters of 
the spectrum. So, it is clear that • measurement of R , and a alone is in- 

sufficient to determine any of the temperatures or densities. However, prior 
knowledge of any one of T. T. s T /T. t n(0 + )/n or n is sufficient to deter- 

0 o p & 

mine all the rest. 

Sometimes, for example, an independent measurement of the electron den- 
sity is available. Then the ratio o/n via Figure 4.9 will give T /T . which 
combines with R via Figure 4.7 to give «(0 + )/n . Finally, since f, depends 

L 4. 

on T. , T JT . and n(0 )/n combine via Figure 4.8 to give T .. 

An assumption about either or n(0 + )/n g will likewise enable all 

the other quantities to be determined. When T. is known it is necessary to 
make an initial assumption above either n(0 + )/n or T /T . to get the other 

quantity from either Figure 4.7 or 4.8. If the value of n(0 + )/n is assumed 

initially and Figure 4.8 is used to reduce T /T. t then this value of T /T . 

0 IS 0 % 

can be used in Figure 4.7 to determine n(0 + )/n . This iterative procedure 
is continued until n(0 + )/n and T /T. converge. 

C 0 v 

One point which has been ignored in the above discussion is the fact 

that two values of n(0 + )/n are possible for a single value of R and T /T . 

as shown in Figure 4.7. However, this ambiguity can be removed by additional 


J * 
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information concerning the structure of the upper atmosphere. At lower alti- 
tudes it is known that the amount of H* is very small. As the investigation 
is extended upward from the Fj layer only (me of the two possible sets of 
values for n(0 + )/n gives a continuous variation of n(0*)/n with height. 

w 9 

If information about n_, TJT. , T or T. is not available and if one 
is to deduce information about temperatures and composition from scattering 
measurements alone, then another independent measurement on the spectrum is 
necessary. A survey of the shape of spectra in Figure 4.6 shows that the 
tail of the spectrum becomes more extended as n(0*)/« decreases. It appears 
that the measurement of the slope of the tail of the spectrum will yield in- 
formation on «(0 + )/n^ and T 0 /r.. The choice for this new parameter is the 
slope of the spectrum at the point where the power spectral density falls to 
one-half its central value and it will be called Sy In order to keep it in- 
dependent of the scales of power density and Doppler shift, can be meas- 
ured on a semi logarithmic plot of the spectrum in which the central value of 
the power spectral density is normalized to zero. Because of this normaliza- 
tion, is then expressed as the number of decibels of change over a fre- 
quency interval corresponding to the Doppler shift to a point - 3 dB from 
the central value of the spectrum. 5^ is independent of the electron and ion 
temperatures and is a function only of T^/2\ and ^(Q*)/^, 

4.3,3 Spectra for a mixture of H* and 0* itith different ion tempera- 
tures. The influence of differing ion temperatures on the spectrum has been 
ignored in the previous discussion. The ion temperature also has been as- 
sumed to be the same for all the ions present. This assumption is valid at 
lower heights where the ions, electrons, and neutrals all have the same tem- 
perature. When T JT. is found to be equal to one using the analysis of the 
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previous secticn one can safbly assume that the Ions have a single tempera- 
ture. But as TjT i Increases above one then the possibility that the Ions 
have separate temperatures must be Investigated because H* and 0* Ions are 
heated at different rates by the thermal electrons. In this section the 
effect of separate ion temperatures on spectra for a mixture of H* and 0* 
are investigated in detail. 

In general the parameters, a , and are functions of each ion 

temperature. In any of these parameters are functions of T. then it is a 

v 

function of both the hydrogen ion temperature (T(H*)) and the oxygen ion 
temperature (T(0*)). The only parameter which was a function of T. was /.. 
Therefore, one expects / ^ to be a function o5 both 1(0*) and T(H*). The 
other parameters may also be functions of T(H*) even though they were not 

functions of T.. This will be discussed later in this section. 

% 

Instead of dealing with both T(0*) and 7(H*) the ratio T(H*)/T(0*) and 
T( 0*) will be used. This simplifies some of the numerical calculations which 
will be carried out later on. In Figure 4.10 the spectra for values of 
n(0*)/n equal to 0.0 to 1.0 are shown for different values of the ratio 
r(H*)/r(0*) . In this figure the ratio T/T. is taken as 2.0. Notice in this 
figure that a point near the peak due to H* ions the curves converge to a sin- 
gle value. For all values of n(0*)/n the power at zero Doppler frequency 
is greater for larger values of T(H*)/T(0*). Also the power decreases less 
rapidly after the H* peak for larger values of T(H*)/r(0*). The shape of 
the spectrum is different for each value of n(0*)/n fi and r(H*)/!T(0*) . There 
are points where the spectra intersect which provide a way of calculating 
n(0*)/n fl . The ratio T(H*)/T(0*) can be calculated from places where the 
spectra differ. Unfortunately most of these parameters also depend on other 

quantities such as 7(0*), T 0*), and n(0*)/n . It is very difficult to 

6 & 


RELATIVE POWER RELATiVE POWER 




{ Figure 4.10 Spectra for different values of n(0 + )/« fl and 

f different values of r(H 1/T(0 + ) where T e /T{ 0 I * * 4 ") 

| (a)- «(0*)/« e - 0, (b)- n(0 + )/«. - 0.3, 

| (c)- u(0*)/« c *• 0,6, (d)- «(0 + )/n^ * 0,9. 
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find parameters of the power spectra which only depend on one quantity such 
as «(0 + )/n . This makes the determination of T(0 + ), T JT (0 + ), !T(H + )/r(0 + ) , 
and n(0 + )/n very complicated. Sometimes the solution is ambiguous due to 
two or more * ys of matching theoretical results with experimental data. 

The parameter which was a function of T., is a function of both ^(O*) 
and !T(H + )/r(0 + ) . As can be seen from Figure 4.10, increases as r(H + )/T(0 + ) 
increases. Another parameter which increases with r(H + )/rCO + ) is o. For 
larger values of J’(H + )/2'(0 + ) amplitude of the spectrum at frequencies near 
zero Doppler shift is also larger. Earlier it was found that as !T e /r(0 + ) in- 
creases, o decreases, so r(H + )/!T(0 + ) has the opposite effect on a as does 
T e /T(0 + ). 

The effect of T(H + )/T(0 + ) on R is not readily apparent from Figure 4.10. 
Both the spectrum maximum and the spectrum at zero Doppler shift changes as 
changes. If one compares the values of R for n(0 + )/n = -/8 
at different values of T(H + )/T(0 + ) one can get the relation for the dependence 
of R on T(H + )/r(0 + ) at « 1, R lQ = S 1 /S Q and at r(H + )/T(0 + ) = 1.5, 

R 1S - (Sj + LS)/{S q + A S) so 


“ 1 - AS 


s i - 5 o 

TV 


Thus, as T(H + )/r(0 + ) increases from 1.0 to 1.5, R decreases slightly. 

4.4 Characteristics of Theoretical Autocorrelation Functions 

Although the analysis of various par inters of the incoherent-scatter 
power spectrum is very useful in determining temperatures and densities of 
electrons and ions, this analysis cannot be directly applied to experimental 
data. Instead, the autocorrelation function (ACF) or the Fourier transform 
of the power spectrum is available. To get the experimental power spectrum, 
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the Inverse Fourier transform of the experimental ACF has to be calculated. 
This leads to larger errors in the determination of temperatures and densi- 
ties of electrons and ions. Fortunately, analysis of theoretical ACF's can 
be performed to determine temperatures and densities of electrons and ions. 

The method of this analysis is discussed in the following subsections. 

The ACF like the power spectrum depends on r(0 + ), TjT{0*) t f(H + )/r(0 + ) 
n(0 + )/n . By identifying certain characteristics of the ACF one can under- 
stand how the ACF can be used to deduce T(0*), T /T(0*) t !T(H + )/!r(0 + ) and 
n(0*)/n , Some of these quantities affect the ACF significantly. Others 
have a slight effect on the ACF. Parameters which are derived from the ACF 
are discussed in the following subsections. These parameters can be used 
to determine the temperatures and densities of electrons and ions. In 
Chapter 5 analysis of the theoretical ACF will be applied to experimental 
ACF's to calculate the temperatures and densities of electron and ions. 

4.4.1 ACF for a single ion. The analysis of the ACF when a single ion 
is present will serve as a basis for the more complicated case of a mixture 
of ions which will be discussed later. The calculation of the ACF has been 
performed numerically using Simpson's rule for integration. The time delay 
step size has been chosen to be at least as short as time delay of experi- 
mental ACF's. The general features of the ACF are illustrated in Figure 4.11, 
which shows ACF's for scatter from ionized gases where the ion present is 
oxygen or hydrogen and for electron-ion temperature ratios of 1.0, 1.4, and 
2.4. All the ACF's have been normalized to one at zero time delay. The posi- 
tion of the first minimum of the ACF, Tj, is smaller for H + than for 0 + , and 
is nearly inversely proportional to the square root of the ion mass. Since 
the shape of the ACF is of particular interest, the abscissa can be divided 
by the square root of the ion atomic mass shown in Figure 4.12. It is clear 
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T (4* lO^sec) 

Figure 4.12 ACF's for 0 + and H + for T e /Tf ■ 1.0, 1.4, 

and 2.4 with abscissa normalized by the | 

atonic mass of the ion present (AO. j 
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from this figure that the shape of the ACT is only slightly dependent on the 
nass of the ion present in the gas. If there were no contribution to the 
spectrua due to electrons in equation (4.2) then the shape of the spectrua 
and the ACF could be independent of the ion Bass. 

Figure 4.11 also shows that the depth of the first ainiaua of the ACF, 
AMIN, is strongly dependent of the electron-to-ion temperature ratio. This 
characteristic of the ACF is due to the dependence of peak-to-center ratio R 
of the spectrum on . As R increases for increasing so does the 

depth of the first minimum of the ACF. To confirm the relationship between 
R and AMIN one would except that they do not depend on the ion temperature. 
Analysis of Figure 4.13. which shows ACF's for different values of T indi- 
cates that the AMIN does not depend on T ^ even though the position of AMIN 
does. Thus, if only one type of ion is present in the gas the electron-to- 
ion temperature ratio could be determined from AMIN. 

Another feature of the ACF which is shown in Figure 4.13 is the depend- 
ence of the first zero crossing of the ACF (t^) on the ion temperature. As 
the ion temperature increases, t q decreases. This is opposite to the depend- 
ence of / ^ on which was discussed in Section 4.2. It was found that / ^ 
depends on 2^ . Here, the dependence of t q on T ^ is t q • 3^ . The first 
zero crossing doesn't depend on T. alone, as can be seen from Figure 4.11 
which shows ACFs for different values of T /T * . It is evident from this fig- 

p X* 

ure that t q also depends heavily on T g /f ^ . As T^/7^ increases, t q decreases. 
So, the functional dependence of is complicated since it both depends on 
and T e /T v 

When only one type of ion is present one can use the dependencies of 
o, t q , and AMIN on an< * n e t0 deduce these quantities. These 
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^ Figure 4.13 ACF's for 0* for different values o 
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various dependencies are illustrated in Figures 4.14 and 4.15 for the ions 
0* and H*. In Figure 4.14 the ACF at its first minimum, AMIN, is shown as 
function of T for temperature ratios froai 1.0 to 2.4, the appropriate 
range of for the upper atmosphere . Since the ion mss has a slight 

effect on the shape of the ACF (see Figure 4.12) AMIN depends just slightly on 
the ion mss. The effect of T. on AMIN is slight and can be Ignored compared 
to the dependence of AMIN on f^/ST.. It is even less significant than the 
effect of changes in the ion mss. Thus, the functional form of AMIN can be 
expressed in terms of a 1 one. 

Both Tq and o depend on but can be determined from AMIN. 

Thus, one can find T . from t- and n from a by using T/T. from AMIN. Al- 
though o is only slightly dependent on the ion mss, t q is strongly affected 
by the ion mss. Therefore, no further complications are introduced by the 
ion mss when using a to find n^. But, t c cannot be used alone to find T 
Either some information about the ion present in the gas must be available or 
t q must he normalized with respect to some other parameter. Sometimes the 
ion present :n the gas is known. For example, at altitudes near the P ea * 
in the electron density it is known that the ionic composition is almost ex- 
clusively 0*. At other altitudes the knowledge of the ion present is less 
certain. In fact, at most altitudes there is a mixture of ions with no one 
ion more than 90 percent of the total ion density. 

Suppose there is only one ion present and it is not known. must thsn 
be norm li zed with respect to some other parameter. A suitable parameter to 
use is t q ; which is a value of t q at a fixed ion temperature of 1000 K. This 
ion temperature is chosen because the ion temperatures in the topside iono- 

t 

sphere are near 1000 K. In Figure 4.15 the normalized Tq is shown as a func- 

4 + t 

tion of for 0 and H . Although depends slightly on the irn mass 


AMIN 



Figure 4.14 The first minimum of the ACF is called AMN. It is 
shown for 0* as a function of T e /T^(.0*) and T{(0*). 
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one can determine T. fairly accurately if TJf* is known from analysis of o. 

9 

If only one type of ion is present and if AMIN, t q , and a can be mea- 
sured then it is possible to determine T ^ and n # using the curves of Fig- 

ures 4.14, 4. IS and 4.9. However if not all these scattering quantities can 
be measured, it may be necessary to make some assumptions about the tempera- 
tures, or perhaps measure in another way. In Figure 4.16 the various com- 
binations of measurements and assumptions are shown by a schematic diagram. 

The case where not all of a, AMIN, and can be measured can be con- 
sidered with the aid of Figure 4.16. Of the three quantities, the most dif- 
ficult to measure is AMIN. Examination of Figure 4.16 shows that the measure- 
ment of Tq and a must be supplemented by a measurement or an assumed value 
of one of the other quantities before anything can be deduced. Any of T., 

T , T /T., or n will suffice. At certain times, the temperature ratio in 
the topside ionosphere can be accurately estimated. For example, Z^/2^ ■ 1 
at night; also, the ion temperature could be obtained from rocket measure- 
ments if an isothermal atmosphere is assumed. Either of these assumptions 
would be sufficient to deduce both temperatures and the electron density 

I 

from measurements of o and Tq. For example, values of the electron density 
from topside sounders could oe combined with o to give T /T.. This value of 

© I* 

I 

T /T. along with analysis of T n would give T. and T . 

£ t U 

Two quantities which are necessary for the determination of temperatures 

I 

and the electron density are t q and o. Without them it is not possible to 
determine temperatures and the electron density. Early measurement of Bowles 
[1962] only contained measurements of o. In that case, assumptions about 

T /T . will give n from o; or an independent measurement of n will give 

# % # 0 

t 

T /T,, When only x rt is measured, one of T .. T , or T / T , must he obtained 

e t 0 t e e t 
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Figure 4.16 Schematic diagram of interrelations between scattering 
observations and characteristics of an ionised gas when 
only one type of ion is present. 
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by independent measurement or assumption and then the others can be obtained 

i 

from Tq. 

4.4.2 ACF for a mixture of two ions at a single ion temperature . The 
interpretation of the scattering from an ionized gas is considerably more 
complicated when two types of ions are present than when only one kind is 
present. The discussion in this section is related to a mixture of H* and 
0 + , the ions present in the topside ionosphere. 

While the shape of the ACF for two different ions are similar for a 
given temperature ratio, the shap* of the ACF for a mixture of the two ions 
is quite different from those of either ion alone. This is evident in Figure 
4.17 which shows ACFs for various mixtures of 0 + and H + for a temperature 
ratio of 1.4 and an oxygen ion temperature of 1000 K. This figure also in- 
dicates that the depth of first minimum of the ACF, AMIN, is now a function 
of n(0*)/n g as well as of T g /T^. This complicates the calculation of T g /T^ 

and n(0*)/n . The functional form of AMIN can be represented by a whole 
0 

family of curves rather than a single line as in Figure 4.14. Figure 4.18 

shows AMIN as a function of n(0 + )/n and T/T. for mixtures of 0* and H + . 

6 6 1 

There are a number of characteristics in Figure 4.18 which either hinder 

or aid the determination of T/T. and n(0*)/n . One characteristic which 

0 % 0 

aids the determination of T /T. and n(0 + )/n is the fact that at either small 

0 % 0 

or large values of n(0*)/rt , T/T. can be determined fairly accurately if 

n(0 + )/« e is known from other measurements. In opposite sense if T^/T. is 

known, Figure 4.17 will give n(0*)/n g if it is small or large. 

One characteristic which hinders the determination of T/T. and 

s x 

n( 0*)/n is the fact the curve of one value of T/T. either crosses the other 

0 0 1 

curves of constant T /T . or becomes almost equal to the curves of constant 

T/T. at a point in the range between n(0*)/n equal to 50 and 70 percent. 
ex s 
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In fact the curve for T * 1,0 crosses all the other curves twice between 
«(0 + )/n c equal .4 and .7. This characteristic means that 2 values of T^/T^ 
can be deduced from AMIN at these points of crossing. Also, AMIN cannot be 
used to find n(0 + )/n or T/T, when n(0 + )/n is between .4 and ,7 since the 
curves cross and reach a maximum. This implies, with the use of Fig. 4.18, 
that when AMIN is greater than -.1 some other parameter must be used to find 
n(0*)/n and T/T.. 

0 e 1 

Another characteristic in Figure 4.18 is the similarity of values of 
AMIN for n(0 + )/n near 0.0 or 1.0. This means two values of n(0 + )/n can 
be determined from a single AMIN using one value of T / T .. Usually this is 
not a problem, since the dominate ion is often known in a particular alti- 
tude region and n(0*)/n can be assumed to be smoothly varying with altitude. 
Also, the ACF is quite different for n(0 + )/n near 0.0 and 1.0 (see Figure 
4.17). 

An important parameter of the ACF which depends on n(0 + )/n , T /T ,, and 
T. is the time delay at the point where the ACF is one half (THALF) . Figure 
4.19 shows the dependencies of THALF on n(0 + )/n , T /T and T.. The depend- 
ence of THALF on n(0 + )n is more pronounced with THALF much less for n(0 + )n 
equal to 0 than 1.0. THALF depends on T. more than it depends on T /T,. 

Yet, the dependence of THALF on T /T. is significant and cannot be ignored. 

0 % 

Since THALF depends on all the quantities that need to be deduced, one must 
combine measurement of THALF with measurement of at least two other para- 
meters such as AMIN, to find n(0 + )n . T/T., and T 

S0% % 

Another important parameter of the ACF is the time delay at the first 
zero of the ACF (TZERO). TZERO depends on T., T_/T., and n(0*)/n as shown 
in Figure 4.20. In this figure, TZERO has been normalized by THALF. Hence- 
forth, the normalized TZERO will be referred to as TZERO. 


THALF 







Figure 4.20 TZERO as a function of w(0 )/n^» 

The thin lines are for Ts • 1000 K. The thi 
lines are for T{ ■ 2000 K. The dotted lines 
indicate discontinuities in TZERO. 
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In Figure 4.20, the dependence of TZERO on n( it noticed at the 
strongest dependence and the dependence of TZERO on as the next strongest . 
The dependence of TZERO on is very slight when n(0 + )/n # is greater 

than 0.60 and less than 0.20. When n(0*)/n tf is between 0.2 and 0.6 there 
is some separation between the curves for different values of The 

curves are discontinuous when n(0*)/n Is between 0.2 and 0.6. The reason 
for this has to do with the disappearance of the minimal of the ACF due to 
H* ions as n(0 + )/n tf increases. The minimum becomes less in magnitude as 
n(0 + )/n increases and increases above zero. Then the minimum due to H + 
disappears altogether and the first minimum is due to 0 + . The minimum changes 
discontinuous ly when this occurs. The place where the curves become disconti- 
nuous depends on the ratio T g /7\ as well as on n(O*)/n 0 . Since the respec- 
tive step sizes for n(0 + )/« e and were 0.1 and 0.2 when these calcula- 

tions were made, the resolution of the curves in the range of n(0*)/n be- 
tween 0.2 and 0.6 is not very good. The resolution could be improved with 
smaller step sizes for either nC0 + )/n tf , or T^/3^ , or both. 

In Figure 4.21, TZERO is shown for n (0 + ) /n g in the range 0.2 to 0.6 with 
step sizes of 0.04 for n(0 + )/» . This figure shows that the curves have bet- 
ter resolution when the step size for n(0*)/n is smaller. Although Figures 

9 

4.20 and 4.21 could theoretically determine n(0*)/n very accurately when 

G 

TZERO has a discontinuity in its experimental measurement, this is dependent 
on the accurate measurement of T/T. and T .. Even though T* and T /? . are 

not accurate. Figure 4.20 still could be used to get a fairly accurate value 

* 

of n(0 )/n g if TZERO is greater than 2. When TZERO is less than 2 and a dis- 
continuity has just occurred n(0 + )/n is between 0 and .54. If TZERO is less 

0 

than 2 and n(0*)/n is near 1 then TZERO can again be used to find n(0 + )n 

0 0 
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If a single ion temperature is used for both ions and there is a Mixture 

of 0* and H + then the parameters TZERO, THALF, and AMIN (Figures 4.18-4.21) 

can be used to find T^, and and o can be used to find n g . 

If o is not measured, cannot be found. If measurement of one of the other 

a 

parameters is not possible then the independent knowledge of one of T^, 

or n(0*)/n mist be had in order to determine the rest of these quantities. 

4,4.3 ACF for a mixture of 0 + and H + with different ion temperaturee. 
This section covers the special case of deducing all the possible quantities 
of a mixture of 0* and H* from theoretical parameters of the ACF. That is, 
n(0 *)/«„, T(0 + ), !T(H + )/r(0 + ), T/TQi*) and n can all be calculated with suit- 
able parameters of the ACF. This is possible theoretically since the power 
spectrum and in turn the ACF depend on all these quantities. The dependence 
of the ACF on T(0*) if T. • 7(H*) - T(0 4 ) has been shown previously in Figure 
4.13. The dependence of the ACF on T /T. and n(0 + )/n has been shown in 
Figures 4.11 and 4.17 respectively. The dependence of the ACF on 2r(H + )/r(0*) 
is shown in Figure 4.22 for the cases n(0*)/n g ■ 0.4 and 0.9 with T g /T( H + ) 
equal to 1.4. These are only two of the many possible cases. Other cases 
will be discussed and illustrated in this section. 

The dependence of the ACF on r(H*)/T(0*) as shown in Figure 4.22 is weak 
but still strong enough to permit determination of r(H*)/r(0*). The para- 
meters THALF, AMIN and TZERO which depend on «(0 *)/m , T/T. , and T. may now 

m e i» w 

depend on T(H*)/r(0*) also. From Figure 4.22 one notices that THALF depends 
on TCH^/rfO*). Since THALF also depends on n(0*)/n g , TJT., and T ^ it is 
hard to represent all the dependencies of THALF on a single graph. THALF is 
shown in Figure 4.19 as a function of n(0 4 )/« # , 3^/^. and in Figure 4.23 
THALF is shown as a function of n(0*)/n g , TW^/TCO 4 ) , and T(0 4 ) with 
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7_/7( H + ) ■ 1.0. Examination of Figure 4.23 shows that the dependence of 
THALF on 7(H + )/7(0 + ) weakens as 7(0*) increases since the curves are closer 
together for the case where 7(0*) ■ 2000 K case than for the 1000 K case. 

Even though the dependence of the ACF and THALF on 7(H + )/7(0*) (Figures 4.22 
and 4.23) indicates that it is theoretically possible to deduce ratios of 
7(H + )/7(0*) if the other quantities that the ACF and THALF depend on are 
known. 

The parameter TZER0 which depends on n(0*)/n and 7. does not depend 
on 7(H + )/7(0 + ) in any significant way when «(0*)/« is greater than 0.6. The 
only places where THALF depends on 7(H + )/7(0*) are at the discontinuities. 
Since THALF does not depend on 7(H + )/7(0*) for most of the range of n(0 + )/n , 
the dependence of THALF on 7(H + )/7(0 + ) will be ignored, just as the dependence 
of THALF on 7/7. is ignored for most of the range of n(0 + )/n . 

Another parameter which depends on 7(H + )/7(0 + ) is AMIN. In Figure 4.24 
AMIN is shown as a function of 7(H + )/7(0 + ) and n (0 + ) /« for 7 /7(H + ) ■ 1.0 
and 7(0 + ) = 1000 K. Fortunately, the dependence of AMIN on 7(0 + ) is very 
slight; so this figure can be used for all values of 7(0 + ) . The dependence 
of AMIN on 7(H + )/7(0 + ) is strongest, for values of n(0 + )/n near 1.0. When 
n(0 + )/n decreases below 0.5 there are values of AMIN which give 2 values 
of 7(H*)/7(0 + ). This is because the curves cross one smother. Thus, other 
parameters must be used to deduce 7(H + )/7(0 + ) when n(0 + )/n is less than 0.5. 
The dependence of AMIN on 7(H + )/7(0 + ) is also somewhat similar to the depend- 
ence of AMIN on 7 /7(H + ) as shown in Figure 4.18. This could cause a problem 
with ambiguous solutions, i.e. two sets of values for 7 /7(H + ) and 7(H + )/7(0*) 
could give the same AMIN. This problem can be eliminated if another para- 
meter can be found which depends on 7 /7(H + ) and 7(H + )/7(0 + ) in much different 
ways than AMIN does. 
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A parameter which depends on T(H + )/TCO*) is ADHALF, which is the deriva- 
tive of ACF at t * THALF. ADHALF also depends on some of the other quanti- 
ties. In Figure 4.25 ADHALF is shown as a function of n(0 + )/n g and T(H + )/ 

r(0 + ) with T /7(H + ) - 1.0 and 7(0 + ) ■ 1000 K. The dependence of ADHALF on 

6 

7(H + )/T(0 + ) is slight compared to the dependence on «(0 + )/n . Usually 

& 

n(0*)/n is determined quite accurately from TZERO, so the dependence of 
ADHALF on 7(H + )/7(0 + ) can be used to determine T(H + )/T(0 + ) . This is not very 
easy because ADHALF also depends on T /T (H + ) . Although the dependence of 
ADHALF of T /T(H + ) is not similar to the dependence on T(H + )/T(0 + ) , as in the 
case of AMIN, it is close enough to cause some problems with ambiguous re- 
sults. Usually this can be overcome with some assumptions based on physical 
reasoning. For example the quantities T /T(H + ) and 7(H + )/T(0 + } are expected 
to vary smoothly with altitude. Since many altitudes are observed, results 
that are discontinuous will be rejected. 

There are many other parameters of the ACF which depend on the quantities 
7(0 + ), n(0 + )/n , T /T(H + ), and 7(H + )/7(0 + ) . Parameters such as the time de- 

o & 

lay at the first minimum of the ACF (TMIN) and the derivative of the ACF at 
the first zero (Aq) could be used to reduce these quantities. Like the 
other parameters discussed previously, TMIN and Aq have complicated de- 
pendencies. They will not be discussed here but they may be used if another 
parameter is not useful in reducing the ion and electron temperatures. 

4.4.4 Surmary. There are five quantities to be determined from analysis 
of the ACF - n , n(0 + )/n , T( 0 + ) , T/T( H + ) , and r(H + )/T(0 + ) . The only quan- 

e & o 

tity which is easy to determine is n(0 + )/n if it is greater than 0.5. It 

can be determined quite accurately from TZERO. The quantity n can be deter- 

mined to within a factor of 2 from o. If T /T CH + ) is known, then n can be 

6 €■ 

determined much more accurately. The 7(0 + ), T /T(H + ), and 7(H + )/7(0 + ) are 

6 
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much more difficult to determine because the ACF parameters depend on all 
these quantities. If either T/T CH + ) or 2 , {H + )/2 , CO + ) is known then the task 
of finding the other quantities is made easier because one dependence is re- 
moved. If all the quantities are unknown then the ACF parameters in conjunc 
tion with physical assumptions have to be used to reduce these quantities. 
Sometimes even this doesn't work, in which case the data is rejected either 
because of an experimental failure causing the problem or the occurrence of 
other physical events which are not of interest in this study. It was found 
that useful information can be reduced from almost all of the data. 
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5. METHOD OF ANALYSIS AND EXPERIMENTAL RESULTS 
5.1 Introduction 

The preceding chapter dealt with different parameters that can be derived 
from theoretical incoherent-scatter autocorrelation functions and spectra. 

In this chapter the method of fitting these theoretical parameters to the 
corresponding experimental parameters will be described. From this one can 
deduce electron, 0 + and H + ion temperatures, and the relative concentrations 
of 0 + and H + ions. The only ions present in the scattering medium are as- 
sumed to be H + and 0 + . 

The radar signals from the Arecibo, Puerto Rico facility are in the 
form of long (2000 ys, 300 km) and short (300 ys, 45 km) pulses transmitted 
alternately in order to obtain both adequate sensitivity at the highest alti- 
tudes and sufficient height resolution for the lower altitudes in the F re- 
gion. The short pulse provides 48 km resolution over the range 230 to 902 km 
while the long pulse provides 108 resolution over the range 340 to 1850 km. 
That is, the short pulse samples the ACF every 48 km and these ACF's are 
averaged over 45 km. The long pulse samples the ACF every 108 km and these 
ACF's are averaged over 300 km. 

The basic design of the Arecibo facility has been described previously 
by Gordon [1964], The transmitter operates at 430 MHz with a peak power of 

2.5 MW. For this work the beam was directed at the zenith. The received 
430-MHz signal is prc messed by converting it to an IF of 30 MHz and passing 
this through a 125 kHz wide filter and then converting to a second IF of 

62.5 MHz which was finally sampled and applied to a digital autocorrelation 
[Hagen and Heu t 1974]. The hardware correlator has been described by Hagen 
and Farley [1973], This hardware correlator forms the autocorrelation 


functions in real time. 
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The autocorrelation functions are sets of 24 points for each of the IS 
altitudes sampled. The time spacing between adjacent points is 8 us. This 
data is stored on magnetic tape so that additional processing can be done at 
a later time. The ACF's on tape are 12 minute averaged ACF's. That is, 
every 12 minutes the integrated ACF's for both the short and long pulses are 
written on magnetic tape. 

5 . 2 Methods of analysis of the ACF 

Two methods of analysis of the ACF are described in this section. The 
first method is the method of fitting ACF parameters using Newton's method 
of solving nonlinear equations. The second method is the method of fitting 
the ACF's using nonlinear least squares analysis; in the second method the 
entire ACF is used. 

5.2.1 Fitting ACF parameters . Information about the temperature and 
composition of the charge particles of the F region can be obtained by 
matching experimental ACF parameters with theoretically generated ACF para- 
meters. The parameters THALF, TZERO, AMIN, and ■ ■ * ADHALF are 

first calculated from the experimental ACF. Then these values can be com- 
pared with values calculated from theoretically generated ACF's which has 
been discussed previously in Chapter 4. As discussed in Chapter 4, these 
parameters are not simple functions of T(0 + ) , r(H + )/T(0 + ) , T /T (0 + ) , and 
n(0*)/n . So in order to derive temperatures and composition from THALF, 
TZERO, AMIN, and ADHALF one must first make a initial assumption about the 
temperatures and composition. Then n[0*)/n is calculated from TZERO. This 
value of n(0 + )/n is then used to calculate T fT{ 0 + ) from AMIN. These values 

0 c 

of T /T(0 *) and n(0*)/n g are then used to calculate T(0 + ) from THALF. !T(H + )/ 
r(0 + ) is then calculated from ADHALF using the previously calculated values 
of T{ 0*) , T /T(0*) and «(0 + )/»t . This procedure is repeated with these new 
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values to solve for n(0 + )/n_, T. /r(0 + ), T(0 + ), and TOH^l/rCO*) until they 

c w 

converge. This procedure is illustrated in Figure 5.1. This is the general 
method for solving for n(0 + )/« . TJT(0*) t T( 0 + ) and T(H*)/T(0*). Sometimes 
this method has to be altered due to problems in the data. 

The numerical method used in deducing the ion and electron temperatures 
and number densities was Newton's method for solving nonlinear equations. 

The equation that was solved is 

/(*) * V e - y t («) * 0 ( 5 -l) 

where y is the experimental ACF parameter, y. is the theoretical ACF para- 
meter and x is the physical quantity one is trying to find. The solution 
for x using Newton's method is found by generating successive approximations 
from the iteration 

x i+l * x i - - 9(pj) (5.2) 

where /' (x») is the derivative of j/.(x). 

The expected error in x is 

Xfc+i - * a *4 - x) 2 (5.3) 

where is between x^ and x. 

One of the main problems is the obvious fact that when n(0 + )/n * 1.0 
no information about the H + temperature can be deduced from the ACF. Even 
though n{Q*)/n e is not often 1.0 it is close enough to 1.0 at altitudes 
around F 2 layer peak so that H + temperatures are almost impossible to deduce 
from the ACF alone. In this case the H + temperature can be derived from the 
shape of the 0 + and electron temperature profiles using results from therma* 
processes analysis i.e. the 0 + and electron temperature profiles can have 
different shapes depending on whether the H* temperature is equal to or sepa- 
rate from the 0 + temperature. This will be discussed in more detail in 


1. Assume Initial Values 


2. Calculate n(0*)/n from 
TZERO 


3. Calculate T fT{ H + ) from 
AMIN e 


4. Calculate TfO*) from 
THALF 


5. Calculate 3’(H + )/2’(0 + ) 
from ADHALF 


6. Repeat 2-5 until all 
values converge 


Figure 5.1 Schematic diagram of the procedure for 

finding T ( 0 ), T(H + ), T , and n(0*)/n. 

& & 
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Section 5.4. 

The deduction of n(0 + )/n from TZERO is sometimes ambiguous due to the 
discontinuity in the theoretical value of TZERO at different values of 
n(0*)/n e for different T^T (0 + ). This was discussed in Chapter 4. Even 
though the experimental data of TZERO has discontinuities also, it is diffi- 
cult to match the theroetical TZERO with the experimental TZERO at the alti- 
tudes where the discontinuity takes place. There are two possible solutions 
to this problem: the first is to try to fit theoretical and experimental 

values of TZERO but working from the other side of the discontinuity; the 
second is to find another parameter from the ACF which depends on n(0 + )/n g in 
a different way than TZERO does. This will be discussed in detail later in 
Section 5.4. 

The parameter used to calculate !T(0 + ) is THALF. THALF depends on !T(0 + ) 
in a simple way but THALF also depends on «(0 + )/« T(H + )/f(0 + ) , and T /T(0*). 

This complicates matters a little but if n(0 + )/n is known !T(0 + ) can be cal- 

6 

culated with reasonable accuracy (200 K or 20%) since the dependence of THALF 
on ?(0 + ) is stronger than the dependence of THALF on T fT{ H + ) . The calcula- 
tion of TfO*) will be discussed in Section 5.5. 

Solving for T fT (0 + ) is simplified if one first ignores the dependence 
of AMIN or. T(H' , ')/!r(0 + ) . Then AMIN only depends on T/T{0*) and n(0 + )/n ; it 
doesn't depend on T( 0 + ) . Now if n(0 + )/w is assumed to be known as it was 
for finding T( 0 + ) from THALF then T JT{ 0 + ) can be found to within 20%. This 
accuracy can be improved and this will be discussed in more detail in Section 
5.5. 

There are two more quantities which need to be found n g and 7'(H + )/f(0 + ) . 
The electron density n g can be found quite easily if the total scattered 
power is known. Data from Arecibo gives the scattered power very accurately 
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and thus the electron density can be calculated quite accurately. The de- 
termination of T(H*)/r(0 + ) is not so precise. This is because the ACF de- 
pends on 7(H 4 ')/T(0^) only weakly. Yet there is a definite dependence in the 
shape of the ACF. Some of the parameters of the ACF do depend on r(H + )/r(0 + ) 
and these will give approximately values for T(H*)/?(0*) . Extensive analysis 
of the ACF which will be discussed in Section 5.S give better accuracy foT 
T(H + )/r(0 + ). 

5.2.2 Fitting ACF’8. Experimental ACF's can be matched with theoret- 
ical ACF's using the usual least square criterion. The 24 experimental 
points are each weighted by the reciprocal of its expected variance p n< 

The physical quantity x is found from locating the minimum of 
24 , 

/(*) - I [(C n - V*))/p„l (5,4) 

n»l 

where G is the theoretical ACF. The condition for the minimum is simply 

3,f/3x ■ 0 (5.5) 

Since the theoretical ACF's are not simple functions of x, numerical tech- 
niques are used to find the best fit. The solution for x using nonlinear 
least squares analysis is found by first numerically differentiating fix). 
Then successive approximations are generated by using 

x k+1 . x k - fix k )/f'(x k ) (5.6) 

If x represents a set of physical quantities then the second terra on the 
right side of equation (5.6) is the solution to the matrix equation 

A&x^ ■ B (5.7) 

where A..^ * (3//3x,.) (3/*/3x^) 

B i - (3 f/9x € )f 
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Equation (5.6) then becomes 
■**♦1 ^ +k 

X ■ X * X 

The variance of the final fitted quantities is given by 
variance (x^) ■ 


(5.8) 


(5.9) 


The matrix is inverted in the course of solving equation (S.7). 

The physical quantities which can be used for x are r(0 + ), T g /T(0*) $ 
T(H + )/3’(0 + ), and n(0*)/n,. So matrix A is a 4 x 4 matrix. Particular care 
must be used in inverting this matrix because sane of the elements of A are 
very small compared to other elements. This gives a determinant of A which 
is close to zero. If the determinant of A is in fact zero then the inverse 
of A does not exist and matrix equation (5.7) has no unique solution. If the 
determinant of A is close to zero then the solution of matrix equation (5.7) 
is very difficult using nonlinear least squares analysis. In this case, 
Newton's method can be used to fit the theoretical ACF to the experimental 
ACF. 

By fitting the entire theoretical ACF to the experimental ACF better 
accuracy can be achieved when deducing the ratio r(H 4 )/?(0*) and also when 
deducing 7'(0 + ) and T/T(0*). The ratio n{0*)/n is not deduced by this 

P P 

method because it can be deduced accurately from TZER0. 

The method of fitting the entire ACF is a more accurate method of de- 
ducing the physical quantities than fitting parameters of the ACF. The rea- 
son is that much more information is available in the ACF than in the ACF 
parameters discussed in Chapter 4. All the physical quantities could be de- 
duced by this method but n(0*)/n # and sometimes the other quantities are cal- 
culated by fitting ACF parameters because fitting the entire ACF is computa- 
tionally slower and requires much more computer resources. The ratio n(0*)/w 
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and some of the other quantities are deduced by fitting ACF parameters when- 
ever the accuracy is good. 

5.3 Ion and Electron Number Densities 

The determination of n(0*)/n g and n g will give the ion and electron 
number densities since n(0 + ) ♦ n(H*) ■ n g in the topside ionosphere. From 
TZERO one can calculate n(0*)/n g if 7(0 + ) is known. Even if T(0*) is not 
biown, n(0*)/n g can be estimated from assuming a typical value for 7(0*). 

For example, if one assumes that 7( 0*) is 1000 K throughout the topside 
ionosphere, one gets what is shown in Figure 5.2 for n(0*)/n g profiles. Also, 
shown in Figure 5.2 is the n(0*)/« profile for 7(0*) ■ 2000 K. Comparison 
between the 7(0*) ■ 1000 K and 7(0*) ■ 2000 K profile for n(0*)/n shows that 
error is not too large in the measurement of n(0*)/n if 7(0*) is not known. 
Thus, reasonable estimates of n(0*)/n can be made even if 7(0*) is not known. 

Profiles of «(0*)/n can also be calculated using different assumptions 
on T g /T^. This is shown in Figure 5.3 with T /2\ equal to 1.0 and 2.0. As 
one can obviously see from this figure the effect of 7 JT, on n(0*)n is very 
small except when n(0*)/n is around 0.5. When n(0*)/n is greater than 0.6 
the dependence of n(0*)/n on 7 JT. can be ignored. The reason for disagree- 
ment in n(0*)/n g measurements around 0.5 is due to the discontinuity in theo- 
retical values of TZERO near n(0*)/n g - 0.5. 

Another parameter which has very little effect on n(0*)n g is 7(H*)/7 (0*) . 
It is fortunate that 7/7. and 7(H*)/7(0*) have very little influence on mea- 
surements of n(0*)/n because other parameters of the ACF depend on n{0*)/n 

G G 

besides 7 /?( 0*), 7(0*), and 7(H*)/7(0*). Thus if the dependence of n(0*)/n 
can be eliminated from the other parameters of the ACF then better measurements 

of 7 /no*), 7(0*), and 7(H*)/7(0*) can be made. 

e 
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Usually the ratio n(0*)/n is found from TZERO which has been described 

previously. If n(0 + )/n is less than .55 n(0*)/n is found from the time delay 

at ACF » 1/2 (THALF). THALF is dependent on T fT{ 0 + ) , T(H + )/r(0 + ) # and T(0 + ) 

8 

besides being dependent on n(0*)/n , so the determination of »(0 + )/n is much 

more complicated than previously. If f(0 + ), T JT( 0 + ), or are 

not known then one must assume values for them. This leads to errors in the 

calculation of n(0 + )/n . These errors can be reduced if T(0*), TJT (0 + ), or 

r(H + )/r(0 + ) are calculated separately. For example, T( 0*) is calculated when 

n(0 + )/n is greater then .55. If !T(0 + ) is assumed to be constant with in- 

creasing altitude, then n(0*)/n g is found with better accuracy where also 

T/T( 0 + ) - 2 , (H + )/T(0 + ) - 1. 

& 

The determination of n is a bit more difficult than finding n(0 + )/n . 

This is because is found from measurements of a which depends on the elec- 
tron to ion temperature ratio. The parameter a also depends on the electron 
Debye length to wavelength ratio, 

a ■ 4 itD /A (5.10) 

where D is the electron Debye length given in equations 4.1 and X is the 
transmitted radio wavelength. The dependence of o on T / T . and o can be ex- 
pressed as [ Moorcroft , 1963] 

a " Ve /(1 + a2)Cl + T e /T i * a2) (5.11) 

where is Thomson scattering cross section given in equation (4.7). Usually 
in the F region a is very small so a only depends on n and T /T*. In the 
topside ionosphere the electron to ion temperature ratio is generally between 
one and three; therefore the electron density can be found within a factor 
of 2 since the measurements of o are quite accurate. 
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The measured power depends on many experimental factors besides depending 

on n and T IT .. If all these factors are known the electron number density 
a e % 

can easily be found, but usually ionosonde data are used to find the peak 
electron number density. (An ionosonde is a radio wave vertical sounder 
which sweeps in frequency from about 1 to 25 MHz.) This is then used to find 
the electron number density at the altitude where the peak power occurs. At 
other altitudes the electron number density can be found by normalizing the 
power by the peak power and multiplying this by the peak electron number den- 
sity found by ionosonde data. In Figure 5.4 are shown some electron number 
density profiles where o is assumed to be zero and T IT. is set equal to one. 
These profiles give an indication of variation of the electron number density 
during a winter day (Feb. 10, 1972). Notice that the peak electron number 
density increases before dawn and reaches a maximum in the late afternoon. 
After this the electron number density profiles decrease slowly during the 
night. The proportion of the 0 + and H + ions also varies during the day. This 
is shown in Figure 5.5 where the temperature of electrons, the 0 + and H + ions 
is equal to 1000 K. In Figure 5.5 there are two set of profiles - one for 
the daytime and one for the nighttime. The daytime profile indicates that 0 + 
is the major ion at all altitudes in the altitude region observed. The 
nighttime profile indicates that H + becomes the predominate ion around 750 km, 
that is, at 750 km the transition between 0 + and H + as the major ion takes 
place. Below 750 km at night 0 + ions are formed through a rapid charge ex- 
change with H + ions with the reaction 

H + 0 0 + H + AF 

The source of H + ions is the inward flow of H + ions from high altitudes. 

During the night H + ions diffuse inward which results in the maintainance of 
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A-THE EFFECT OF T(0 + ) NOT INCLUDED 
B-THE EFFECT OF T(0+) INCLUDED 
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Nighttime and. daytime profiles of 
n(0 + )/n $ for Feb. 10, 1972. There 
are profiles where the effect of the 
ion temperature is included (A) and 
profiles where the effect of the ion 
temperature is not included (B). 
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substantial ion densities in the F region. During the day inward diffusion 
of H* is not important due to the large photoionization rate of 0. In fact, 
during the day there is either a slight outward flow of 0* ions or diffusive 
equilibrium exists. This is the reason why 0* is the predominate ion from 
230 to 1200 km during the day. 

Also in Figure 5.5 the effect of the ion temperature is included in the 
calculation of nighttime and daytime profiles of n(0 + )/n . As mentioned previ 
ously, TZER0 depends slightly on T., So T. and n(0 + )/n are calculated from 
THALF and TZER0 respectively. THALF also depends on T/T. which can be calcu- 
lated from AMIN. The profiles of n(0 + )/n where the effect of T is included 
should be more accurate because the effect of all the physical quantities are 
included. The profiles of Figure 5.5 are almost the same except that the 
transition between H + and 0* occurs at a slightly higher altitude in the pro- 
file where the effect of r(0 + ) is included. The difference in the profiles 
of Figure 5.5 indicates that the errors in n(0 + )/n are less than 5 percent. 

The effect of T /T . on electron number densities is shown in Figure 5.6 
e i 

for daytime profiles of n . Profile 1 is for the case where T/T. is assumed 

0 0 u 

to be one. The other profiles include the effect of T/T. where T /T. is 
calculated from AMIN or from fitting the ACF. The difference between pro- 
file 1 and profile 2 or 3 is greatest around the peak of the electron density 
profile. This i; because the ratio T/T. is closer to two at altitudes near 
the peak. At other altitudes the effect of T/T. is reduced either because 

6 7 * 

n is much less than its value at the peak or T/T. is near 1.0. Also shown 

3 3 % 

in Figure 5.6 is a electron number density profile where T /T. has been cal- 
culated from nonlinear regression analysis. This method is described in 
Section 5.2.2 and is the method where the experimental ACF's are fitted to 
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theoretically generated ACFs. A library of theroetical ACFs is created for 

various values of T . and T /T * and the experimental ACFs are fitted using 

nonlinear least squares method. The values of T ^ and 3P /r^ corresponding 

to the best fit to the theoretical ACF are found from this method. The 

value of T IT. is then used in the calculation* of electron density profile 3 
0 x 

given in Figure 5.6. This electron density profile is very similar to pro- 
file 2 which indicates that the errors of electron density measurements are 
small. Although the profiles given Figures 5.6 agree fairly closely, above 
806 km there is a problem because the measured total scattered power is much 
smaller than it is around 400 km. Even though one expects the total scat- 
tered power to be less because of smaller electron densities, the reduced 
total scattered power causes larger uncertainities in the measurement of 

T /T.. Thus, electron densities are less accurate. This can be overcome by 

e x 

longer pulses which were used. 

The electron density profile for long pulse signals is shown in Figure 
5.7. There are three profiles shown which have been calculated using the 
three methods used for Figure 5.6. Profile 1 is the case T JT. ■ 1. In 

c? 

profile 2 T /T . is calculated from AMIN. By fitting experimental and theore- 
tical ACF's one gets T IT. used in profile 3. The largest difference between 
the profiles occurs below 666 km where T JT . is close to 2. Above 666 km the 
profiles are approximately equal. Profiles 2 and 3 are about equal from 342 
to 1206 km. This indicates the accuracy in values of TJT . is good because 
two different methods of finding TIT. give approximately the same results. 

The expected error of the electron densities given in Figure 5.6 and 
Figure 5.7 depends on the accuracy by which T JT. can be measured. The 
error of T IT. will be discussed later in Section 5.4. The error of the 
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Figure 5.7 Profiles of n for Feb. 10, 1972 f: 

9:07 to 9:17 Xm for long pulse dati 

where the effect of T /T . is shown 
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electron density also dopends on the error in the measurement of the total 
scattered power which depends on the number of correlated sables taken. As 
more samples are taken the error in total scattered power decreases. Usually 
enough samples are taken so that the error in the power is very small (less 
than 10’ 4 ). Therefore, the error of the electron density only depends on 
the error of T fT.. Typical values for the error of T/T. are less than 
25 percent. So the error of the electron density is less than 25 percent. 

5.4 Ion and Electron Temperatures 

While the determination of n(0*)/n and n are relatively quite simple, 

9 9 

the determination of the ion and electron temperatures is considerably more 
complicated especially if separate ion temperatures for 0* and H + are cal- 
culated. There are a number of parameters of the theoretical ACF which de- 
pend on T /T(Q*), T(H + )/T(0 + ) , and T(0 + ) . Therefore, one tries to choose 
parameters which depend on only one or just two of quantities T /T(0 *) , 
T(H + )/r(0 + ), and r(0 + ). 

Fortunately, /WIN depends on T /T(0*) and r(H + )/r(0 + ) and does not de- 
pend on r( 0 + ). AMIN also depends on «(0 + )/n which can be calculated from 
TZERO. When n(0*)/n is between 0.4 and 0.6 temperature ratios deduced from 

9 

AMIN are unreliable and another parameter of the ACF must be used to deduce 
temperature ratios. Usually AMIN is very useful in deducing temperature 
ratios for daytime because only at high altitudes does n(0*)/n drop below 0.6. 
The ion and electron temperatures discussed in this section apply mostly to 
daytime results. Daytime is concentrated on because the temperatures during 
the day are much larger than at night. Also during the day the ion and elec- 
tron temperatures vary more with altitude than at night. 

At nighttime H* becomes a dominant ion at much lower altitudes; so 
n(0*)/n & drops below 0.6 at lower altitudes than during the day. At night 
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the ion and electron temperatures seem to be fairly constant with altitude. 

In Figure 5.8 ion and electron temperatures are shown for nighttime on Feb. 

10, 1972 where a single ion temperature is assumed. This night was charac- 
terized by low magnetic activity. The temperatures were calculated from 
comparing theoretical and experimental values of AMIN and THAU. The electron 
and ion temperature profiles shown are almost identical at all altitudes. 

They can be considered the same since experimental errors are greater than 
the difference between the electron and ion temperatures . 

In Figure 5.9 ion and electron temperature profiles are shown for Feb. 

10, 1972 from 9:07 to 9:19. These profiles were calculated from AMIN and 
THALF using the method described in Section 5,2.1 where a single ion tempera- 
ture was used. Notice that the electron temperature profile is significantly 
greater than the ion temperature profile. The ion temperature seems to stay 
relatively constant up to about 500 km. Then the ion temperature increases 
steeply up to about 700 km. Above 700 km the ion temperature increases more 
slowly but there is a definite temperature gradient. The electron temperature 
increases rapidly from 230 to about 400 km. Above 400 km the gradient of the 
electron temperature gradually decreases until the gradient of the electron 
temperature is about 1 K/km at 1000 km. 

The ion and electron temperatures can also be calculated by fitting the 
experimental ACF to a theroetical ACF. A library of theoretical ACF’s used in 
this analysis is given in Appendix D along with the computer program which 
created the library. In this analysis the ratio n(0*)/n is first calculated 
by matching theoretical and experimental values of TZER0. Then the experi- 
mental ACF is fitted to a theoretical ACF by minimizing the square of the 
difference of the two ACF's. The results of this analysis are shown in Figure 
5.10, 5.11, 5.12, 5.13, and 5.14. 
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In Figure S.10 profiles of 7(0*) and T g are shown for Feb. 10, 1972 from 
9:07 to 9:29 AM. This figure is for two consecutive 12 -minute time periods. 

The circles represent data for the first time period. The data for the se- 
cond time period are given by X's. Data from both long and short pulses are 
given. The data above 900 km are not very accurate (±200 K) due to the lack of 
an extensive library of ACF’s for ion temperatures greater than 2000 K. The 
accuracy of the profiles below 900 km is approximately ±50 K. The accuracy 
of the profiles below 900 km is much better because of the extensive library 
of theoretical ACF's from 7( 0*) ■ 800 to 2000 K. Also, one can see that the 
accuracy below 900 km is good because of the close agreement of data of the 
two time periods. 

Comparison of the results of the methods of fitting ACF parameters (Fig- 
ure 5.9) and fitting the entire ACF (Figure 5.10) indicates a close agreement 
of thj profiles of 7(0*) and T g . Therefore, the results for 7(0*) and T g are 
very reliable. 

When n(0*)/n is approximately 1.0, 7(H*) is very difficult to deduce 
from experimental data. 7(H*) can be found from analysis of the profiles: of 
7(0*) and T g . Profiles of 7(0*) depend on 7(H*) because 7(0*) has a different 
shape for different circumstances. For example if a single ion temperature 
is used then 7(0*) becomes equal to the neutral temperature at a higher alti- 
tude than that when each ion has different temperatures. Inspection of Figure 
5.10 indicates that 7(0*) approaches a neutral temperature of 1000 K at about 
500 km. If 7(0*) were equal to 7(H*), 7(0*) would approach 1000 K at about 
600 km or higher. So one can conclude 7(H*) must be significantly larger than 
7(0*) in altitude region from 500 to 600 km because a larger 7(H*) increases 
7 (0*) . 


TEMPERATURE (°K) 


Figure 


.10 0 + ) and T & for Feb. 10, 1972 from 

9:07 to 9:29 AM where fitting the 
experimental ACF to theoretical ACF's 
is used to find T( 0 + ) and T g . Data 
for two consecutive time periods are 
shown. Circles are for 9:07 to 9:19 
x's are for 9:19 to 9:29. 
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In Figure 5.11 the ratio T(H*)/T(Q*) is shown for Feb. 10, 1972 from 
9:07 to 9:27 AM. The circles are for the first time period, X's are for the 
second. Both long pulse and short pulse data are shown in this figure. The 
lower altitudes don't have data points for the first time period because 
«(0 + )/n was nearly 1.0 at these altitudes during this time. During the next 
time period nC0 + )/n g was slightly below 1.0 which allowed the calculation of 
2 , (H + )/!T(0 + ) . Figure 5.11 shows that 7’(H + )/2’ (0 + ) has a maximum near 550 km of 
about 1.45 and drops off to 1.0 at low and high altitudes. 2 , (H + )/2 , (0 + ) is 
very large from 300 to 700 km. The fact that T(H + )/r(0 + ) has these large 
values indicates that much more heat is being transferred to H + from the elec- 
tion than to 0 + . The theory of thermal processes with regard to ion and elec- 
tron temperatures is discussed in Chapter 6. 

In Figure 5.12 the profile of n(0 + )/tt g for Feb. 24, 1972 is shown. This 
figure includes both long pulse and short pulse data from 7:57 to 8:24 AM 
where n(0 + )/n e is calculated from TZERO. The circles represent data for 7:57 
to 8:10. The X's represent data for 8:10 to 8:24. The morning of Feb. 24, 
1972 was characterized by disturbed magnetic activity. From 300 to 600 km 
n(0 + )/« g is nearly equal to 1.0. Long pulse data (300 km altitude sample) 
give values slightly less than 1.0 at 342, 450, and 558 km which allows the 
calculation of T(H + )/!F(0 + ) at these altitudes. 

The profiles of T(0 + ) and T for same data given in Figure 5.12 are 
shown in Figure 5.13. Above 900 km the accuracy of the profiles is approxi- 
mately ±300 K. This is due to the lack of an extensive library of theoreti- 
cal ACF's for T(0 + ) greater than 2000 K. Below 900 km the error in the pro- 
files is less (±100 K) than above 900 km but the error is larger than in Fig- 
ure 5.10. This is because the data in Figure 5.12 are earlier in the morning 
than Figure 5.10. At night the ion and electron temperature is approximately 
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Figure S.12 Profile of n( 0 )/rt^ for Feb. 24. 1972 
from 7:57 to 8:24 AM. Data represents 
two consecutive time periods. Circles 
are for 7:57 to 8:10. *'s are for 8:10 

to 8:24, 
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Figure 5, 


3 !T(0 + ) and T e for Feb. 24, 1972 from 

7:57 to 8:24 AM. Data represents 
two consecutive time periods. Circles 
and *'s are the same as in Figure 5.12 
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constant with altitude and much less than during the day (see Figure 5.8). 

So, at dawn the ion and electron temperatures must change as a result of 
solar radiation creating photoelectrons which in turn heat the ambient elec- 
trons. One can expect rapid changes immediately after dawn with more gradual 
or slower changes into the morning. Therefore, the heat transferred to the 
electron and ions could be changing more at 8:00 AM than at 9:00 AM. Even 
though the changes are slight this probably is the reason why the errors of 
Figure 5.13 are larger than those of Figure 5.10. 

Although the errors of Figure 5.13 are larger than those of Figure 5.10, 
Figure 5.13 still gives a good picture of the profiles of T( 0 + ) and T . In 

« y 

fact the profiles of Figure 5.13 are similar to the profiles of Figure 5.10. 
T(0 + ) in both Figures 5.10 and 5.13 is approximately 1000 K up to 500 km. 

From 500 km to 700 km T(0 + ) increases rapidly to 2000 K. Above 700 km the 
gradient of T(0 + ) decreases gradually until it is about 1 K/km at 1200 km. 

T in both Figures 5.10 and 5.13 starts at 1000 K at 230 km and increases 
steadily to 2500 K at 600 km. At 900 km T is approximately equal to 3000 K 
in both Figures 5.10 and 5.13. 

In Figure 5.14 the ratio r(H + )/T(0 + ) is shown for the same data used for 
Figures 5.12 and 5.13. The maximum of T(H + )/2’(0 + ) occurs at about 550 km and 
is approximately 1.4. There are substantially large values cf T(H + )/!r(0 + ) 
from 450 to 770 km. As in Figure 5.11, ?(H + )/T(0 + ) again falls off to 1.0 at 
low and high altitudes. 

5. 5 Smmapy 

The results given in this chapter apply to data taken at the Arecibo 
Observatory in Puerto Rico. There are two methods used to reduce the data: 

1) fitting the parameters of the ACF and 2) fitting the entire ACF. These 
methods are described in Section 5.2. There are five physical quantities 





Figure 5.14 F(H + )/r(0 + ) for Feb. 24, 1972 from 

7:57 to 8:24 AM. Data represents 
two consecutive time periods. 
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which can be reduced from this analysis. They are n , n(0 + ), T(0*) t r(H + ), 

& 

and T where n is assumed to equal n( 0 + ) plus n(H + ). The quantities n a and 
T(0*) can be found directly from this analysis. The other quantities are 
found by first computing n(0*)/n , T(H*)/f(0*), and T JT{ 0 + ). 

In Section 5.3 the results of n and n(0 + )/n are discussed and it was 

0 & 

determined that n and n(0 + )/n can be found very accurately. Figures 5.2 

& 6 

and 5.3 show that n(0 + )/« could be computed to within +.05 even if T(0+) is 
not known. If T(0 + ) is known then n(0 + )/n could be calculated with errors 
les.s than .01. The determination n g is less accurate than the determination 

of n(0 + )/n because n depends on an accurate measurement of T /T(0 + ). In 

& & & 

the topside ionosphere T/T(0*) is generally between one and three, so n 
can be found within a factor of 2. Usually T /T( 0 + ) is calculated to within 
±,1 or less so n g can be computed with a relative error less than 4 percent. 

In Section 5.4 results of ion and electron temperatures are discussed. 

It was determined the profiles of r(0 + ) and T were reliable and fairly accu- 
rate because of the general agreement of the profiles for both methods used 
(see Figures 5.9 and 5.10). Using a library of theoretical ACF's and fitting 
the experimental ACF to the best theoretical ACF it was determined that !F(0 + ) 
and T g could be calculated to within ±100 K (Figures 5.10 and 5.13). The fact 
that the results of one 12 -minute time period are almost exactly the same as 
the results of the next 12-minute time period (Figures 5.10 and 5.13) indi- 
cates that the accuracy of these results is very good. The general agreement 
between Figures 5.10 and 5.13 and between Figures 5.11 and 5.14 suggests that 
these results are typical of morning data because two different days are re- 
presented. These days include both quiet (Feb. 10, 1972) and disturbed (Feb. 
24, 1972) magnetic activity. 
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Figures 5.11 and 5.14 show the ratio 3’(H + )/T(0 + ) can have significantly 
large values in the altitude region between 350 and 700 km. This means that 
T(H + ) can be significantly larger than T(0 + ) in this altitude region with 
the maximum difference occurring near 550 km. Although the accuracy of the 
values of T(H + )/J’(0 + ) given in Figures 5.11 aid 5.14 isn't too good it is 
certain that !T(H + ) is larger that T(0 + ) in the altitude region between 350 
and 700 km. The larger values of T(H + ) has an effect on the r(0 + ) profile 
because !T(0 + ) approaches the neutral temperature at 500 km rather than at 
600 km. Analysis of the ion and electron temperatures based on the energy 
balance equations is given in Chapter 6. 
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6. ION AND ELECTRON TEMPERATURES FROM THE ENERGY BALANCE EQUATIONS 
6.1 Introduction 

During the past IS years some progress has been made in describing the 
ion temperature and electron temperature structure in different regions of 
the magnetosphere. By equating the electron-ion heating rate with the loss 
of energy of the ions to neutrals, Hanson [1963] was the first to show that 
the oxygen ion temperature profile depends heavily on the altitude. His re- 
sults show that the ion temperature increases gradually from the neutral gas 
temperature at 300 km towards the electron temperature above 800 km. Dalgamo 
et al. [1963] investigated the effect of solar ultraviolet radiation on the 
rate of heating of the ambient electrons and found that the electron tempera- 
ture reaches a maximum temperature near 200 km at noon in a quiet ionosphere. 
Geisler and Bowhill [1965] developed a theory of nonlocal heating of the 
topside ionosphere and the protonosphere by photoelectrons which escape the 
F 2 layer. They found that the electron temperature has a stjep temperature 
gradient at altitudes less than 1000 km for sunspot maximum conditions. The 
effect of oxygen cooling on ionospheric electron temperature has been re- 
ported by Dalgamo et al. [1968] who found that the inclusion of electron 
cooling by excitation of the fine-structure levels of atomic oxygen removes 
the discrepancy between Thomson scatter data and theoretically predicted 
electron temperatures above 320 km. The influence of the ion thermal conduc- 
tivity has been studied by Banks [1967b] and Sanatani and Hanson [1970]. 

Their calculations of the ion temperature resulted in an ion temperature sig- 
nificantly less than the electron temperature at high altitudes. In the 
studies by Banks [1967b] and Sanatani and Hanson [1970] the electron tempera- 
ture was assumed to be constant above a certain altitude. 
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Theoretical profiles of temperature for different model neutral atmos- 
pheres and different ion and electron density models are computed from the 
basic steady state energy balance equations for ions and electrons described 
in chapter three. In the following sections results of these calculations 
are discussed. Also, in this chapter the effect of photoelectron heating on 
electron temperatures will be examined. The importance of thermal conduction 
in both the electron and ion energy balance equations will be studied too. 
Other topics studied in this chapter are the -"feet. of temperature gradients 
at the upper boundary and the effect of various elastic and nonelastic energy 
transfer processes on electron and ion temperatures. 

6 2 numerical Methods 

The electron, 0 + , and H + energy balance equations are coupled differen- 
tial equations which are second order in a (height). These equations cannot 
be solved analytically, so numerical methods must be used to find the solu- 
tion. In order to find the suitable method among the possible methods avail- 
able, one employs a trial and error procedure. A method which works on one 
equation under certain conditions may not work on the other equations. For 
example, a triple diagonal mesh method (which will be explained later) is 
effective in solving the ion energy balance equation, but not the electron 
energy balance equation. Thus, one must search for the best method to numeri- 
cally solve each energy balance equation. 

6.2.1 The triple diagonal mesh method. The numerical method known as 
the triple diagonal mesh method is used to solve the ion energy balance equa- 
tion when thermal conduction is included as an energy transfer term. It was 
used by Banks [1967bl to solve for the ion temperature profile under restric- 
tive conditions. In the following analysis, a modification of the basic 
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procedure of Banka [1967b] is used, 

Banka [1967b] assumed a single ion temperature and an electron tempera- 
ture that was constant with altitude. Also, he assumed a zero temperature 
gradient as the upper boundary condition. These assumptions are modified in 
this study so that the 0* and H* ion temperature profiles can be found. An 
electron temperature profile from the electron energy balance equation is 
used. Also, the upper boundary condition is modified so that a variable 
temperature gradient can be used. 

The steady- state ion energy balance equation is 

ZT . 

-sin 2 I |- -gji - Pfe-i) - L(i-n) - C(i-j) (6.1) 

where I is the magnetic dip angle, K. is the ion thermal conductivity, P(e-i) 
is the energy transfer between electrons and ions, Df.'-n) is energy transfer 
between ions and neutrals, C(i-i) is the energy transfer between one ion and 
other ion species, i is a subscript for the ion species (G* or H + ), and J is 
a subscript denoting the other ion species. 

The ion thermal conductivity K. for a single ion gas is 

K. • 4.6 x 10 4 (7’. 5 /A .)**/ 

*t % T 

[1 ♦ 3.22 x loV, 2 /*.) Z n ? n ] eV cm" 1 sec" 1 K _1 (6.2) 

1 If U 

Substituting 0 » T. 7 ^ 2 in the energy balance equation gives 
a 2 * 

■2-4 » /.(a,0) (6.3) 

32 

where ^(*,0) « -7 (P(e-f) - L(i-n) - Ctf-fty {2F.sin 2 l) 
and F. - 4.6 x lO 4 /^.’ 4 [1 ♦ 3.22 x 10 4 (J\ 2 /n.) I n Q n ] 

t Z i. if 71 U 


. J 
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If one assumes that 

0(a) ■ ^(a) + A(») (6.4) 

where 0 (a) is known from a previous iteration or initial conditions and 
A (a) is a difference quantity, equation (6.3) can be put in a numerically 
solvable form after linearization of /j(a, 0). Linearization of /j(*,0) gives 


/*! («,e) ■ / 2 (a, 0 ) ♦ / 3 (a, 0 ) 


(6.5) 


where »nd are function which depend on the energy production 

2 

and loss rates and factors of the order of A and higher are neglected. Using 
and (6.5) with (6.3) one gets 


2 2 

—4 ■ / 2 (*,0) ♦ A(** 0 )A - 

3a z 2 3 3a^ 

Applying central difference formulas to (6.5) results in 


Vi * Vf * Vi ' 


( 6 . 6 ) 


(6.7) 


where B. ■ -2 - 6 

D j " ^2 (a ^ ^ fy*l “ * fy-l' 

6 is the altitude step size (5 km), and the *s refer to integration steps 
extending from j ■ 0 to a maximum value, j'max. 

Since (6.3) is a icond order differential equation one needs two bound- 
ary conditions. The first boundary condition is that the ion temperature is 
equal to the neutral temperature at an altitude of 200 km. This gives 
A q * 0. The second boundary condition is that the ion temperature profile 


has a temperature gradient at an altitude of 1200 km of 37\/3a|p 00 . Typical 
experimental values of 37\/3a|p 0 g lie between 0 and 2 K/km. This gives 

5/7 


/max ♦ 1 


A. ♦ 7/2 
jinas, 


g .' 3f./3a .~ nn 

y jmax t 1200 


( 6 . 8 ) 
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Since (6.7) is in the triple diagonal form, the solutions for the A^ 
are found by first calculating mesh coefficients Kj and L ^ for j ■ 0 to 
Jmax according to 


»i • 4 0 

~~ sr - 


* 

j B. - L 


. 1 
r l " 


(6.9) 




7 j-i 

Then one works downward from j ■ Jmax to j ■ 1 using the relation 


Vi " *M 




( 6 . 10 ) 


with 


. ^max “ fymax+1 * ^/max + ^i^ s jl200 fymax^ 

Vx.l (1 * 


( 6 . 11 ) 


The iterative solutions for 8 (and T are found by using equations 
(6.7) - (6.11) to find A^ which is added to the previous value of g The 
initial function g is found from the temperature profile which neglects the 
thermal conduction effects. The iteration process is continued until all the 
A •' s are essentially zero. The final function g is the equivalent to 8 which 

4? 

gives the steady- state ion temperature profile. The convergence of the iter- 

-4 


y*fl 


for all the j*s reaching 10 


ation process is rapid with the ratio |A. 
after four repetitions of the process. 

The profiles for each ion temperature are found by solving for 7(H*) 
using the triple diagonal mesh method where 7(0*) is set equal to the ion 
temperature where a single ion temperature is assumed. Then, the triple 
diagor.cl mesh metnod is used again with 7(H*) given from the previous step 
to solve for r .'(0*). The procedure is continued until there is no change in 
either the 7(0*) or the T(H*) profile. It generally takes four repetitions 
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of this iteration process to get convergence of both 0 + and H* temperature 
profiles. 

6.2.2 Modified Newtoti 'e method. A modified Newton's method is used to 
solve for the ion temperature profiles when thermal conduction effects are 
ignored. Thus, the ion temperature is calculated using a formula of the form 

x Ul = x i ~ afix 

where x is the ion temperature, a is a modification factor, f(x) is the right 
side of equation (6.1), and fix) is the derivative of f(x) with respect te x. 

To find the 0 + and H + temperature profiles, it is first necessary to find 
the single-ion temperature profile. Next, the H + temperature profile is 
found, followed by the 0 + temperature profile. Finally, the process of 
finding the H + temperature profile and then the 0 + temperature profile is 
repeated until the respective profiles converge. 

6.2.3 Integration of the energy balance equation. The electron energy 
balance equation is integrated numerically to solve for the electron tempera- 
ture profile. The integrated steady state energy equation is 


dl f » 

sin I K e E J (.Q - L) dz 


( 6 . 12 ) 


where Q is the heat input through photoelcctron heating of electrons, L is the 

heat lost by Coulomb, elastic, and nonelastic collisions with ions and neutral 

particles, and K is the thermal conductivity of the electron gas given in 
* ^ 

equation (3.86). 

The upper limit of the integral in equation (6.12) is infinity. This is 
the result of using the assumption of open field lines or flat earth. In 
other words, the field lines are parallel out to infinity. This is a very 
good approximation in the altitude region below 1200 km. The integral is 
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not actually integrated to infinity. The upper boundary condition gives the 
value of the integral from 1200 km to infinity. This is done by setting the 
value of the integral from 1200 km to infinity equal to the left side of 
equation (6.12) where the upper boundary condition is on dT g /da. 

The lower boundary on T is that thermal conduction is unimportant below 1 

6 | 
200 km. That is, thermal conduction is negligible compared to heating by } 

photoelectrons or cooling by neutral particles. 

The electron temperature profile is found by adopting a trial value of 

T at 1200 and integrating equation (6.12) downward in a stepwise fashion 

(in 5 km steps). The temperature gradient, dT /da, is determined from equa- 

tion (6.12). The energy equation is used to calculate the second derivitive 

of T . /dfa 2 , where the value of dT /da from equation (6.12) is used. The 

0 e e 

electron temperature at the next lowest altitude is found by using a Taylor 
expansion 


<£Z\ 


V» ' Aa) ■ V 3) - la 


As 


d?T 

1 a ._2 


Lb 


(6.13) 


l s * d*r 

where Aa ■ 5 km. An appropriate value of T at the upper boundary will give 
a profile which satisfies the lower boundary conditions (see Section 6.4). 

The procedure for finding the electron temperature profile is summarized in 
Figure 6.1. 

6.2.4 Time relaxation method . The time relaxation method is sometimes 
used for investigation of diurnal changes in ion and electron temperatures. 

It could also he used to calculate steady- state ion and electron temperatures. 

The time relaxation method as described by Da Rosa (.1966) involves solv- 
ing the energy balance equation which can be written in the form for ions as 
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1. Input T g at 1200 km 

2. Integrate the heat 
source, loss terms 


3. Find dTjdz 

e 

4. Find d 2 T f da 2 

e 

5. T (a - As) = 

/(a) - Aa - h ^ T e Aa 2 
e da da 2 


6. Continue 1-5 until lower 
boundary condition is 
satisfied 


Figure 6.1 Schematic diagram of the procedure for finding the 
electron temperature profile. 
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2 2 

II - i-J ♦ t^a.t.mH) + (6.14) 

where i/(a,t,IP) * (y n e^ -1 sin 2 * ^ 

- (y^fc)* 1 | sin 2 


V(a,t,r) - Cjn^r 1 (P(e - t) - L(€ - »)) 

Central difference formulas can be used to reduce equation (6*14) to 

r* #, , * r> * f * . * * * 2f • * + r* 


. £/. . > - W,~ r .*4 

At t,j (A«r 


+ V., . 
i**0 


T. , , - T. . 

Jlhi LuL 

A* 


♦ V. . 


(6.15) 


where i is the subscript for the height step. At; 
j is the subscript for the time step, At; 

and U. . * I/(tAs, j'At, T.), etc. 

*• *0 * 

The unknown T. . , can easily be found since all the other quantities in 
1*0 + 1 

equation (6.15) are known. The height and time steps have to be chosen care- 
fully so that the solution of the difference equation converges to the solu- 
tion of the differential equation. Sometimes even if Aa and At approach zero 
there is no convergence (the values of T oscillate wildly with each iteration) 
The equation is then said to be unstable. The condition of stability as given 
by Da Roea [1966] is 


max (V. .) -- -* 

%t0 (Aar 


1 

I 


(6.16) 


This condition imposes strict restrictions on the steps Aa and At, and makes 

it impossible to solve the physical problem, as verified by the following 

4 -3 

example: if the electron number density is 10 cm and the temperature is 

1000 K, then with a height step of 20 km, the maximum time step is 0.1 seconds 


181 


Even time intervals of one hour would take far too many iterations and also 
would render the solution meaningless by the accumulation of truncation and 
round-off errors. 

Different schemes have been formulated to circumvent these difficulties, 
but care must still be taken in choosing the time step. Methods other than 
the time relaxation method were found to be easier to apply to this study, 
and had fewer difficulties. 

If one is interested in the thermal response to changing electron num- 
ber density or varying photoelectron-thermal electron heating one would prob- 
ably take the extra care in employing the time relaxation method. Fortunately, 
this study involves the solution of the steady-state ion and electron energy 
balance equations which can more easily solved by methods described in Sec- 
tions 6.2.1, 6.2.2, and 6.2.3. 

6.3 The Effect of the Electron and Ion Number Densities on the Electron 

and Ion Temperature Profiles 

In this section several different models of the electron and ion number 
density profiles corresponding to different atmospheric conditions will be 
used to calculate electron and ion temperature profiles. In the two cases 
considered in this section, parameters other than number densities will be 
held constant. In the neutral atmospheric model used, the exospheric tempera- 
ture is equal to 1000 K (CIEA [1972]; see Appendix A). Also a thermal -electron 
heating rate at a solar zenith angle, x* of 0° from Mantas [1973] will be 
used (see Figure 3.3). The heating rate above 1000 km is assumed to have a 
constant value set equal to the heating rate at 1000 km. 

The electron temperature profile is calculated by integrating the elec- 
tron energy balance equation which includes thermal conduction. The boundary 
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condition at 1200 km Is that the gradient of the electron temperature is 
1 K/km. The lower boundary condition is that thermal conduction is unimpor- 
tant at the altitude of 120 kn. Below 120 km, the electron temperature pro- 
file is controlled by energy loss to the neutral gases, so that the electron 
temperatures are nearly equal to neutral temperatures . 

The ion temperature profile is calculated using the triple diagonal 
mesh method. Thermal conduction is included in the ion energy balance equa- 
tion. As in the case of the electron temperature profile, the upper boundary 
is applied at 1200 km and is that the gradient of the temperature is 1 K/km. 
The lower boundary condition is determined by the fact that thermal conduction 
is unimportant at altitudes below 200 km. The ion temperature is controlled 
by energy loss to the neutral gases at altitudes below 200 km. Thus, the ion 
temperature becomes eqv.al to the neutral temperature near 200 km. 

It is necessary to consider two cases depending on the electron and ion 
number densities. In the first case the effect of the magnitude of the elec- 
tron number density on the temperature profiles is considered. Electron and 
ion temperatures (Figure 6.3) have been calculated with the electron number 
density equal to one, one half, and two times that shown in Figure 6.2. The 
0 + plus H* number density shown in Figure 6.2 equal the electron number den- 
sity of Figure 6.2. When the electron density is twice and one half that 
shown in Figure 6.2 m( 0 + ) is found by calculating n(0*)/n from Figure 6.2 
and multiplying this by twice or one-half r»(H + ) is found by subtracting 
n(0 + ) from n , 

When the electron number density is smallest, the electron temperature 
is largest. Also, the separation between the electron and ion temperatures 
increases as the electron number density decreases. This is the result of 
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Figure 6.3 T e and T i profiles for x ■ 0* and different electron 
density profiles. 1 refers to the profile in Figure 
6.2. 2 refers to 2 tines that profile. % refers to 

0.5 that profile. 
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two factors. First /the heating of the ions depends on both n g and 
as n g decreases and T g increases the heating of the ions remain* relatively 
unchanged. Second, the effect of the electron density on the ion temperature 
profile is not as great as its effect cm the electron temperature profile. 
While the electron density determines the magnitude by which electron energy 
is lost to other particles, the electron density is important in both ion 
energy production and loss. 

In the second case the effect of ion number density profiles on electron 
and ion temperatures is considered (Figure 6.5). In Figure 6.4, the ion 
number density profiles used to calculate the temperature profiles of Figure 
6.5 are shown. Two examples are used to show the effect on electron and ion 
tenperatures when 0* dominates at all altitudes and when H + becomes the pre- 
dominant ion at a low altitude. These two examples illustrate the effect of 
Coulomb collisions on electron and ion temperatures. 

Since H* has a larger energy transfer rate than 0*, the ion temperatures 
are higher when H* is the dominant ion than when 0* is the dominant ion. 

This is indeed the case, as is shown in Figure 6.5. Also, since the electron 
number density is essentially the same at altitudes near 1200 km, the electron 
temparature profiles A and B are almost identical. Therefore, the species 
of ion present has very little effect on the electron temperature profile 
at 1200 km. 

Most of the effect of changing the ion density profiles is felt in the 
ion temperature profiles. Since the ion density profiles change from night 
to day there will be changes in the ion temperature profiles from this factor 
alone. 

The 0* ion density is dominant at all altitudes between 200 and 1200 km 


ALTITUDE (km) 




ALTITUDE (km) 



TEMPERATURE (*K) 

Figure 6.5 Ion and electron temperature profiles for different ion 

number density profiles. The A and B temperature profiles 
were calculated using the respective A and B ion density 
profiles of Figure 6.4. 
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during the day. This corresponds to profile A to Figure 6,4. At night there 
is no photoionization source for 0* ions, but toward flow of H* ions and ac- 
cidentally resonant charge exchange acts as a source of ionization for 0* 
ions. Thus, there are substantial 0* ion densities at night, although H 4 
becomes important at lower altitudes. This corresponds to profile B to Fig- 
ure 6.4. The nighttime profiles of the electron and ion temperatures are 
discussed in the next section. 

6.4 Changes for the Ion and Electron Temperature Profile from Eight to Dag 

The electron and ion energy balance equations can be solved for the case 
where no solar radiation is -received. This implies that photoelectrons are 
no longer being produced. Before dawn and after dusk the ambient electrons 
can be heated by photoelectrons which originate in the conjugate ionosphere. 
This occurs if the conjugate ionosphere is still sunlit. During the night 
there are no photoelectrons, so by neglecting the heating of ambient electrons 
by photoelectrons, nighttime ion and electron temperature profiles can be 
calculated. Throughout this section the neutral atmosphere model is CIRA 
1972 with T m ■ 1000 K. 

In Figure 6.6, the effect of zero heating of the ambient electrons on 
the electron temperature profile is shown. This case also uses the ion 
density profiles shown in Figure 6.4 (profile A). The electron temperature 
profile has been calculated using thermal conduction in the energy balance 
equation with upper and lower boundary conditions as discussed in Section 
6.3 using the method described in Section 6.2.3. A fixed ion temperature 
is used. The most significant effect of zexv Seating of ambient electrons is 
that the electron temperature profile becomes equal to the ion temperature 
profile at 400 km and below as shown in Figure 6.6. There is no maximum 
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Figure 6. 


6 Electron temperature profiles for night and day where the 
upper boundary condition is ■ 1 K/kra for both night 

and day. The ion number densities used to calculate T 
are shown in Figure 6.4 A. 
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of the electron temperature at 2S0 tat, as there is in the daytime electron 
temperature profile. At higher altitudes the electron temperature is less 
at night than during the day f but it is still greater than the ion tempera- 
ture. This is the result of heat flux by thermal conduction from the higher 
altitudes. The biggest differences between the electron and ion temperatures 
occur between 500 and 700 km. This is because while ion-neutral energy trans- 
fer controls the ion temperature profile at higher altitudes, electron -neutral 
energy transfer becomes important at a lower altitude than the ion -neutral 
energy transfer. 

The electron temperature profile for the case where the ion densities 
are given by profile B in Figure 6.4 is shown in Figure 6.7. The lower and 
upper boundary conditions are the same as the previous example. In this case 
the ion temperature profile also is calculated using the method described in 
Section 6.2.1. As in the previous case the electron temperature is equal to 
the ion temperature below 400 km. The difference between the ion and electron 
temperatures is largest near 550 km. The ion temperature is still less than 
the electron temperature up to 1200 km. The electron temperature is equal to 
the ion temperature only in the region below 400 km. 

During the night it is possible that the temperature gradients at the 
upper boundary may be reduced from those during day. In Figure 6.8 the upper 
boundary conditions for the ion and electron temperatures are that the tem- 
perature gradients are zero and 0.5 K/km respectively. The lower boundary 
conditions are the same as the previous example and the ion densities are 
given by profile B in Figure 6.4. The ion and electron temperatures in Fig- 
ure 6.8 are reduced more than 1000 K at altitudes near 1200 km than those in 
Figure 6.7. Also, in Figure 6.8 the electron temperature profile is approxi- 
mately equal to the ion temperature profile at all altitudes. The electron 
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Figure 6. 


8 Electron and ion temperature profiles 
for nighttime where the upper boundary 
conditions are 3!T /3a ■ ‘0.5 K/km and 
3T e /3a ■ 0. The ion density profiles 
are shown in Figure 6.4 B. 
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and Ion temperature profiles increase slowly from 1000 K at 400 km to 1500 K 
at 1200 km. Therefore, if the temperature gradients at the upper boundary 
are zero, the electron and ion temperatures increase only 500 K above the 
neutral temperature (Figure 6.8); whereas, if the temperature gradients at 
tho upper boundary are 1 K/km, the electron and ion temperatures increase 
more than 1500 K above the neutral temperature (Figure 6.6 and 6.7). 

6.5 The Effect of Inelaetio Collieione on the Electron Temperature Profile 

Inelastic collisions of electrons with neutral particles are important 
as sinks of heat for electrons at altitudes below 400 km. Inelastic colli- 
sions were discussed in Section 3.7. Electrons can be cooled through rota- 
tional and vibrational excitation of N 2 and 0 2 . Another important sink of 
heat for electron is the fine structure excitation of 0. These inelastic 
collisions control the electron temperature profile at the lower altitudes 
in the F region. They gradually become unimportant with increasing altitude 
due to decreasing neutral densities. 

The importance of the various electron heat transfer processes is shown 
in Figure 6.9. These energy transfer rates were calculated using the CIRA 
[1972] neutral atmosphere model with T ^ ■ 1000 K, the ion density profiles 
given in Figure 6.4 profile A, and the electron and ion temperature profiles 
given in Figure 6.10 profile B. Curve 1 in Figure 6.9 is the cooling due to 
electron-ion Coulomb collisions. Curve 2 is the cooling due to electron- 
neutral elastic collisions. Cooling through vibrational excitation of N 2 
and 0 2 is shown in curve 3. Cooling through rotational excitation of N 2 and 
0 2 is shown in curve 4. Curve 5 represents the cooling through fine structure 
excitation of 0. Curve 6 is the heating of thermal electrons by photoelec- 
trons. In curve 7 energy transfer through thermal conduction is shown. 
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Figure 6.9 Electron energy transfer rates for all the 

heat transfer processes: 1 - cooling through 

electron-ion Coulomb collisions, 2 - cooling 
through electron-neutral elastic collisions, 

3 - cooling through vibrational excitation of 
N 2 and O2, 4 - cooling through rotational 
excitation of N 2 and 0 2 , 5 - cooling through 
fine structure excitation of 0, 6 - heating of 
thermal electron by photoelectrons, 7 - energy 
transfer by thermal conduction. 
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Figure 6.9 illustrates that thermal conduction is the most important 
energy transfer process at altitudes above 400 km. Photoelectron heating of 
the ambient electrons and cooling through Coulomb collisions of electrons and 
ions are the only other significant energy transfer processes above 400 km. 
Below 400 km there are many important energy transfer processes. Photoelec- 
tron heating of the ambient electrons is important below 400 km while cooling 
through Coulomb collisions loses its importance at altitudes i/Clow 250 km. 
Figure 6.9 shows that elastic collisions are not important at any altitude. 
Inelastic collisions are important below 250. Cooling through fine-structure 
excitation of 0 and through vibrational excitation of and O 2 are the im- 
portant inelastic cooling mechanisms. 

The importance of inelastic collisions in controlling the electron 
temperature profile at altitudes below 400 km is illustrated in Figure 6.10. 
Ion density profiles found using profile A of Figure 6.4, the CIRA [1972] 
neutral atmosphere model with 2*^ ■ 1000 K, and the thermal -electron heating 
rate from Figure 3.3 with x * 0° were used to calculate the electron tempera- 
ture profiles shown in Figure 6.10. Inelastic collisions act to reduce the 
peak in the electron temperature profile near 250 km, as can be seen by the 
fact that profile B is less than profile A in Figure 6.10. Profile A is the 
electron temperature profile without including the effect of inelastic colli- 
sions in the electron energy equation, whereas profile B does include the 
effect of inelastic collisions. Thermal conduction is also important in re- 
ducing the peak of the electron temperature profile in both profiles A and B. 
As can be seen from Figure 6.9 photoelectron heating is the important source 
of heat causing the peak in the electron temperature. Electron-ion Coulomb 
collisions channel much of this heat away from the electrons at altitudes 
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Figure 6.10 Electron temperature profiler for cases where 
inelastic collisions are not included in the 
energy balance equation (profile A) and 
inelastic collisions are included (profile B) . 
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above 300 km. However, below 300 kn thermal conduction and inelastic colli* 
sions must transfer this excess heat to other altitudes or other particles. 
Since inelastic collisions are not sufficient to absorb all the heat, there 
will be temperature gradients. Thus, a small peak in the electron tempera- 
ture profile occurs. 

6.6 Different Ion Temperatures 

The importance of considering different 0* and H + ion temperatures has 
been examined by Banks [1967a] for the case where thermal conduction is not 
included in the energy balance equation. He found that H + ions can have tem- 
peratures 250 K greater than 0* ions. In this section the separate ion tem- 
peratures are calculated with thermal conduction included in the energy bal- 
ance equations using the method described in Section 6.2.1. Several different 
cases are considered. Two sets of ion density are considered along with 
changes in the upper boundary condition. 

The lower boundary condition for H* and 0 + temperatures is ■ T n at 
200 km. The electron temperature's lower boundary condition is that thermal 
conduction is not important at 120 km. The thermal -electron heating rates 
are from Figure 3.3 with x * 0°. 

In Figures 6.11-6.13 the effect on the ion temperature profiles of chang- 
ing the upper boundary is shown. The ion number density profiles used are 
shown in Figure 6.2. The electron temperature profile is held constant in 

these calculations. In Figure 6.11, 37./ 3a ■ 0 where j is the subscript for 

J 

both H* and 0% The H* temperature is greater than the 0* temperature because 
H* has a larger electron-ion heating rate than 0*. The difference between 
the 0 + and H + temperatures is largest near 500 km. This is because ion- 
neutral cooling is larger for 0 + than H + . The 0* temperature is controlled 
more by ion-neutral cooling at the lower altitudes, so that 0* temperature 
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Figure 6.11 Temperature profiles for 0 + and H+ where the 

upper boundary condition is 37y/3s = 0. The ion 
number densities used to calculate 7(0+) and 7(H+) 
are shown in Figure 6.2, 
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Figure 6.12 Temperature profiles for 0* and H + where the 
upper boundary condition is 3^/33 * 0.5 K/km 
and the ion number densities are from Figure 6.2 
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Figure 6.13 Temperature profiles for 0* and H + where the 
upper boundary condition is 3 r t */3s "1.0 K/km 
and the ion number densities are from Figure 6.2. 
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is close to the neutral teaperature up to about 500 fas. 

The ion temperatures are such lets than the electron teaperature at the 
higher altitudes in Figure 6.11. There is a gradient in the electron teapera- 
ture and no gradient in the ion teaperature . A large aaount of heat is being 
transferred in the electro) gas through conduction, but very little in the 
ion gas. In fact, the ion teaperatures only increase by about 700 K above 
the neutral teaperature, while the electron teaperature is more than three 
times the neutral teaperature at 1200 kn. 

In Figure 6.12 the ion teaperature profiles is shown for the case where 
the upper boundary condition is 32^ 3* « 0.S K/ka. In this case the dif- 
ference between the 0* and H* teaperature profiles is larger than where 
3 Tj/te ■ 0. The largest difference occurs at about 5(h) ka. The ion teapera- 
ture profiles also increase to larger values in Figure 6.12 than in Figure 
6.12 than in Figure 6.11. This is because of h«t being transferred to lower 
altitudes through thermal conduction and thereby raising the teaperatures. 
Thus, as the teaperature gradients increase more heat is transferred from the 
high altitudes, and the temperatures increase significantly. 

The case where the upper boundary condition is 3T^/3a ■ 1 K/ka is shown 
in Figure 6.13. In this figure, the difference between the 0* and H* tempera- 
ture profiles is even larger than the previous two cases. As before, the 
largest difference occurs near 500 ka. There are significant differences 
even up to 1200 ka. Both the 0* and H* temperatures increase to larger values 
at 1200 ka than they do when 3J\/3* ■ 0 and 0.5 K/ka. This is due to more 
heat being transferred from above the upper boundary to lower altitudes. 

More heat is transferred to lower altitudes by thermal conduction if teapera- 
ture gradients are larger because thermal conduction depends on temperature 
gradients. 
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Figures 6.11, 6.12, end 6. IS shoe the effect of different upper boundary 
conditions on profiles of 7(0*) and 7(H 4 ). These profiles were calculated 
based on an electron teaperature profile which was the sane in each figure. 
This is not a self-consistent calculation because the electron teaperature 
profile is »cted by changes to 7(0 + ) and 7(H + ). Figures 6.11, 6.12, and 
6.13 still provide useful information on the effect of the upper boundary 
on ion teaperature profiles. 

Using Figures 6.11-6.13, H*-to-0* teaperature ratios, 7(H + )/7(0 + ), can 
be calculated. These ratios are shown to Figure 6.14. This figure shows 
that ?(H*)/?(0*) is largest in the altitide region nsar 500 ka. Curve A 
represents T{H*)/T(0*) for Figure 6.11 where 32V 3 a ■ 0. Curves B and C are 
for Figures 6.12 and 6.13, where 3 7^/3a ■ 0.5 and 1.0 K/ka respectively. 

Curve C has the largest values of 7(H 4 )/7(0 + ). The teaperature ratio 
T(H 4 )/T(0 4 ) is large in curve C even at altitudes near 1200 ka. 7(H 4 )/7(0 4 ’) 
is much smaller at 1200 km in curves B and C. Since curve C represents the 
case 37^/3* ■ 1.0 K/ka, this means that more heat can be conducted downward 
through the H* ion gas. Thus the H* ion temperature is larger than the O* 
ion temperature. The peak in ?(H*)/?(0 4 ) is due to the relatively constant 
teaperature gradient of ?(H 4 ) with altitude. The 0* teaperature decreases 
more rapidly and approaches the neutral teaperature at about i»00 ka. 

The ion density was found to he aportant in controlling the ion and 
electron teaperature in Section 6.3. The ion density is also important in 
controlling the difference between the H* and 0* ion temperatures . The ef- 
fect of the ion density profiles A and B of Figure 6.4 on 0* and H* tevpera- 
ture profiles is shown in Figures 6.15 and 6.16 respectively. Also calculated 
along with O* and H* temperature profiles is the electron teaperature profile. 
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The upper boundary condition for the ion and electron temperature is that the 
temperature gradient at 1200 km is 1.0 K/k*. The CIRA [lp72] neutral atmos- 
phere model with T m - 1000 K was used. The thermal -electron heating rates 
were taken from Figure 3.3 with x ■ 0*. 

The ion density is important in controlling the magnitude of the differ- 
ence between H* and 0* temperatures. If H + or 0* densities are large 
(-lx 10* cm" 3 ) at 1200 km (Figure 6.15 and 6.16) then !F(H + )/r(0 + ) is re- 
duced over the case where n(0*) ■ 1 x 10 5 cm’ 3 (Figure 6.13), There are still 
significant values of r(H + )/T(0 + ) in the region near 500 km as shown in Figure 
6,15 and 6.16, The magnitude of r(H + )/T(0 + ) in Figure 6,16 is less than that 
in Figure 6,15 because H densities are greater for Figure 6.16, The larger 
H* densities allow more heat to be transferred to 0* by Coulomb collisions. 
Therefore, the H* temperatures are reduced and the 0 + temperatures increase. 

The importance of the electrons as a heat source for ions is shown in 
Figure 6.17. In this figure the electron temperature is much greater at al- 
titudes above 500 km than in Figure 6.16. The electron upper boundary con- 
dition in Figure 6.16 is 32^/a* | l200 ■ 1 K/km and in Figure 6.17 it is l.S K/km. 
If the temperature gradient is larger, the temperature also are larger. 

Since the electron temperatures increase, the H + and 0* temperature also 
increase. The H + temperature increases more than the 0* temperature because 
H* has a larger electron energy transfer rate than 0 + , The difference between 
H* and 0* temperatures increases greatly from Figures 6.16 to 6.17. Therefore, 
as the difference between electron and ion temperatures increases, the dif- 
ference between H* and 0 + temperatures also increases. 

6.7 The Effect of the Neutral Atmosphere cm Ton and Electron Temperature 
Profi lee 

The neutral atmosphere can be important in determining the ion and elec- 


ALTITUDE Own) 



208 


tron temperature profiles. In particular, different temperature profiles can 
result from neutral number densities that differ from one neutral atmosphere 
model to another. The neutral atmosphere is important at lower altitudes 
where the major heat losses for ions and electrons are through elastic and 
inelastic collisions with neutral particles. 

In Figures 6.16, 6.18, 6.19, and 6.20, ion and electron temperature pro- 
files are shown for the cases T m ■ 1000 K, 500 K, 700 K, and 1500 K respec- 
tively, The lower boundary conditions is that the electrons 0 + and H + , and 
neutral temperatures are equal at 120 km. For each neutral temperature, the 
upper boundary condition for electrons, 0 + , and H + temperatures is the gra- 
dient of temperature with respect to altitude is 1 K/km. Ion densities used 
throughout this section are given in Figure 6.6 profile B. The thermal- 
electron heating rates are from Figure 3.3 with x * 0*. The electron tem- 
perature profile is calculated using the method described in Section 6.2.3 
The ion temperature profiles are calculated using the method described in 
Section 6.2.1. 

Comparison between Figures 6.16 and 6.19 indicates that the ion and elec- 
tron temperatures profiles change very little at the altitudes above 800 km 
as the neutral temperature changes from 700 K to 1000 K. This is because 
the dominate heat transfer process above 800 km is thermal conduction. Above 
800 km neutral number densities are very small and cooling of ions and elec- 
trons is accomplished through thermal conduction of heat to lower altitudes. 
Since the thermal conductivity is independent of the neutral atmosphere, the 
ion and electron temperature profiles are independent of the neutral atmos- 
phere above 800 km. 

Below 800 km the effect of the neutral atmosphere begins to determine 
the shape of the ion and electron temperature profiles. For T m • 1000 K, the 
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Figure 6,18 T(0*), r(H + ), and T g where the upper boundary 
condition is 3T^/3s - 1.0 K/km and ion number 
densities are from Figure 6.4 B. The neutral 
atmosphere model is from CIBA [1972] with 
T - 500 K. 
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.19 T(0 ), T(H ), and T g where the upper boundary 
condition is 3T^/3a * 1.0 K/km and ion number 
densities are from Figure 6.4 B. The neutral 
atmosphere model is from CIRA [1972] with 
T • 700 K. 
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0* ion temperature decreases more rapidly than the H* ion temperature, leaving 
a difference between 0* and H* temperatures. As T m decreases, the difference 
between 0* and H* temperatures decreases. As the neutral temperature de- 
creases, the neutral number densities also decrease. As a result, the 0* 
temperature profile and H* temperature profile are much closer to each other 
for smaller neutral temperatures. Thus, the importance of neutral collisions 
with ions and electrons is pushed to lower altitudes for smaller nsutral tem- 
peratures . 

Although the 7(H*) and 7( 0*) do not change at altitudes above 800 km 
when 7" ■ 1000 and 700 K (Figures 6.16 and 6.19), 7(H*) and 7(0*) change when 
7„ * 500 K and 1500 K. In Figure 6.18, (7 w - 500 K) and 6.20 (7^ » 1500 K) 
7(H*) and 7(0*) are much smaller than in Figure 6.19. For 7^ ■ 500 K, 7(0*) 
and 7(H*) are nearly equal from 200 to 1200 km. For 7^ ■ 1500 K, 7(0*) and 
7(H*) are much different especially near 700 km. When « 500 K, 7(0*) and 
7(H*) are so close that it is very difficult to distinguish different tempera- 
tures. Also when T m ■ 500 K, 7(0*) is approximately 1000 K from 400 to 1200 
km. When 7 ■ 1500 K, 7(0*) is nearly equal to T from 200 km to 700 km. 

The electron temperature doesn't change drastically as 7(0*) and 7(H*) 
changes when T m • 500, 700, 1000, or 1500 K. It is approximately the same 
in for each neutral thermosphere model. This is expected because thermal 
conduction is more important for electrons than it is for ions. 

6.8 Sumary and Comparison with Chapter S Results 

In this chapter electron and ion temperatures were found to depend on 
many factors. The electron number density and the 0* and H* number densities 
are important (see Figure 6.3 and 6.5). Also, the neutral thermosphere model 
is an important factor (see Figures 6.16, 6.18, 6.19, and 6.20). The time 
of day can determine 7 below 500 km (see Figures 6.6, 6.7, and 6.8). 
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Inelastic collisions have an effect on the 7 # profile at altitudes below 400 
kn (see Figure 6.10). Thermal conductivity and the upper boundary condition 
are also very important in determining the shape of the T g , 7(0*), and 7(H*) 
profiles (see Figure 6.11, 6.12, 6.13, and 6.17). 

In Figure 6.9 the importance of the various electron heat transfer pro- 
cesses is shown. Thermal conduction is the most important electron energy 
transfer process at altitudes above 400 km. Even though the thermal conduc- 
tivity of 0* and H* is much less than the electron thermal conductivity, the 
thermal conduction of ions is still important in keeping 7(0*) and 7(H*) sig- 
nificantly below 7 . Since the thermal conductivity of 0* is smaller than 
0 

the thermal conductivity of H*, less heat flows through 0* than through H*. 

The electron-H* energy transfer rate is larger than the electron-0* 
energy transfer rate, and 7(0*) and 7(H*) are weakly coupled. Therefore, 

7(H*) is generally larger than 7(0*) (see Figures 6.11, 6.12, 6.13, 6.15, 

6.16, 6.17, 6.18, 6.19, and 6.20). Thermal conduction and the upper boundary 
condition control the magnitude of the difference between 7(H*) and 7(0*) . 

As the ion temperature gradient increases from zero to one K/km the peak dif- 
ference between 7(H*) and 7(0*) . Larger values of 7 result in more heat 
transferring to H* and thus a larger difference between 7(H*) and 7(0*) (see 
Figure 6.17). 

Some of the results of this chapter can be compared with experimental 
results of Chapter 5. The ion densities are given by experimental results 
of Figures 5.3, 5.6, 5.7, and 5.12 and by the empirical model of Figure 6.4. 
Experimental daytime electron temperatures of Figures 5.9, 5.10, and 5.13 
can be compared to theoretical electron temperatures of Figures 6.3, 6.5, 6.15, 
6.16, and 6.17. Nighttime electron temperatures are shown in Figures 5.8 and 
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6.8 the 0* and H* temperatures are given by the experimental results of 
Figures S.9, 5.10, 5.11, 5.15, and 5.14 and by the theoretical results of 
Figures 6.15, 6.16 and 6.17. 

The ion densities of Figures 5.3, 5.6, 5.7, and 5.12 are consistent with 
the ion densitie.. shown in profile A of Figure 6.4 because Figure t.4 is based 
on experimental results. The magnitude and the location of the peak of the 
0* density are approximately the same in both Figures 5.5 and 6.4. At alti- 
tude below 500 km 0* is the dominant ion. At altitudes near 1200 km nCO*)/^ 
is about 0.5. The agreement between Figures 5.5, 5.6, 5.7, and 5.12 and Fig- 
ure 6.4 A is good enough so that the ion and electron temperatures based on 
these figures can be compared. The errors for theoretical ion and electron 
temperatures due to any disagreement in the ion densities would be less than 
the differences between profile 1 and 2 of Figure 6.3 which is less than 100 
K for 7\ and 200 K for T g . In Figure 6.5 the effect of the ion species cm T g 
is shown to be unimportant. 

The electron temperature profiles in Figures 5.9, 5.10 and 5.13 have a 
temperature gradient at the upper boundary of about 1.0 to 1.5 K/km. In Fig- 
ure 6.15 and 6.16 T g has been calculated based on the assumption that the 
temperature gradient at the upper boundary is 1.0 IC/km. Figure 6.17 shows T g 
when the temperature gradient is 1.5 K/km. The values of T in Figures 5.10 
and 5.13 are somewhat greater than those in Figures 6.15 and 6.16 but less 
than those in Figure 6.17. This indicates that the theoretical T g profiles 
are consistent at the upper boundary with the experimental T profile since 
the experimental temperature gradient is between 1.0 and 1.5 K/km at the 
upper boundary. The experimental and theoretical T g profiles don't agree 

below 400 km. The theory predicts 7 to be approximately equal to f(H*) and 
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7(0*) from 300 to 400 ka, and 7 # has a relative maxi m um at about 240 km. How- 
ever, the experimental T § profile becomes equal to 7(H*) and 7(0*) below 230 
ka and T g has a steep teaperature gradient from 300 to 400 ka. This region 
below 400 ka needs to be investigated in future research to find the reason 
for the discrepancy between experiaental and theoretical T g profiles. There 
are several possible reasons for this discrepancy - 1) the heating of the 
aabient electrons by photoelectrons is not accurate below 400 ka, 2) the elec- 
tron density model below 400 ka is not accurate and actual experiaental elec- 
tron densities should be used instead, and 3) there are other sources of heat 
for electron which have not yet been identified as important. 

In Figure 5.8 nighttime electron and ion temperature profiles are shown. 
Although the profiles only extend to 650 km one can still compare then with 
theoretical calculations of T and 7. given in Figures 6.6, 6.7, and 6.8. 
Figures 6.6 and 6.6 don't agree with Figure S.8 because T g is significantly 
greater than TV. from 400 to 650 km in Figures 6.6 and 6.7. However, Figure 
6.8 is similar to Figure 5.8. The T g profile is approximately equal to the 7 ^ 
profile in both Figures 5.8 and 6.8 and the temperatures are fairly constant 
with altitude in both figures due to zero or near zero temperature gradients 
at the upper boundary. 

The agreement between experimental and theoretical profiles of 7(0*) and 
7(H*) is quite good (see Figures 5.9, 5.10, 5.11, 5.13, 5.14, 6.15, 6.16, and 
6.17). Beth the experimental and the theoretical profiles of 7(0*) are close 
to the neutral temperature of 1000 K up to about 500 km where 7(0*) starts to 
increase rapidly. The teaperature gradient at the upper boundary is approxi- 
mately 1.0 K/kra for both experimental and theoretical 7(0*) profiles. 

Comparison between experimental and theoretical 7(H*) profiles cannot be 
made directly because Figure 5.11 and 5.14 give 7(H*)/7(0*) instead. 
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Calculation of values of 7(H*)/T(0*) for Plguras 6.15, 6.16, and 6.17 give 
values of about 1.0 at low and high altitudes which Is the sane for Figures 
5.11 and 5.14. The peak value given by Figures 5.11 and 5.14 Is located at 
550 kn. This is also about the sane location of the peak of 7(H*)/7(0*) 
calculated from Figures 6.15, 6.16, and 6.17. The values of TQf)/T(0*) at 
550 km calculated from Figures 6.15, 6.16, and 6.17 are about 1.3, 1.2, and 
1.5 respectively. The maximum of ?(H*)/7(0*) in Figures 5.11 and 5.14 is 
about 1.3 and 1.4 which are consistent with those of Figures 6.15 and 6.17. 

It is important that the experimental anc Theoretical ion temperature 
profiles agree because it means that all the important energy transfer pro* 
ceases are included in the ion energy balance equations. There is no reason 
to look for additional energy sources or sinks for ions because there is no 
discrepancy between experimental and theoretical ion temperatures. 

The experimental and theoretical electron temperatures agree above 400 
km. However, below 400 km there is a discrepancy which needs to be investi- 
gated by further research of electron energy transfer processes. 
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CIRA 1972 Neutral Atmospheric Models 
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APPENDIX C 

The following is a calculation of the nonanalytic tera in the power 
spectra of incoherent -scatter returns, 
u/n 


Hli/ 

* ft 


exp y dy 


■M, 

' a 


2 

(i ♦ r ♦ Jr * )«fe 


I ♦ + + 

5 


- g- jFjd/2; 3/2; to 2 /n 2 ) 

2 2 

where (1/2 ; 3/2; w /D ) is the degenerate hypergeometric function or Kummei 
function given in Abtvmntita and Stegun [1964], 

Gradshteyn and Ryahik [1965] gives an expression for jFj as 

so ^ „ J 

jFjCl/2; 3/2; w 2 /n 2 ) - exp (w 2 /n 2 ) E C-» 2 /n 2 ) 

2 2 2 2 

This summation is useful when o> /fl is not large. When to /(l becomes 
large the summation will not converge. Therefore, one has to use an asymptotic 
expansion, which is given by [Abramouita and Stegun , 1964] 

2 n 2 2 2 
.^.(1/2; 3/2; u> 2 /fl 2 ) ■ — -y exp (w 2 /^) (1 ♦ -2 *■ + 3 -2y 
1 1 2u> ■>*> 


2b> 

4 


Zb> 


♦ 15 


n 


2 « 


♦ 105 1 » 

2b) 
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APPENDIX D 

COMPUTER PROGRAMS USED IN THE CALCULATION OF EXPERIMENTAL 
AND THEORETICAL ELECTRON AND ION TEMPERATURES 


In this appendix a list and description of some of computer programs 
used in Chapters S and 6 is given. Programs MAC, NOINE, and ACNAL were used 
to calculate the experimental electron and ion temperatures and densities 
given in Chapter 5. The program used to calculate the theoretical electron 
and ion temperatures is also described and given in this appendix. 

The following program calculates a library of theoretical ACF's and 
stores them on tape. Subroutine DISP calculates the electron and ion disper- 
sion functions. Subroutine CFS calculates the cosine factors for the cosine 
transform of the power spectrum. Subroutine ACF calculates the autocorrela- 
tion function. The power spectrum is calculated and stored in array W. The 
variables Ml, M2, and M3 refer to the masses in kg of 0 + , H + , and the elec- 
tron respectively while Tl, T2, and T3 represent these particles temperatures. 
Constant factors necessary for the calculation of the power spectrum are given 
by C, C2, C3, and C4. 

The power spectrum is calculated for a range of values of T(0 + ), 

T /T{ 0 4 ), 7’(H + )/7’(0 + ), and «(0 + )/n (B4). After the ACF (A) is calculated it 

0 c 

is sent to the printer and stored on magnetic tape for use by other programs. 
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PROGRAM NAC< I NPUT , OUTPUT , TAPE? . ACA , TAPES* I NPUT , PUNCH . TAPES* AO A ) 
COMMON/CFA/CFC 24,301) 

DIMENSION H<30n.F<3Bt>,G<301>*A<24> 

DIMENSION H(301 ) 

DIMENSION X<9> 

REAL L,K,H2,M3,P0,K1,PI,M1 


PI=3. 141392654 
K*1.380622E-23 
P0=8.854185E-t2 
Ml=2.655913E-26 
M2=t . 672614E-27 
M3*9 . 109558E-31 
L=. 69719186 
NACF*0 
E=2.0E11 

C=L**2V<32.*PI**2.tK> 

C2=C*PI 

C3 S (1.6021917E-19>«2/PO/K 
C4*<L/4./PI>»2. 

04*1. 0E3 
TS=8 . E-6 

CALL CFS< D4.TS) 

READ< 5, 100 )ITO 
READ( 5, 100 )ITS 
READ< 5, 100 >IT1 
FORMAT* 14) 

T1=ITS*FL0AT< I TO— 1 )*IT1 
REWIND 7 
READf 5, 100 )IF1 
UPITE<6,189)IT0,ITS,IT1,IF1 

FORMAT* IX. 7HN0. INC,I4,7HST TEMP, I4.8HTEMP INC,I4,8HN0 FILES, 14) 

REA0C5, 100 )ITRS 

READ< 5, 100 )ITR1 

REAQ( 5, 100 )IBS 

READ<5, 100)161 

WRITE<6> 110)ITRS, ITRI, IBS , IBI wui . tjlN 

FORMAK 1X5HTR ST,I4,6HTR INC.I4,5HB4 ST,I4,6HB4 INC. 14) 

IF(IF1 .EQ.0)GO TO 85 

00 84 IFIL*1 , IF1 

DO 84 IFK=1 ,96 

READ* 7 >A 

CONTINUE 

CONTINUE 

X1=M1/T1 

CALL DI$P<X1,D4,F,C> 

DO 19 I T = 1 ,4 

T2=T1+FLQAT( IT-1 >*T1/10.0 
X2=M2 T2 

CALL DISP< X2,04,G, C ) 

DO 23 J=1 , 4 

T3=< ITRS-K J-l XITRI >*T1^10 
X3=M3'T3 

CALL DISP<X3,D4,H,C> 

DO 56 18=1,6 

84=<IBS+<IB-i)*IBI V1000 
NACF=NACF+1 

WRITE'’ 6, 104 )Tl , 72, T3, B4.NACF 
FORMATC 1X,3F8.1»F82,I5) 

W1=0 0 

DO ?0 IK=1 , 301 
F1=F'' IK) 

F2=S* IK) 

P 7 s H * IK ) 

F8=' r -84 )*$0PKC2*X2)*EXP<’-C*X2 *W1**2 > 

P5=-F8iWl 

F9 = S0RKC2IX3)*EXP<-C*X3 *W1**2 > 


I 
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I 

! F€«-F9*Wt 

F7*B4 *$QRT<C2*XI>*EXP<-C*X1 tmt*2 ) 

F4a-F7»H1 

38 Al»<l,+C4*C3*E*<B4*Fl/m<i.-B4>*F2/T2>>**2 *F9 

A2«C4*C4*<C3*E*< F4/T1* F3/T2>>**2 *F9 

A3»C4*C4*<C3*E/T3>**2 *<F3**2 +F€**2 >*<F?+F8> 

A4*< 1 >C4*C3*E*< F3/T3+F1*B4/Tt+F2*< 1 . -B4 >/T2 ) >**2 
A5«C4*C4*<C3*E*<F6/T3+F4 /TI+F5 /T2>>**2 

I4< ItO=< AUA2+A3 >/< A4+A5 ) 

UUW1+D4 
70 CONTINUE 

CALL ACF< W. A/D4) 

WRITE<6/103>A 
56 CONTINUE 

23 CONTINUE 

19 CONTINUE 

END 


SUBROUTINE DISP<X1.D4/F/C> 

DIMENSION F< 301 ) 

Z=C*2 
K7=l 
C5=l 0 

70 77 1=1/301 
N1«FL0AT< 1-1 >*D4 
IF(U1 HE 0.0 >G0 TO 82 
Fl*l .0 
GO TO 8b 

IFCK7.LT.6P)G0 TO 83 

F 1 =- 1 X /C5-3 W /C5»*2- 15 . /C5M3-105 . /C5**4-945 . 'C5**5 

GO TO S6 
hK1=1 0 
AF2=1 0 
F1*O.0 
Gl=8 0 

DO 61 K7=l .61/2 

AK1=<2 0 * FLOAT( K7 >-l . 0 >*AK1 

F1=AF2 'AK1+F1 

AF2 S U41*U1*Z*X1)*AF2 

Akl='s2 0 * FLOAT< K7 >+1.0 >*AK1 

F1=-AF2'AK1+F1 

AF2= <UUimZ*Xl>*AF2 

IF'. hB$^F1-G1 VF1VLT. 1 .0E-8)GO TO 69 


>*AK1 


G 1 =F l 


CONTINUE 
GO TO 83 
Fl=l 0 -FI* 

p -: I >=f l 

CONTINUE 

RETURN 

END 


(uituitzm) 


*■ 
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SUBROUTINE CF$< D4,T$> 

COMNQN/CFA/CF< 24,381) 

DO 94 1-1,24 
141*0.0 

CF<I,1 >*i .0000000 
HO*FLOAT<I-l>*T$ 

DO 97 K9*2,298,2 
W1=W1+D4 

CF< I,K9)-4.*C0S<N1*H0> 

M1-N1+D4 

CF<I,K9*1>*2.*C0S<W1*H0> 

CONTINUE 

CF< 1 , 300 >*COS( < I41+D4 )*HO >*4 . 0 
CF( 1 , 381 >*COS< < H1+D4T2 . >*H0 > 

CONTINUE 

RETURN 

END 

SUBROUTINE ACF<W,A,D4> 

COMNON/CFA/CF<24,301 ) 

DIMENSION W< 301 ), A< 24 > 

DO 94 1*1,24 

141*0.0 

H0=FLOAT< 1-1 >*8.E-6 
A( I >*I4< 1 > 

00 97 K9*2,298,2 

A<I )=A< I )+CF< I , K9 >*H< K9 >*CF< I , K9+ 1 >*M< K9+ 1 > 


CONTINUE 

A( I >*A< I >+CF< 1 , 300 >tW< 300 >+CF< I .* 301 >*H< 301 > 
IF( I NE 1 JAC I )-A< I )/A< 1 > 

CONTINUE 

RETURN 

END 
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Below is a list of the physical quantities for the library of ACF's. 

Each file contains 96 ACF's. For files 1.34 7(H*)/?(0*) goes from 1.0 to 1.3 
in steps of .1. For files 35-37 ?(H*)/7(0 4 ) equals 1.4-1. 7 in steps of .1. 


FILE 

nfoVn , 

rco + ) 

T e / T ( 0 *) 

1 

0.0 

m 

0.5 

800 

1.4 

- 1.7 

2 

0.0 

m 

0.5 

1000 

1.4 

- 1.7 

3 

0.0 

- 

0.5 

1200 

1.4 

- 1.7 

4 

0.0 

. 

0.5 

1400 

1.4 

- 1.7 

5 

0.0 

• 

0.5 

1600 

1.4 

- 1.7 

6 

0.0 

• 

0.5 

1800 

1.4 

- 1.7 

7 

0.0 

- 

0.5 

2000 

1.4 

- 1.7 

8 

.55 

. 

.80 

800 

1.4 

- 1.7 

9 

.55 

• 

.80 

1000 

1.4 

- 1.7 

10 

.55 

• 

.80 

2000 

1.4 

- 1.7 

11 

.55 

m 

.80 

2400 

1.4 

- 1.7 

12 

.55 

m 

.80 

2800 

1.4 

- 1.7 

13 

.85 

m 

.975 

800 

1.4 

- 1.7 

14 

.85 

- 

.975 

1000 

1.4 

- 1.7 

15 

.85 

m 

.975 

1200 

1.4 

- 1.7 

16 

.85 

- 

.975 

1400 

1.4 

- 1.7 

17 

.85 

- 

.975 

1600 

1.4 

- 1.7 

18 

.85 

m 

.975 

1800 

1.4 

- 1.7 

19 

.85 

m 

.975 

2000 

1.4 

- 1.7 

20 

.85 

- 

.975 

800 

1.8 

- 2.1 

21 

.85 

- 

.975 

900 

1.8 

- 2.1 

22 

.85 

. 

.975 

1000 

1.8 

- 2.1 

23 

.85 

m 

.975 

1100 

1.8 

- 2.1 

24 

. 8 $ 

• 

.975 

2000 

1.0 

- 1.3 

25 

.85 

- 

.975 

2500 

1.0 

- 1.3 

26 

.85 

- 

.975 

3000 

1.0 

- 1.3 

27 

0.0 

- 

0.5 

800 

1.0 

- 1.3 

28 

0.0 

- 

0.5 

1200 

1.0 

- 1.3 

29 

0.0 

- 

0.5 

1600 

1.0 

- 1.3 

30 

0.0 

• 

0.5 

2000 

1.0 

- 1.3 

31 

.55 

• 

.80 

800 

1.0 

- 1.3 

32 

.55 

m 

.80 

1200 

1.0 

- 1.3 

33 

.55 

- 

.80 

1600 

1.0 

- 1.3 

34 

.55 

m 

.80 

2000 

1.0 

- 1.3 

35 

.85 

- 

.975 

1200 

1.5 

- 1.8 

36 

.85 

- 

.975 

1400 

1.5 

- 1.8 

37 

.85 

- 

.975 

1600 

1.5 

- 1.8 

38 

1.1 

D 

800-1500 

1.4 

- 2.5 
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Program NOINE calculates n(0 + )/n # by fitting theoretical values of TZERO 
to experimental values of TZERO. This pro gran also calculates by fit- 

ting theoretical values of AMIN to experinental values of TZERO. A library 
of theoretical ACF's is read from magnetic tape. Experinental ACF's are 
also read from magnetic tape. Theoretical and experinental values of TZERO 
are found in subroutine AZH. In subroutine FAMIN theoretical and experi- 
mental values of AMIN are found. Subroutine RRAF calculates the array of 
theoretical values of AMIN. Newton's method (see Section S.2.1) is used 
to calculate n(0 + )/« e . This is done in subroutine FINDX. is also 

calculated by Newton's method. The program outputs the results of n(0 + )/« e 
and T g /T i for the 15 altitudes of both short and long pulse data into file 
XFILE for later use by program ACNAL. 

The array AC contains the library of theoretical ACF's. The arrays 
TZT and AMINT contain the theoretical values of TZERO and AMIN. X(l) repre- 
sents n(0*)/n g . X(2) represents 


2 SO 


I14iIAPE»-XFItE, TAPES- INPUT, 


01 HENS I ON ACE<24> 

REMIND 7 
REN I NO 8 
READ<5, 103)IREC 
F0RHATU4) 

IF< IREC.EQ.8XiO TO 231 

00 238 1*1, IREC 

DO 229 >1,15 

READ< 8 )ACF 

READ< 8 )SN 

READ< 8 )EN 

CONTINUE 

00 1 1-1,29 

DO 1 L-1,96 

REAO( 7 >ACF 

DO 2 1-1,24 

READ(7)< AC(N,1,I >,H-1,24> 

URITE(6, 188X AC<N, 1, 1 >,N-1,3> 

REMIND 7 

DO 281 1-1,6 

00 281 L-1,96 

REAO( 7 )ACF 

DO 3 1-1,24 

RE ADC 7 )( AC< N, 2, 1 ), H-l ,24 ) 

URITEC6, 188 )(AC(N,2, t >,H-1,5) 

REM I NO 7 

DO 4 1-1,33 

00 4 L-1,96 

REAO< 7 >ACF 

00 5 1-1,24 

REA0(7XAC(n,3,I),N«l,24> 

MPITE(6,188 )(AC(H,3,24),H>1,9> 

REMIND 7 

DO 6 1-1,9 

00 o L— 1 , 96 

READC 7 )ACF 

DO 7 1*1,24 

READC 7 X ACC H, 4, 1 >,N-1,24> 

MRITEC 6,188 X AC<H,4, 1 >.H-i,5> 

REMIND 7 

DO 8 1-1,23 

DO 8 L-1,96 

READCDACF 

DO 9 1*1,24 

READ' 7 X ACC H, 5 , 1 >, H-l , 24 > 

WRITE' 6 • 180 XAC<H,5,1>,H-1,5> 

REMIND 7 

DO 18 1-1,12 

DO 18 L-196 

PEAD<7)ACF 

O' 1 H 1-1,24 

PEA0< 7 X ACC H, 6, 1 ),H-1 ,24 ) 

WRITE' 6, 108X AC<H,6,1 ),H-1,5> 

REMIND 7 

DO 12 1*1,19 

00 12 L-1,96 

READ' 7 )ACF 

DO 13 1*1,24 

READ' 7 X AC' M, 7, 1 >,H-1.2<*> 

WRITE' 8, 188 X AC'H,7,1>,H-1,5> 

DC 131 1-1,7 

DO 131 J-1,24 

ac' i , : , j >-i 0 

FORMAT' IX, 12F10.5) 
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READ IN EXR ACF 


CALL RRAF<AHIHT,AC> 

C FIND NO+/HE 

DO 78 IFI*1,10 
X< 1 )■ . 075 
X< 2 >*1 .4 
DO 67 JH*2, 16 
READ< 8 >ACE 

HRITE<6# 109)JH,< ACE<H),H*1,5> 

109 FORMAT< I4,2X,5F10.3) 

IFUFI EO. 2. AND, JH.GT.ll )CO TO 66 
101 FORNATC IX, 16F6.2) 

CALL AZH< TZ< JH ># TH( JH ># ACE > 

CALL FAHIN<AHIN<JH),ACE> 

WRITE< 6 , 100 >TZ< JH >, TH( JH >, AHIH< JH ) 

TH< JH )*10.$< TH< JH)-t .3) 

IF(X( 1 ).LT. . 85 >GO TO 31 
DO 20 1*1,6 
DO 19 M«l,24 

19 REX<M)*AC(M,6> 1 ) 

20 CALL AZH<TZT( I ),THA,REX) 

XH1*.975 

XLl* 850 

CALL FINDX< XM1,XL1,X,TZT,TZ< JH>) 

IF(X( 1 ) . CT . XL1 >CO TO 60 
IF< X< 1 ) .LT . XL1 )X< 1 )*XL1 
C FINO NO+/NE IF 33<H0^NE< .00 

DO 21 1*1,6 
DO 22 11*1,24 
22 REX<M)*ACU1,4, 1) 

21 CALL AZH<TZTa>,THA,REX> 

31 IF<X<1> LE. 526)G0 TO 32 

XM1* . 80 
XL1-.55 

CALL FIN0X<XM1,XL1,X,TZT,TZ<JH>> 

IF<X<1).GT. 526)G0 TO 60 
C FINO NO+/NE IF 0.0<NO+/HE< .3 

DO 24 1*1,6 
DO 23 M*l,24 
REX< M >*AC( M,2, 1 ) 

CALL AZH< TZT( I ) . THT< I > , REX > 

THT< I >*'10.t<THT< I >-l .3) 

XL1*0 . 0 
XM1* . 5 

CALL F I NDX< XN 1 , XL 1 , X , THT , TH< JH > ) 

IF(X< 1 ).CT. . 524)X( 1 >*324 
I F< X< 1 ) . LT . 0 , 0 >X< 1 >«0 . 0 
READ IN THE REST OF AH I NT FILE 
FIND TR FROH AMIN 
IF<XU).GE. .83)1X1*3 
IF( X< l ) . LT . . 85.AND.X( 1 ).GT. .526)1X1*2 
IF( X< 1 ) . LE . .526)1X1*1 
IB4=(X< 1 >-XLl )/( XM1-XL1 >/.2*t.49 
IF(IB4.CT 6 >IB4*6 
K2*l 73 

IF<X< 1 >.GT . .83 >X2*2. 15 

DO 65 ITER*i , 10 

IF< X( 2 ) LT . 1 .0 >X( 2 )*1 . 0 

IF<X'2) CT X2)X<2)»X2 

ITR«<X(2)-1 0)*10 *1 .49 

X*1 

IF<’ ITR EO 1 )K*-1 

DX2=^ AMINT*’ 1X1 , ITR, 184 )-AHINT< 1X1, ITR-K, IB4))*K*10. 
Y, ( 2 )*X< 2 AMINT< 1X1 , ITR. IB4)-AHIN< JH))/0X2 
IF' ABS' AMINT' 1X1 , ITR, 164 >-AMIN< JH )) LT .001 )CO TO 66 
€5 CONTINUE 

I?'X'2> LT 1 005^(2 >*1.05 

66 UPITE'6,104>X 

67 WRITER 9, 104 >X 


23 

24 
52 


r 

C 

60 


104 FORMATC F8 . 5 , F8 . 5 ) 

READC 8 >SN 
REA0C8)EN 
WRITEC 6, 101 )SN 
70 CONTINUE 
END 

SUBROUTINE RRAFC AMIN, AC) 

DIMENSION ANINC3,12,6),ACC24,7,24>,REX<24> 

DO 641 K=l,3 
DO 62 1=1,4 
DO 62 J= 1 , 6 
K3=CK-1)*2+1 
J3=J+C 1-1 >t6 
DO 61 11=1,24 

61 REXCM>=ACCM,K3,J3> 

62 CALL FAMINCAMINCK, I, J),REX) 

DO 63 1=5,8 

DO 63 J=1 , 6 
K3=CK-l>*2+2 
J3=J+CI-5>*6 
DO 621 M=1 , 24 
621 REXCM)=ACCN,K3,J3> 

S3 CALL FAMINC AMIN<K, I,J),REX> 

541 CONTINUE 
REWIND 7 
DO 632 1=9, 12 
DO 632 J= 1 , 6 
J3=J+CI-9>*6 
DO 631 M=1 , 24 

631 REXC M )=ACC M , 7 , J3 ) 

632 CALL FAMINC AMINC 3, 1 , J ), REX ) 

RETURN 

END 

SUBROUTINE AZH(TZ,TH,R) 

DIMENSION R( 24 ) 

DO 19 1=2,24 
IFCRCD.LT. .5 )GO TO 20 

19 CONTINUE 

20 TH=( RC I )-( R< I )-R( I- 1 > >*FLOATC I >- . 5 >/< RC I- 1 >-RC I ) ) 
DO 21 11=1,24 

IFCRC II ).LT.0.0)GO TO 22 

21 CONTINUE 

22 T2=FLQAT< 1 1 >+C RC 1 1 > )/< RC 1 1 - 1 >-RC ID) 

TZ=TZ TH 

RETURN 

END 

SUBROUTINE FAMINC AMIN, R) 

DIMENSION RC 24 ) 

DO 1 1=2,23 

IFCRCD.GT.RC I-1DG0 TO 2 

1 CONTINUE 

2 CALL P0L3C A1 , A2, A3, 1 ,R ) 

Xh=-A2''2. /A1 

AMIN=A1*XA**2+A2*XA+A3 

RETURN 

END 

SUBROU T I NE F I NDXC XM 1 , XL 1 , X , TZT , TZ ) 

DIMENSION XC 2 ), TZTC 6 > 

DO 1 ITER=1 , 10 
IFC XC 1 > . LT . XL1 )XC 1 )=XL1 
IFCXC 1 ‘.GT.Xril )XC 1 )=XM1 
IB4=CX< 1)-XL1 )/CXMl-XLi)/.2+1.49 
K= 1 

IFC IB4.EQ 1>K=-1 

[>,:=( TZTC IB4 >-T2TC IB4-K ) )«/C XMl-XLl V.2 
>'3=C TZTC IB4 )-TZ > 

/.' 1 >=;•« 1 >-X3/DX 

IF' ABS'’ X3 > LT . 01 )G0 TO 2 
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CONTINUE 

CONTINUE 

RETURN 

END 

SUBROUTINE POL3<B0,B1,B2, J,R> 

DIMENSION R<24> 

SUBROUTINE TO FIT A THIRD DECREE POLNOHIAL TO ACF 

Jt-J-2 

J2«J-1 

J3=J 

B 1 ■< R< J2 >*FLOAT< J3M2- Ji M 2 HR< J3 >*FLOAT< J 1 M2- J2M2 HR< J 1 XFLOAT 
1< J2M2-J3M2 ) )/< FLOAT< J2*< J3M2- J1M2 HJ3*< JIM2-J2M2 >♦ Jl« J2M2 
1-J3M2))) 

BB«< R< J1 >-R< J3 >-Bl*FLOAT< JI-J3 ) >✓< FLOAT( J1M2-J3M2 ) ) 

B2=R< J1 )-B0*FLOAT< J1M2 >-Bl*FLOAT< Jt ) 

RETURN 

END 


Program ACNAL calculates rCO + ) , T /T(Q *) , and 7’(H + )/r(0 + ) where n(0*)n 
from program NOINE is used. The method used is the method of fitting the 
theoretical ACF to the experimental ACF by minimizing the square of the dif- 
ference between the two ACF’s. 


J* 
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,^2i? A SnS£N8 L iIli£y T i9y TPUT ' XF ^E,TAPE5»IHPUT. TAPES, 

1 1 SEES r,°yi?yi ' ™ P |9 ,XF 1 LE , XF I LE I , TAPE 1 0*XF I LE 1 , ACFF , TAPE7«ACFF > 

COMMON ACC 24, 3, 96) 

2 \ SiHi ISH >>ACF<24) *HTC 16)* HTS< 1 6 ) , ENC 1 6 > * YC 1 > 

g I f ENS 1 ON TZU 6 > * THC 1 6 > * TZTC 4*6), REXC 24 ) , SNC i 6 ) , ACEC 24 ) 

REWIND 7 
REWIND $ 

READC5* 183)IREC 
F0RMATCI4) 

IRC IREC,EQ,8)G0 TO 231 
DO 230 1*1 * IREC 
DO 229 J=l, 13 
REAOC S >ACF 

FQPHATC IX* 12F7.4/1X, 12F7.4) 

continue 

REAOC S >SN 
REAOC S)EN 

F OR M A TC 1 X/ 1 X / 1 X * 1 P8E9 , 1 / 1 X * 1 P8E9 . 1 / 1 X , 8F6 . 2/ 1 X * 8F6 . 2 ) 

CONTINUE 
CONTINUE 
XL4*1.0 
XM4=1 .3 
F0RMATCFS.3) 

READ IN EXP. DATA 
DO 75 IFI=i>2 
DO 72 JH=2, 16 
REWIND 7 
REAOC 8 )ACE 
READ IN XC1)*X<2) 

REAOC 9* 104 )XC 1 )* X< 2) 

WRITEC6,104)XC1 ),XC2) 

UR I TEC 6- 107)(ACECt1),n*l*5) 

FORMAT*' 1X.5F10.5) 

FOPMATC 2F8. 5 ) 

IFC IFI . EQ . 2 , AND . JH . GT . 10)GO TO 72 

IFCXC 1 > , GE . .825)1X1*3 

IFC XC 1 >.GT. .526.AN0.XU >.LT. .825)1X1*2 

IFCXCl.'.LE. .526)1X1*1 

GO TOC 4 1,51, 61 >,IX1 

CONTINUE 

Il=XCl >*10,+i,49 
REWIND 7 
DO 42 1=1-26 
DO 42 L=1 * 96 
READC 7 >ACF 
DO 43 1=1*3 
DO 43 L=i • 96 

RE ADC 7 >c ACCM • I *L )* M=1 .24 ) 

CALL COMP< ACE*UAR* 13, 12* 14* 1 1 ) 

UAR1=UAP 

,v3)=800 +C 13-1 >*400 
NC 2 '=1 0+CI2-1)* 1 
XC4^*1 0+C 14-1 )t . 1 
DO 44 L=i * 96 

READC 7 )C ACC M. 1 > l >* M*1 * 24 ) 

REWIND 7 

REWIND 7 

DO 441 1=2,3 

30 441 1=1.96 

READC 7 *C aCC M, I , L >* M=1 , 24 > 

CALL COMPC ACE* UAR* 13* 12* 14,11 ) 

IFCUAPl LI UAR )G0 TO 451 
UAR1=UAR 

- ' 2 i= 1 4+>' 12-1)4.1 

IF' 13 EO 1 )X<2 ) = 1 CHC 12-1 >*. 1 

X'3>*?08 +f 1 3-2 >*2O0. 

IF- 13 EO 1)X<3>*2000. 

4 =1 0 + ' 14-1 tt 1 


526 ! AND . XC 1 ) . LT . .825)1X1*2 
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451 DO 46 1*1,3 
DO 46 L*t , 96 

46 REA0<?><AC<H,I,L>,H«1,24> 

CALL COMPC ACE, UAR, 13, 12* 14*11 > 
IFC UARl . LT.UAR)GQ TO 461 
UAR1=UAR 

XC4>=l.0+CI4-l>*.l 
XC2>=1 .4+C 12-1 )*. 1 
XC3>=1200.+CI3-1>*.1 
461 DO 47 1=2,3 
DO 47 L*l,96 

47 READ( 7 X ACC M, I , L >» M=l , 24 ) 

CALL COMPC ACE, UAR, 13,12,14,11) 
1 FCUARl.LT. MAR >GO TO 71 

XC 4 >=1 . 0+< 14-1 >* . 1 
XC3>=1600.+C 13-1 >*200. 

XC 2 )=1 . 4+C 12-1 >*200 . 

IF< 13. E0. 3>XC3)*1200. 

48 GO TO 71 

C READ IN ACF FOR 55<N0*/NE< .85 

51 REWIND 7 

1 1 =c XC 1 )- . 55 >*20 .+1.49 
DO 511 1=1,7 
DO 511 L=1 ,96 
511 READC7)ACF 

DO 52 1=1,3 
DO 52 L=i , 96 

52 READC 7 >( ACC M, I , L >, 11=1 , 24 ) 

CALL COMPC ACE*UAR» 13,12,14,11) 
UAR1=UAR 

X(2>=1.4+CI2-l>*.l 
XC3>=600.+< 1 3- i >*200 
IFC 13 . EQ . 3 >X< 3 >=2000 . 

X'4>=1 0+C 14-1 >* 1 


DO 521 L=i ,96 

521 READC 7 X ACC M, I, L >, 11=1, 24) 

DO 531 1=1,18 
DO 531 L=1 , 96 
531 READC 7>ACF 

DO 53 L=l,96 

53 READC 7 >< ACC M • 3 , L ) , M= 1 , 24 > 

CALL COMPC ACE,UAR, 13, 12, 14, 1 1 ) 
IFC UARl . LT . UAR )GO TO 540 
UAR1=UAR 

XC 4 >=1 . 0+C 14-1 >*.l 

XC2>=1 4+< 12-1 )*. 1 

IFC 13 . E0 . 3 >XC 2 >=1 0+C 12-1 >*. 1 

XC3>=2400 +CI3-1 >*400. 

IFC 13 EQ . 3 >XC 3 >=800 . 

540 DO 541 1=1,3 
DO 541 L=l, 96 

54 1 READ' 7 X ACC M , I , L > , M= 1 , 24 > 

CALL COMPC ACE , UAR ,13,12,14,11) 
IFC UARl. LT. UAR >GO TO 71 

XC 4 >=1 0+C 14-1 >*. 1 
XC 3 >=1200 . +C 13-1 >*400. 

■« 2 >= 1 0+C 12-1 >* .l 
55 GO TO 71 

61 REWIND 7 

I1=CXC1>- 35 >*40 +1.49 
IF' II GT 6 >1 1=6 
IFM1 LT 1)11 = 1 
DO 62 1=1-12 
DO 62 L = ! - 96 

62 READ' 7 >ACF 
DO 63 1=1-3 
DO 63 L = 1 - 96 
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*3 REAO< 7 )< AC< M, I , L ), M-l , 24 ) 

CALL COMP< ACE, UAR, 13, 12* 14/ II ) 
UAR1-UAR 

X<2>=1 .4+< 12-1 >*. 1 
X< 4 )*1 , 0+( 14-1 )* , 1 
X(3}=800.+C 13-1 >*200. 

DO 64 1-1,3 
DO 64 L=l,96 

64 PE AD( 7 >< AC< M , I , L > , M* 1 , 24 ) 

CALL COMPC ACEiUAR, 13,12,14,11) 
IF<UARl.LT.VAR)GO TO 641 
UAR1-UAR 

X< 3 )=* ! 400 +( I 3- l >*200 . 

;« 2 >*!.-»+<! 2- 1>*.1 
::;4>=i.0+(i4-i>*.i 
641 DO 65 1*1,3 
DO 65 L=1 .* 96 

65 PEAD(?KACCM.I.L>.H*1.24> 

CALL COMP< ACE, UAR, 13, 12, 14,11) 

X F< MAP 1 ,LT.UAR)GO TO 651 

UhR 1 =UAR 

:<(3>=800 .-K 1 3-2 >*100 
IF( 13 ,EQ . 1 >X< 3 >=2000 . 

Xf 2 >=1 8+( 12-1 )* . 1 
IF( 13 . EQ . 1 )X( 2 >=1 . 4+< 12-1 )*. 1 
X(4>=1 0K 14-1 >< . 1 
651 DO 66 1=1,3 
DO 66 L=1 , 96 

66 READ<?XAC(M,I,L),M=1,24> 

CALL CuMR< aCE.UAR, 13,12,14,11) 

IF' UAR 1 LT OAR >GO TO 661 
:»rtPl=UAP 

y.< 3 >=1000. +( 13-1 XI 00. 

IF' 13 EO , 3 )X< 3 >=2000 . 

X< 2 >=1 . 8+< 12-1 >< . 1 
IFU 3 . EQ . 3 >X< 2 >=1 . 0*< 12-1 >* . 1 
X(4>«1 ,0+< 14-1 X. 1 
66] DO 67 1=2,3 
DO 67 L=1 • 96 

67 READ( 7 X AC< M • I , L ), M=1 ,24) 

CALL COMP( ACE, UAR, 13, 12, 14, II > 

IF< UARl LT l'AR>GO TO 71 

X( 3 >=25O0 . -K 13-2X500. 

IFf 13 EQ 1 W 3 >=1000 
X<2>=1 ,0+( 12-1 X. 1 
IFC 13 EQ 1 >X(2>=1 8MI2-1X.1 
U<4 >=1 0+C 14-1 >* 1 

71 CONTINUE 

MRITE<6. 100>JH,X< 1 > , X< 3 > , XC 2 > , X< 4 > 
1O0 FORMAT' 14. F8 2 , FS . 0, 2F8 . 2 > 

72 WRITE' 10,111 ><X< IXM), IXM=l,4> 

111 F0RMAKF7 4,F7 3F7 0,F?.3> 

READ* S>SN 
READ 1 ' 8 >EN 
WPITE‘1 6, 105 >SN 
105 FORMAT' IX, 16F7.3) 

’’5 CONTINUE 

FWp 

SUBROUTINE CQMP( ACE, UAR, 13, 12,14.11) 
COMMON AC< 24,3,96) 

DIMENSION ACF(24;,ACE'24> 

UAR 1=10. 

DO 51 11=1,3 

00 51 J= 1 , 4 

DO 51 K K = 1 , 4 

L=I l +24*' l'k -1 ) + 6 *' J-l ) 

DO 40 1=1,24 
40 mCF' I )=AC' 1 , 1 1 1 L > 
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sz uar«o . e 

ACF( 1 )*1 ,0 
DO 7 K*2*24 

7 UAR«UAR+<ACE<K>-ACF<IO>**2 

IF<UAR.GT.UAR1 )GO TO SO 
UAR1*UAR 
I2 a J 
13=1 1 
I4=KK 

50 CONTINUE 

51 CONTINUE 
UAR=UAR1 
RETURN 
END 


The following program calculates profiles for T ( 0 )» r(H ), and T @ . 

The triple diagonal mesh method is used to calculate y(0 + ) and T(H + ) . T 

& 

is calculated by integration of the energy balance equation. Arrays D, 

T, and S represent T , T( 0 + ), and T(H + ), respectively. The neutral 
atmosphere model is contained in arrays R, S, and P. Arrays W and Q 
represent the model used for photoelectron heating of thermal electrons. 
The ion densities are given in array N. 


* 








fif'' a 1.-JBW ■ V * ,*!»»!:« <• ■ 


1 '- T t 'h i vo.ES.U* -, O.FS.f-ES.0.FC..C.FC.O.rC.fi.ri:.e.F6.0.ES.0iFs.0 

t Mil P r.,0* FT', 1 . )F "•*. 2 

V) IMM E'f '-iri.u:-0t ]*RSCS.$3M<[ 20 1 ]• Lt 201]. P$£ 6. 21 7 )»Nt 2. IS 1 3. SC 201 J 

40 HIM NT 1 1 0 3 . C C £ 1 7 ] 

50 LORD HrtTR 17. R 
60 60 HUB I860 

70 PRINT "FILE=")17. T N=" . Rt 1 » 93 3 

8S LORD HATH £2*N 

80 FOR 1-17 TO 201 STEP 10 

100 Li = I*5-t 1 15 

110 ..1= 1/2-7. 5 

120 WRITE 1 15, 201 LI. - t 1. 13. PC 2. I 3. PC 3. I ]. PC 4. I]. Pt 5. 1 3. PC 6. 13. NCI . J ]. NC 2. J ] 

136 NEUT [ 

140 LOflC DATA 10. D 

150 N7=H8=1 

160 LORD I'RTfl 18+K0.T 
170 Ljfll: BRTfl 17+K0.S 

ISO T?*SC 201 ]t<2''?> 

190 PRINT " T<C + ) PROFILE " 

200 PRINT L. 1 T ( H+ ’ T<0+) T£ 6 D KI LI CCO+.H+) PE-I' 

2l8 Kl 1 J-lE 1 3*6 

220 F MR [ i =205 TO iii5 STEP 5 

2 90 j“( LI -195 )-'5 

...... I—.Ll 115). '5 

250 T6-CE I ] 

260 T7 = sc J 3 + 1.0. 20571. 236) 

£70 TFI-Tf I 3+ '.'0.285714286) 

£88 GuSUB 1270 

£90 GO BUB 450 

800 L £ .1 3= 1 • B-Lt .1-1 3 ■ 

8 ! 0 I 1 J 3=- I' -K .-H ■* [ 1 1 

:_'j i" ini l: ;o*."ttL. 300 then 350 

0.0 WMTC ' .15 * 8 40 'Ll T ?.T0.T6.2«£:.B.I.C J 3 • L [ J J > C 1 *■ ' T3-T7 ) » S4* < T6-T8) 
i i h M ■ i ' F5 . l'- . f if, 0. _F5. 0. 7F9. 1 
'on in- r , ■ 

. - .* ■_*• i sio 

:■ 0 i :• s ihei -oo 

8T0 *'■ :M' ') 1 T[ 0 - 


411.) .'OF' II 1-1 T F lt:*K'..T 

4 o ,..ii o'- ;02 0 

4 : o hi. 1 9 - t [ £ * j l . ’ ; ro tr. 5E-04 then 160 

4 00 OiOUS 580 
440 PET URN 

450 F [80 0+ HEHT LOSS ONE PRODUCTION 

4 ) ’ T: 01 -. ',§1;6*N£p5. 5*02+2. 8*H5>*03 
400 p-> •: 1 , 1 , Ol*O:*.i.0R< ’8 + T4) )*0 
4sU P 1 - 1 0 . 4*H*C 3 *SQP.( T 4 ) 

490 84-4. 3*03*E^S0R(To+3)*lE+07 

5*10 0 1*3. 3E+10*H3*O3 * 1 6*'. T8 - 16+T7M1 . 5 > 

510 E : - ' ‘:.4*£ + $l *2 + P6*2 + P 1 *£ + *: P6 ‘ 9 T8 + T4 ) >* ( T8-T4 * + C 1*2 ) 
520 B = '.BfCl*'.T£-T7)*3/aS-'16+T7>)*0.52474TSt2.5 

5) 0 6=-2-25*B*lE-07 

5 10 !'■. ■ ; -'.] -F't.i-F 1 )*< To-T4 )+C 1 *< T8-T7)> 

550 I'*0.5£*25*<D+S4*' T8-T6>)*lE-07-TC J+l 3+2*TC J3-TC J-l 3 


5 M 1 " ‘ i-'T TFNFTF H ,.fE RESULTS 
F-'M : 0, r i *1*1. 4 800 < 0 . i -)O0 
54|i m I.. 100.5 00*500*4500 
Goo 0 21 ’.OO* 101*)* 10o. l£00 

6; 6 "o’.'' J OO* 5 00* 5 :tj * 4500 
e. 01 VR'll : 4500, 100, loo, 1£O0 

6)0 I REEL 1> 1*0* 1' 

il-Jo 1 )P L : 1 oo TO 12. ,0 6 i FP 10*3 

650 PLOT )0O*L1*-1 

660 ilF'LOI 0.-0, 8 

620 LUBE. ■ *. 1 .7* 1, 7- 0. 1 'Ll 

680 NFF’T .,1 

6 MO Flip . ■ - i Oij .. IQ 4 .j uO 6 1 1 F 500 
7T*0 F; OT TO* 58* -1 
~10 F FLO) -2.5*0 


ORIGINAL T*. 
™ °oor orr-itry 


•"It rtf F R - n ire K. 

L *p t .i . - 1 

<3+ 1 .7.P; l 1 

ro s? 

i K 

■ i . * 1. 1 * ~ 


;■ ! it 

. i • «i 


HEIGHT KM" 
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plot tt.u.-l 

NEXT 1 
FEN 

FOP 1=1 TO EOT 
l1*I*5+195 
T' = ST 1 It I 


0:0 To 1 2 u S ! EP 


UT IT 11,1.1' 


FEN 

RETURN 

PEN HERE 7 ' H + > z R nF !LE IS CHLCULRTEIi 

lord date i?»ko.t 

LORD DATE 18+Ku-S 
PRINT " T<H+'‘ PRuPILE" 

PRINT " LI TCHfi TCO+> TE B 
PC 1 ’ LI ! 7-0 

FOR „1=lu: 10 1100 SI EE’ 5 
J= i L 1 - 1 95 T -'5 
: - > L 1 - 1 1 5 1 ■ 5 
■ -o J 

>: . . ?••.' ... to 
• : ! T 1 . . U 7 P * . o i ■ 


LI C( 0 +tH+) pe-: 


; c l: 1 l . 

! :.i • i > • 1' U )■ 1 ’ : ij 

11 ■ * 1 -'‘'O'# . Too T -iE.n !.:00 

1 ! :■ I : 1' • ■■■•< 1. • ■; • l..-:- D-l [ t J. LI J M 

loo u •. i i 

I .10 .0 :• • nil 

. i : 'HEN 1. 'j- 

1 O . O i f L ' : 

! . *10 ,■ T : 0 0 

1 „ ‘ 0 ; ■ . -U h 7 


i .■ V 1 : *. 

1 ■ 1 ■ . I 1 


•1 :~tll 140 


.'Mi HO.' 

: F J _# I Hi ■ J i > IHtll 1440 

= 10 T . LCT' NE 3 J si o- i . M i ■ 

H 7 - 1 0 t 1 LCT' NT 1 > S' j ■ + L G T o N [ 1 1 J 2 + 1 T> 
OOTO 1460 
o l. i i 
H - Ml 1 • J .. J 


A, 


bJc 

%% 

'V 


. i ! ; ! ' 1 ’ ' l 2 '■>: ] • * T [ ... UU H 20 1 ] + * L 2 UU 5 *N 7 *T[ 200 3 t ■ 
i • . ii! 1* L . *i ! J 
r • 7 ; F •! 

: * : : : i : ■ i* i ; ♦ : : 

. i . l : 4 * i_i i j 

l l I ] III M ‘ n: N i‘ ' 

■ J 

'« c . ‘ ■ I 

' UK 1 — i i TO 200 i EP So 

_ 1 - 1 * 5+ 1 9S 

- 2 - 1 * 5 + 24 * 

1 ♦ j+ 29 . 

' 2 - Ti i + io’it’ 2-7 


0 ’ L 1 ■ L i I i ■ 1 ! i • U ‘ L 2 » L L 1 * 1 0 J ft I + l u J ■ f 2 » L J * T 4 * T 3 » L 4 » 7 5 « TU 
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1740 PEM N* HEFT LOSi HNo PRODUCTION 
1750 52*' J. 1 +N2 +2. 8*08+5. 5*H5 )*H3 
1 780 F*P= ’ 0 . 34+H3+SQR c T7 ) > *0 
1770 R t = > 1 .4*H«H3'*SClRa?+T4> 

1780 S4*77*H3+E . SQRC T S+3 9* 1E+87 

1790 C lO. 3E + 1C*H3*07 <l€*<Te •■16+TD + 1.5) 

1800 B = ‘;$4*2+S2*2+P6*2+R1 *2+<P6 - T7+R1 ='<T?+T4> >#CT?-T4 >+Cl*2> 

1810 R*»,B*C 1»iT?*T8>»-3 •(T8/164T7»>*0. 13/7.'T?42.3 

1870 E: = -7-25+B*lE-e? 

1830 5 = - •82+P8tRl )♦< 1 7-T4'4Cl*<T7-T8)J 

1840 D=0. 1 7*25* ( D+S4‘» T7-T6T >*lE-07-TC J+l 3+2*TC J 3-TC J-l ] 

1 7.0 PC TURK 

1880 01*0.00888 

1870 32*0.01955 

1880 D -' = 0,0279S 

1098 Kl*8.62E-85 

1900 LUfiO DHTF 1 • 14 

1910 rOR 1*1 TC 178 

1980 IF I 2-INK I='2> THEN 1950 

1930 OC 1 KHt I 8+0.51 J 

1940 bOTO 1980 

1950 Qt I I -21+WC I - 2+1 ]K2 

1 9(30 NEXT I 

1970 FOP 1=179 TO 21 T - 

1980 at I 1 = QC 178 J 

1990 NEXT I 

8000 FOP 1=1 TC 89 

80 10 FOR j-l TC f 

-020 F[ J. 1*2-1 >Pt J> I 1 

2080 NEXT J 

.:imh ttr 1 i 

• oP ,-2 1C 178 . : L'P 2 

. •' • :F 1 TC t. 

,•..‘-0 IP i • then 210. 

•' ■ . J : o r ■. t.GT ; i J ■ I - 1 1 ■ + LGT ' Pt J* 1 ♦ 1 1 ' ' 2 ' 

8 .i -M 1 ■ ! 71 !0 

2100 ' T i. i !•-■ p[ j. 1 -1 !•+'[ J. 1 + 1 1 > ■ 2 
7: io ns.:-: * j 

2 1 80 ir’ ' I 
iji/y FuP 4=89 TO 98 
2140 FOR 9=1 TC t. 

2150 PC )• * 1 0-7 1 3 KRl J. n 
2180 he;-: I J 
2)70 NEXT I 

8180 *■ OF. ; ■■ . 7 7 Tu 70; - TEP 10 

7 = 0 'P :!-■•! TO !*!■ 

2.- 00 1 v'R ‘=1 1C t 
. . !■> 7 - : hFN 2240 

7. VI I- . Ml. ti ' LG T ■ PI J. I) ‘-LCT-Pt J. 1 + 101"*' I - 1 1 ' 10 + LGT < Pt J» 1 1) > 

- ■ 1 ' . . ■ s o 

.. .0 , ■ ; : i= Ft I. 13 f 1 J. l + lo 1 •*■. I-1 1 ■ 70+PC Jf I) 
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2. - 0 i‘ 7 + II 
7.70 of::: I 
.7:0 Ft I '.‘F 0 

...--I i -jP ‘O-i TO 1.5 ;• t £F 0.5 
.7 00 = .: • • PI .17 1 

8 ; i o . i + 8>.i0 . n r 3 
.11 i, r s T **-TC 2C 1 1 + 1 2 :• • 

7 r-’ 21 ~ j=D 8 

2.40 0 ■ L • I 01 =0 

.' ' 10 P ■ = ; t 1 ’.- ■ 5 ■•NO 

. > .• iM'l 15.8: 90' ..-kO. Tt 201 )»■ 2 7 08 

. ' i ‘ • V.F ,77 1 -5 TO ! 0 4 TEP •? 
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. -4 »--Ji 
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2bl0 FOR' 1 TC 26 
;..\V! T > • A I : ‘ . 

£bJU GO SUB 3220 

2s40 F 

2650 FOR IH«1 TO 3 

2660 T6=1'5*0.<HM1»50*6 

2670 G0SU6 322E 

£660 IF Ft2<Gt</3 THEN 2710 

2630 IF U >Z flNI F *6<0 THEN £780 

£700 NEXT Ml 

2710 IF F ; 0 THEN 2740 

2720 X?*T6 

2730 GOTO 2rS0 

£74g ;-':J=Tfc 

£750 D$=B5t- X8-X7 >*G nBSG £ 

£760 IF X3-X7-:: FIND F»G<0 THEN 2730 

£770 NEXT N 
£".;0 tl 1 i E‘, 

27=41 GOTO 2870 
2B00 IF D5>DS+!00 THEN 2920 
2610 IF D5.-T7-10 THEN £860 
£620 1 1,* ■ lit J + l .+T7>/£ 

£630 GO SUB 3220 

2540 IF C4:0Ei:*E-F then 2920 

£850 GOTO 2580 

£ 8 o ij s = U 9 = U 8 + F * 5 - • 2 

£6 '0 IF INK fLIV 300;#<Ll:' '308 THEN 2900 

2380 WRITE • l5>2&90)Ll.T7.Iim.C4.F.Fl+F2.F3+F4,Q[ I KE,L9.K9tM 
9Rsn FORMAT F5.0. F6 . 0 . F5 . 0 > E9 . 1.4E9. 1 . F7. 1 .F7. 2.F4.0 

£900 NEXT Ll 
2910 GOTO 9000 

£4£ft WRITE” 15»2&90>I 1 ■ T7> T6> W1 *3. 22'4T6T£> F • F 1 +F£i F3+F4 ■ OC I J«E»L9.K9 

£990 IF D5--T7 THEN 2?60 

£940 1:7 = 08 
£350 GOTO 7’ 9 76 
£1F-J £8=119 


£9:8 j IF N5E.-D6-D7' 1- - 02 THEN 3020 

£940 ":|i 9 9 90 

3000 OF I TE •. 15« 3C'10> .1 >Tr*T0.T6>M 
3010 FORMAT F5.0! 3F6, 0.F4.0 
302O 3T0PE DfiTfi 19+i n.E 
.90 30 ; 9 S I ! 9 ] 60 

3040 ! F * HE 3 ( DC £17 1-29 29-, IE-03 THEN 3060 

3050 GOTO 2300 

3060 ( 4=15=1 

3070 NEXT NS 

3080 ST OF 

3090 PEli CALCULATE EMERGV LOSS AND PRODUCTION TERMS 
3 100 $1=4. 8 + 0 3 -- 1 m 00 


3iio 33-rr+H‘j loo© 

3U'0 S>Kr?E“l5*N2 
3 I 30 34- J , 2 1 E - 14*02 
’■ 1 4 *. . \ -4r-n*o 
■ /..r r - E - 12 - H 
. >0 * E-13+H5 

;:r:i • ; -10+C2 

io*t 2 

3 1 ■'■'H i : ! r 1 j *\_ 4 

3 200 S 2 - 3 . 4 3 E - e ?*0 * E 1 5 + 3 *E XF < - D 1 / K J T4>+EXP<- V 3 •- \ . I ' T 4 > > 
3210 Ft TURN 
3220 31=1 

T '■ -i«Oi THEN .:25U 
3..40 ='!-■■! -t/’P* -■ ?E + 04 ’ » 1 2' 

■■/50 =11 ■ 07*t 2 4 E*E F ■ 2 . * IE- 03*- T6-2T.OO • »* 0 \ 

•>o ■ o , e :p - 1-2 * i t*>-e::P‘ * r <2 » 1 T4-- 

\.:\1 '4 .n, ■» fi : * 1 E IP - S3 M 7* •«£:•' F •- P> M T4 *-» , -'4 
. • o ifTi. V * Ii I * < ' K i - D 1 1 1 T 6 * - E l IF' • - D 1 \ i I 4 * 1 

- * i - ♦ • ■ F ■ T t . ■ * V i 

: 0.1 ■ * ■ ’ • : . T-E“ 02* . FT* ■ 


2 "1.- - jf.E-03*1Y. 

3 , *■ 0. 4 K 5 TH£ = 2340 

; : . : c . i o4 * T 

• c -o :• o. then ro 

.' in ■ ; i . 

" • ’ : : I - n ’ F 7 - 

• «. o : ; 1 411 * u. * J . * ' ♦ * 0 ♦ I- . 4 * H*: * *• . c- 

-0 * 4 1 - * • * 1 E- 12 i 

4ir. 48! :c 4 l . 43* 1-R- TOt*. . ■ 1 + 3 . 22* 7 t » 2 * U 1 ■ 

■4!.i ■ i * - 1 ■* . ♦ ;•.»/< to *b suR'Tti::* 

-i. .* * - ! »' ! «4 - * t * ' o- : r'. *• - . "<*( * 8 ijR T f*. ♦ 4 ♦ Q 3 + 8?' * i 1 * • V ' 

•i i •* « F . • • T p. - ’ -J ■ *E ; PF r t- 

• : * f * ■ I »' - ' 4 ♦ ■ i ’R T t * V 

I'}* F 1 f-c- F 4 


i 4n 


242 




■ rti i. o. *> i. r.«. i joo 

V4 m mx i •. ii:'0.0.05« i 1 1 . 

4 'v r. i i • uo-xoo. : joo 

U4i .h..:. l.i.'t.n.o?- ; • 1.?? 

i * t . > H..!; !. ?:• 100. . 00. 1200 

- ■'•. l' h&‘ i. • * • t • 1 • 0. 1 1 

- *f :--.-oc to i.ini ctep ioo 

.i‘4n i ,.o: i.li.-i 
;*.? o .plot -4.;. -e.3 
> .i iT i! '433 THEM 
TO I LOT -0.5.0 
35S0 LHBEL «,*> 1.5. t. 5.0. t)LI 
3590 NEXT LI 
3000 i’OP 1=1 TC 201 

3610 f-LO T t[ I ]T t 3- 7> C L 1 1> > 2 T ) i I *5+ 1?5> -2 
3620 NEXT I 
3630 ‘TOP 


| 

i 

5 

j 

| 





243 


REFERENCES 

Abramowitz, M. and I. A. Stegun [1964], Handbook of Mathematical Function* 

U. S. Government Printing Office, Washington, D.C. 

Banks, P. M. (1966a), Collision frequencies and energy transfer for electrons. 
Planet. Space Soi . 14, 1085-1102. 

Banks, P. M. (1966b), Collision frequencies and energy tranfer for ions. 

Planet. Space Soi. 14, 1105-1122. 

Banks, P. M. (1966c), Charged particle temperature and electron thermal 
conductivity, Am. Geophye . 22, 577-587. 

Banks, P. M. (1967a), Temperature coupling of ions in the ionosphere, Planet. 
Space Sai. IS, 77-93. 

Banks, P. M. (1967b), Ion temperatures in the upper atmosphere, J. Geophye. 

Ree. 72, 3365-3385. 

Banks, P. M. and G. Kockarts [1973], Aevonamy : Port B, Academic Press; New 
York and London. 

Banks, P. M. , A. F. Nagy, and W. I. Axford [1971], Dynamical Behavior of 
thermal protons in the mid-latitude ionosphere and magnetosphere, 

Planet. Space Soi. 19, 1053-1058. 

Bowhill, S. A. [1962], The formation of the daytime peak of the ionospheric 
?2 layer, J. Atmoe. Terr. Phye 24, 503-519. 

Bowles, K. L. [1958], Observations of vertical incidence scatter from the 
ionosphere at 41 Hc/sec, Phye. Rev. Lett., Vol. 1., So. 12 H623 1/2- 
H523 2/2. 

Bowles, K. L. [1961], Incoherent scattering by free electrons as a technique 
for studying the ionosphere and exosphere: Some observations and theoret- 

ical considerations, J. Ree., RBS, 6SD 1-14. 


244 


Bowles, K. L. [1962], Profiles of electron density over magnetic equator ob- 
tained using the incoherent -scatter technique, MBS Teeth. Mote 199, 1-12 
National Bureau of Standards, Boulder, Colorado and the Instituto 
GeofLico del Peru, LIMA. 

Brachman, R. T., W. L. Fite, and R. H. Neynaber [1958], Collisions of electrons 
with hydrogen atoms III, Phye. Rev. 119, 1157-1161. 

Breigg, E. L. and C. C. Lin [1966], Excitation of the spin sailtiplets of the 
ground state of oxygen by slow electrons, Phye. Rev. 161, 67-79. 

Buneman, 0. [1962], Scattering of radiation by the fluctuations in a non- 
equilibrium plasma, J. Geophye. Res. 67, 2050-2053. 

Cauffman, D. P. and D. A. Gumett [1971], Double-probe measurements of 
convection electric fields with the Injira-5 Satellite, J. Geophye. 

Fee. 76, 6014-6027. 

Chapman, S. and T. G. Bowling [1970], Mathematical Theory of Monmifom 
Gases, Cambridge Univ. Press, London and New York. 

CIRA [1972], Co8par International Reference Atmoephere. 

Cohen, R. S., L. Spitzer, and R. Routley [1950], The electrical conductivity 
of an ionized gas, Phye. Rev. 80, 230-238. 

Cooper, J. W. and J. B. Martin [1962], Electron photodetachment from ions 
and elastic collision cross sections for 0, C, Cl, and F., Phye. Rev. 

126, 1482-1488. 

Dalgamo. A. and T. P. Degges [1968], Electron cooling in the upper atmos- 
phere, Planet. Space Sai. 15, 125-127. 

Dalgamo, A. and R. J. W. Henry [1965], Electron temperatures in the D 
region, Proc. Roy. Soc. A 286 , 521-530. 

Dalgamo, A. and R. J, Moffett [1962], Electron cooling in the V region, 


Planet. Space Sai. 8, 439-441. 


245 


Dalgamo, A., M. B. Me Elroy, and R. J. Moffett [1963] , Electron temperatures 
in the ionosphere, Planet. Spaa « Sot, 11, 463-484. 

Dalgamo, A., M. B. NcElroy, and J. C. Walker [1967], The diurnal variation 
of ionospheric temperatures Plant t, Spaat Sat, IS, 331-345. 

Dalgamo, A., M. B. McElroy, M. H. Rees, and J. C. G. Walker [1968], The 
effect of oxygen cooling on ionospheric electron tesperatures , Plantt. 
Spaat Sat. 16, 1371-1380. 

Da Rosa, A. V. [1966], The theoretical time dependent thermal behavior of 
the ionospheric electron gas, J. Geophye . Res. 71, 4107-4120. 

Desloge, E. A. [1962], Exchange of energy between gases at different temper- 
atures, Phye. Fluids S, 1223-1225. 

Dougherty, J. P. and D. T. Farley [1960], A theory of incoherent scattering 
of radio waves by a plasma. Proa. Ray. Soe. A 2S9, 79-99. 

Drukarev, G. [1946], (hi the mean energy of electrons released in the ionization 
of gas, J. Phys. 10, 81-84. 

Englehardt, A. G., A. V. Phelps, and C. G. Rish [1964], Determination of 
momentum transfer and inelastic collision cross sections for electrons 
in nitrogen fusing transport coefficients, Phys. Rev. j 35 , 1566-1574. 

Evans, J. V. [1962], Diurnal variation of the temperature of the F region, 

J. Geophys. Res. 67, 4914-4920. 

Evans, J. V. [1969], Theory and practice of ionospheric stv. y by Thomson 
scatter radar, Pvoo. IEEE 57, 496-530. 

Fejer, J. A. [I960], Scattering of radio waves by an ionized gas in thermal 
equilibrium, Can. J. Phye. 38, 1114-1133. 

Fejer, J. A. [1961], Scattering of radiowaves by an ionized gas in thermal 
equilibrium in the presence of a uniform magnetic field. Can. J. Phys. 

39, 716-740. 


246 


Geisler, J.* E. [1967], On the limiting daytime flux of ionization into 
the protonosphere , J. Geophys. Res. 72, 81-84. 

Geisler, J. E. and S. A. Bowhill [1965], Exchange of energy between the 

ionosphere and the protonosphere, J. Atmos. Terr. Fhys. 27, 1119-1146. 

Gleeson L. J. and W. I. Axford [1967], Electron and ion temperature varia- 
tions in temperate zone sporadic -E layers. Planet. Space Sci. IS, 123-136. 

Gordon, W. E. [1958], Incoherent scattering of radio waves by free electrons 
with applications to space exploration by radar, Proc. IRE 46, 1824-1829. 

Gordon, W. E. [1964], Arecibo ionospheric observatory. Science, 146 , 26-30. 

Gradshteyn, I. S. and I. M. Ryshik [1965], Tables of Integrals, Series, and 
Products, Academic Press, New York. 

Hagen, J. B. and D. T. Farley [1973], Digital correlation techniques in 
radio science. Radio Sci. 8, 775-784. 

Hagen, J. B. and P. Y. Hsu [1974], The structure of the protonosphere above 
Arecibo, J. Geophys. Res. 79, 4269-4275. 

Hagfors, T. [1961], Density fluctuations in a plasma in a magnetic field with 
applications to the ionosphere, J. Geophys. Res. 66, 1699-1712. 

Hake, R. D. and A. V. Phelps [1967], Momentum transfer and inelastic cross 
sections for electrons in O 2 , CO, and (X^, Phys. Rev. 158, 70-84. 

Hall, L. A., K. R. Damon, and H. E. Hinteregger [1963], Solar extreme ultra- 
violet photonflux measurements in the upper atmosphere of August, 1961, 
Space Res. Ill, 745-771. 

Hall, L. A. and H. E. Hinteregger [1970], Solar radiation in the extreme 
ultraviolet and its variation with solar rotation, J. Geophys. Res. 

75, 105-111. 


247 


Hall, L. A., W. Schweiger, and H. E. Hinteregger [1965], Improved extreme 

ultraviolet absorption measurements in the upper atmosphere, J. Geophya . 
Fes 70 , 105-111. 

Hall, L. A., J. E. Higgins, C. W. Chagnon, and H. E. Hlnteregger [1969], 

Solar cycle variation of the extreme ultraviolet radiation, *T. Geophya . 
Rea. 74 , 4181-4183. 

Hanson, W. B. [1963], Electron temperatures in the upper atmosphere. Space 
Rea. Ill , 282-302. 

Hanson, W. B. and R. Cohen [1968], Photoelectron heating efficiency in the 
ionosphere, J. Geophya. Rea. 73, 831-840. 

Hanson, W. B. and I. B. Ortenburger [1961], The coupling between the protono- 
sphere and the normal F region, J. Geophys. Rea 66, 1425-1435. 

Hanson, W. B., S. Sanatani, L. H. Brace, and J. A. Findlay [1969], Thermal 
structure of an Alouette II topside profile as deduced from rocket 
measurements , J. Geophys. Res. 74, 2229-2239. 

Hasted, J. B. [1964], Phyaica of Atomic Collisions Butterworths , London. 

Hinteregger, H. E. [1965], Absolute intensity measurements in the extreme 
ultraviolet spectrum of solar radiation, Space Sci. Rev. 4, 462-497. 

Hinteregger, H. E. [1970], The extreme ultraviolet solar spectrum and its 
variation during the solar cycle, Ann. Geophya. 216, 547-554. 

Hoffman, J. H. [1967], A mass spectrometric determination of the composition 
of the nighttime topside ionosphere, J. Geophya. Res. 72, 1883-1888. 

Hoffman, J. H. , W. H. Dodson, C. R. Lippincott, and H. D. Hammack [1974], 
Initial ion composition results from the Isis 2 satellite, J. Geophya. 
Rea. 79, 4246-4251. 

Huxley, L. C. H. and R. W. Crompton [1962], The motions of slow electrons 
in gases, in Atomic and Molecular Processes, Academic Press, New York. 


248 


Knudsen, W. C. and G. W. Sharp [1966], Ion temperature profile in the top* 
side ionosphere, J. Geophys. Rea. 71, 4099-4105. 

Lane, N. F. and A. Dalgamo [1969], Electron cooling by vibrational excitation 
of O 2 , J. Geophys. Rea . 74, 3011-3012. 

LeDoumeuf, M. and R. K. Nesbet [1976], Electron-impact excitation of fine- 
structure levels of atomic oxygen, J. Fhye. B . 9, L241-L243. 

Mantas, G. P. [1973], Electron collision processes in the ionosphere, Aevon . 
Rep. No. 54, Aeron. Lab., Dep. Elec. Eng., Univ. 111., Urbana-Champaign. 

McClure, J. P., W. B. Hanson, A. F. Nagy, R. J. Cicerone, L. H. Brace, 

M. Baron, P. Bauer, H. C. Carlson, J. V. Evans, G. N. Taylor, R. F. 
Woodman [1973], Comparison of T and T. from Ogo 6 and from various 
incoherent scatter radars, J. Geophys. Res. 78, 197-205. 

McDaniel, E. W. [1964], Collision Phenomena in Ionized Gases, Wiley, New York. 

Moorcroft, D. R. [1963], On the power scattered from density fluctuations 
in a plasma, J. Geophys. Res. 68, 4870-4872. 

Moorcroft, D. R. [1964], Determination of temperature and ionic coupling 

by incoherent scatter from the ionosphere and magnetosphere, J. Geophys. 
Res. 69, 955-970. 

Nagy, A. F., P. Bauer, and E. G. Fontheim [1968], The nighttime cooling of 
the protonosphere, J. Geophys. Res. 73, 6259-6274. 

Nagy, A. F., L. H. Brace, G. R. Caugnan, and M. Kanal [1963], Direct measure- 
ments bearing on the extent of thermal nonequilibrium in the ionosphere, 
J. Geophys. Res. 68, 6401-6412. 

Nagy, A. F., E. G. Fontheim, R. S. Stolarski, A. E. Butler [1969], Ionospheric 
electron temperature calculations including protonospheric and conjugate 
effects, J. Geophys. Res. 74, 4667-4676. 


249 


Nicolet, M. [1953], The collision frequency of electrons in the ionosphere, 

J. Atmos . Terras. Fhys. 3, 200-205. 

Pack, J. L. and A. N. Phelps [1961], Drift velocities of slow electrons in 
helium, neon, argon, hydrogen, and nitrogen, Fhys. Rev. 221, 798-806. 

Pineo, V. C., L. G. Kraft, and H. W. Briscoe [1960], Some characteristics 
of ionospheric backscatter observed at 440 Mc/s, J. Geophys. Res. 65, 
2629-2633. 

Rees, M. H., R. A. Jones, and J. C. G. Walker [1971], The influence of 

field-aligned currents on auroral electron temperatures. Planet. Space 
Sot. 19 , 313-325. 

Rees, M. H. , J. C. G. Walker, and A. Dal gam o [1967], Auroral excitation of 
the forbidden lines of atomic oxygen. Planet. Space Sd. IS , 1097-1110. 

Renau, J. [1960], Scattering of electromagnetic waves from a nondegenerate 
ionized gas, J. Geophys. Res. 66, 3631-3640. 

Rishbeth, H. and D. W. Barron [1960], Equilibrium electron distributions in 
the ionospheric F 2 -layer, J. Atmos. Terr. Fhys. 18, 234-252. 

Rohrbaugh, R. P. and J. S. Nisbet [1973], The effect of energetic oxygen 
atoms on neutral density models, J. Geophys. Res. 78, 6768-6771. 

Rycroft, M. J. and S. J. Burnell [1970], Statistical analysis of movements 
of the ionospheric trough and the plasmapause, J. Geophys. Res. 75, 
5600-5604. 

Salpeter, E. E. [1960a], Electron density fluctuations in a plasma, Fhys. 

Rev. 120, 1528-1535. 

Salpeter, E. E. [1960b], Scattering of radio waves by electrons above the 
ionosphere, J. Geophys. Res. 65, 1851-1852. 

Sanatani, S. and W. B. Hanson [1970], Plasma temperature in the magnetosphere. 


J. Geophys. Res. 75, 769-774. 


250 


Schulz, G. J. and J. T. Dowell [1962], Excitation of vibrational and electronic 
levels in O 2 by electron impact, Phys . Rev 128 , 174-177. 

Schunk, R. W. and P. B. Hays [1971], Photoelectron energy losses to thermal 
electrons, Planet . Space Sci. 19, 113-117. 

Smith, P. A. [1968], The determination of vertical distributions of plasma 

temperatures and composition from satellite measurements, J. Atmos. Terr. 
Phys. 30, 1203-1209. 

Smith, K. , R. P. McEachron, and P. A. Fraser [1962], Effect of virtual 

excitation of the 2s state on the elastic scattering of electrons by 
atomic hydrogen, Fhye. Rev. 125, 553-557. 

Spencer, N. W. , L. H. Brace, G. R. Carignan, D. R. Taeusch, and H. Niemann 
[1965], Electron and molecular nitrogen temperature and density in the 
thermosphere, J. Geophye. Res. 70, 2665-2698. 

Spitzer, L. [1962], Physics of Fully Ionized Gases, 144, Interscience 
Publishers. 

Stubbe, P. and S. Chandra [1970], The effect of electric fields on the F- 
region behavior as compared with neutral wind effects, J. Atmos. Terr. 
Phys. 32, 1909-1919. 

Swartz, W. E. and J. S. Nisbet [1972], Revised calculations of F region 
ambient electron heating by photoelectrons, J. Geophys. Res. 77, 

6259-6261. 

Tambc, B. R. and R. J. W. Henry [1976], Low-energy scattering of electrons 
by atomic oxygen, Phys. Rev. A 13, 224-228. 

Walker, J. C. G. [1968], Electron and nitrogen vibrational temperature in 
the j?- region of the ionosphere. Planet. Space Sci. 16, 321-327. 


251 


Walker, J. C. G. [1971], Ionospheric effects of auroral electric field, in 
The Radiating Atmosphere, ed. B. M. McCormack, Reidel Publ . Co., 

Do rdrech t - Ho 1 1 and . 




